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With  the  increasing  development  of  our  coastal  areas,  it  becomes 
necessary  to  understand  the  long  and  short  term  behavior  of  the 
dynamic  beach-dune  system  in  order  to  arrive  at  sound  engineering 
decisions  regarding  future  use  of  these  regions,  and  to  protect  exist- 
ing development  threatened  by  beach  erosion.  One  aspect  which  would 
greatly  assist  these  efforts  is  the  ability  to  accurately  predict  the 
amount  of  dune  recession  which  would  occur  during  extreme  storm  events. 

This  research  has  been  directed  toward  a  semi -empirical  solution 
through  the  use  of  a  small-scale  movable-bed  physical  model.  Based  on 
the  principles  of  hydraulic  similitude,  a  new  model  law  has  been  derived 
and  verified.  The  verification  involved  the  reproduction  of  an  actual 
storm  event  to  a  degree  of  realism  never  before  attempted,  and  the 
resulting  erosion  in  the  model  almost  exactly  reproduced  that  which 
occurred  in  the  prototype. 

In  order  to  determine  the  effects  of  wave  period,  wave  height, 
surge  level,  surge  duration,  and  dune  height,  these  parameters  were 
systematically  varied  in  41  experimental  model  tests.  Beginning  with 
an  equilibrium  profile,  the  surge  level  was  increased  over  a  typical 
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time  span  and  then  maintained  at  the  peak  level  until  near-equilibrium 
conditions  prevailed.  Qualitative  comparisons  between  the  resulting 
profiles  indicated  that  it  is  necessary  to  employ  a  time-dependent 
surge  level  increase  in  order  to  more  accurately  simulate  an  actual 
storm  event.  In  addition,  it  was  found  that  each  set  of  conditions 
appeared  to  form  an  equilibrium  barred  storm  profile  which  could  be 
shifted  relative  to  the  prestorm  profile  to  obtain  a  sediment  balance 
between  erosion  and  deposition,  and  thus,  determine  dune  recession. 

Analysis  of  the  experimental  data  provided  a  means  of  calculating 
the  equilibrium  barred  storm  profile  in  terms  of  the  wave  period, 
breaking  wave  height,  surge  level,  and  profile  sediment  grain-size. 
In  most  instances,  the  empirical  expressions  used  in  determining  this 
profile  have  physical  interpretations  in  terms  of  the  incoming  wave 
energy  flux  per  unit  width. 

Finally,  a  method  was  developed  which  corrects  the  maximum  dune 
recession  in  cases  when  equilibrium  conditions  are  not  reached. 
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CHAPTER  I 
INTRODUCTION 

Nature  has  provided  delicate  measures  at  the  land  and  sea  boundary 
on  a  sandy  coast  to  maintain  a  balance  of  forces  under  ever  changing 
wave  conditions.  While  the  sloping  beach  and  beach  berm  provide  the 
first  "line  of  defense"  for  the  absorption  of  normal  wave  energy,  the 
coastal  dunes  are  the  last  zone  of  defense  in  the  advent  of  storm 
(hurricane)  waves  and  the  accompanying  storm  surge  which  succeeds  in 
overtopping  the  berm.  In  essence,  the  beach-dune  system  constitutes 
the  natural  coastal  structure  which  effectively  adjusts  its  shape  in 
the  form  of  erosion  or  accretion  to  cope  with  any  weather  condition. 
Left  in  its  natural  state,  this  system  functions  in  a  most  efficient 
"give  and  take"  fashion  to  provide  a  coastal  defense  mechanism. 
However,  in  recent  years,  increasing  concentration  of  population  and 
industry  in  coastal  areas  has  resulted  in  developments  which  encroach 
upon  the  dynamic  beach-dune  zone  and  therefore  interfere  with  the 
natural  shore  processes  described  above.  The  recent  disastrous 
destruction  of  water  front,  man-made  structures  and  the  magnitude 
of  sediment  required  by  nature  (in  the  form  of  beach-dune  erosion) 
in  adjusting  its  defense  posture  against  Hurricane  Eloise  in  Bay  and 
Walton  Counties  on  the  northwest  Florida  coast  serves  as  a  good  ex- 
ample of  what  could  and  did  happen  to  structures  placed  within  this 
dynamic  zone  (Chiu,  1977;  Frank,  1976). 


-2- 


Realizing  that  the  ever  increasing  development  of  Florida's  beaches 
was  interfering  with  natural  shoreline  processes  and  that  the  Florida 
coastline  is  in  a  high  risk  hurricane  area  (Neumann  et  al . ,  1978),  the 
legislature  established  a  Coastal  Construction  Control  Line  Law  in  Sep- 
tember 1971  (Purpura,  1972).  The  purpose  of  the  law  was  to  prevent 
construction  or  excavation  seaward  of  a  control  line  for  the  protection 
of  the  beach-dune  system  and  for  the  protection  of  upland  structures 
from  storm  activity.  The  establishment  of  such  an  important  line  de- 
pends on  a  number  of  factors,  however.  One  of  the  most  important  fac- 
tors to  be  considered  is  the  landward  erosion  distances  of  the  beach- 
dune  system  that  could  be  expected  during  a  severe  storm.  A  capability 
of  predicting  beach-dune  erosion  distance,  caused  by  storm  activity, 
with  a  reasonable  degree  of  accuracy  would  not  only  be  an  enormous 
benefit  to  the  State  of  Florida,  but  would  also  address  itself  to 
this  great  problem  on  all  the  world's  sandy  coasts. 

The  investigation  of  any  specific  problem  in  the  field  of  applied 
engineering  usually  proceeds  in  the  following  manner: 

1)  Identification  of  the  problem  and  its  important  parameters, 

2)  Development  of  the  physics  of  the  problem  and  formulation  of 
suitable  assumptions  and  simplifications  based  on  observations 
and  common  sense, 

3)  Development  of  a  mathematical  theory  from  the  simplified 
physics  of  the  system  including  all  the  relevant  parameters 
(predictive  in  nature  and  quite  often  highly  empirical),  and 

4)  Theory  verification  using  known  prototype  data  or  by  small- 
scale  model  testing. 

A  review  of  the  currently  available  literature  on  the  subject  of 
surf  zone  dynamics  quickly  leads  to  the  conclusion  that  a  purely 
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theoretical  approach  to  the  problem  of  dune  erosion  based  on  the 
Navier-Stokes  equations  is  presently  beyond  our  grasp.  Among  other 
things,  this  approach  would  require  an  understanding  of  wave  trans- 
formation within  the  surf  zone  and  its  associated  time-dependent  ver- 
tical water  particle  velocity  distribution,  the  turbulent  fluctuations 
and  energy  losses  beneath  a  breaking  wave,  bottom  friction  effects, 
sediment  transport  (both  suspended  and  bed-load)  throughout  the  surf 
zone,  and  consideration  of  random  seas.  At  present,  no  suitable  solu- 
tion exists  for  any   of  these  problems.  However,  experience  tells  us 
that  nature  does  indeed  have  "order,"  and  that  many  natural  phenomena 
can  be  accurately  analyzed  through  empirically  derived  expressions 
relating  the  important  parameters. 

Faced  with  this  insurmountable  theoretical  "wall"  the  investiga- 
tion of  the  dune  erosion  problem  necessarily  becomes  one  of  a  parametric 
study  involving  empirically  verified  relationships  with  the  inclusion 
of  known  physics  when  possible.  A  successful  solution  of  this  type 
would  not  only  provide  a  suitable  engineering  solution,  but  would 
also  help  further  our  understanding  of  the  basic,  underlying  physics 
which  control  this  process. 

There  have  been  two  basic  approaches  to  the  problem  to  date,  the 
gross  simplification  of  the  governing  physics  involved  in  dune  erosion 
and  the  use  of  small-scale  model  testing  to  aid  in  the  parameteriza- 
tion of  the  process.  These  methods  are  discussed  in  Chapter  II,  but 
it  should  be  noted  here  that  none  provide  results  in  agreement  with 
field  measurements.  This  lack  of  verification  is  a  result  of  over- 
simplification of  the  process  to  the  point  that  several  important 
observed  aspects  have  been  omitted,  such  as  offshore  bar  formation 
and  time-dependent  surges. 


Because  current  methods  of  predicting  dune  recession  during 
severe  storms  are  unsuitable,  there  remains  an  immediate  need  for  the 
development  of  an  accurate,  predictive  engineering  solution  of  the 
problem  to  aid  coastal  zone  planners  in  the  orderly  development  of 
the  coastal  communities  without  hindering  the  natural  "storm  defense" 
provided  by  the  beach-dune  system.  This  is  the  primary  goal  toward 
which  this  research  was  directed. 

Important  Factors  and  Assumptions 

Our  present  knowledge  of  the  dune  erosion  process,  derived  mainly 
from  observations  combined  with  common  sense,  make  it  possible  to  list 
the  important  parameters  involved: 

1)  Primary  dynamic  factors: 

a)  the  storm  surge  water  level, 

b)  the  storm  surge  duration,  and 

c)  the  incoming  wave  period  and  wave  height  and  their 
spectral  distribution. 

2)  Primary  geometrical  factors: 

a)  the  configuration  of  the  prestorm  dune,  beach,  nearshore 
and  offshore  zones,  and 

b)  the  grain-size  distribution  on  the  beach. 

No  predictive  model  at  present  incorporates  all  of  these  parameters. 
In  fact,  most  models  are  of  a  steady-state  nature,  which  is  a  condi- 
tion seldom  reached  during  the  short  duration  of  a  storm  event.  This 
indicates  that  any  predictive  model  must  give  transient  beach-dune 
response. 
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Secondary  factors  of  importance  are  given  as: 

1)  Strong  onshore-directed  winds, 

2)  Oblique  wave  attack, 

3)  Heavy  rains, 

4)  Beach  and  dune  vegetation,  and 

5)  Interaction  with  structures. 

While  the  secondary  factors  li'Sted  certainly  have  some  effect 
upon  the  amount  of  erosion  which  could  occur,  it  is  felt  that  the 
small  refinements  resulting  from  inclusion  of  these  factors  into  the 
analysis  are  not  warranted  due  to  the  difficulties  involved. 

In  pursuing  a  solution  to  the  dune  erosion  process,  the  following, 
almost  universally  applied,  assumptions  have  been  invoked  in  this 
study: 

1)  Natural  beaches  with  sand-sized  sediment  distributions; 

2)  Fairly  uniform  offshore  depth  contours  and  straight  beaches; 

3)  No  interactions  with  coastal  structures  or  tidal  inlets; 

4)  A  two-dimensional  onshore-offshore  sediment  transport  situ- 
ation involving  no  net  sediment  losses  or  gains  to  the  system, 
implying  an  absence  of  longshore  currents  (which  field  data 
indicate  is  a  reasonable  assumption  for  severe  storm  condi- 
tions; and 

5)  Sufficiently  high  coastal  dune  as  to  prevent  wave  overtopping 
at  the  peak  surge  level. 


Objectives 

The  overall  objective  of  this  study  has  been  to  develop  the 
capability  of  predicting  the  magnitude  of  beach-dune  erosion  as  a 
result  of  given  storm  conditions.  More  specifically  the  objectives 
include  the  following: 

1)  Develop  and  verify  the  necessary  small-scale  movable-bed 
modeling  relationships  for  use  in  an  experimental  test 
series, 

2)  Determine  the  characteristics  of  the  poststorm  barred  profile, 

3)  Evaluate  the  effects  of  storm  surge  rise  duration, 

4)  Compare  and  quantify  the  relative  importance  of  the  main 
factors  listed  above,  and 

5)  Evaluate  the  effect  of  storm  duration  on  the  transient 
erosion  condition. 


CHAPTER  II 
PREVIOUS  RESEARCH  AND  EXPERIENCE 

The  specific  problem  of  being  able  to  predict  the  amount  of  beach- 
dune  erosion  during  severe  storm  conditions  has  drawn  the  interest  of 
investigators  from  the  fields  of  coastal  engineering,  geology,  and 
land  use  planning  as  well  as  local,  state,  and  federal  agencies.  And 
while  the  destructive  potential  resulting  from  dune  erosion  is  well 
recognized,  descriptions  of  the  phenomenon,  even  in  the  simplest  quali- 
tative evaluation,  provide  little  insight  into  the  actual  physical 
process. 

At  present  there  are  no  verified  methods  for  the  prediction  of 
dune  erosion  during  severe  storms,  even  under  the  limiting  assumptions 
of  a  two-dimensional  process  without  bar  formation,  net  erosion  due  to 
longshore  currents,  or  dune  overwash. 

There  are  three  basic  approaches  to  finding  a  solution  for  a 
hydraulic  process: 

1)  A  complete  physical  analysis  based  on  the  equations  of  fluid 
flow  and  sediment  transport  when  suitable  simplifying  assump- 
tions, based  upon  observation  and  common  sense,  are  utilized; 

2)  A  numerical  representation  of  the  governing  equations  and 
boundary  conditions  which  are  solved  in  an  iterative  fashion 
and  verified  with  extensive  field  data; 
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3)  The  verification  and  operation  of  a  small-scale  model  of  the 
process  carried  out  with  the  hopes  that  the  dominating  forces 
have  been  successfully  scaled  down. 

The  complete  physical  analysis  of  the  problem  is  presently  out  of 
the  question  due  to  the  almost  total  lack  of  understanding  of  the 
fluid  flow  details  within  the  surf  zone.  Likewise,  numerical  methods 
are  limited  by  the  lack  of  detailed  prototype  measurements  before, 
during,  and  immediately  following  a  major  storm  event.  Modeling  in  a 
small-scale  wave  facility  appears  to  be  the  most  practical,  but  this 
too  has  its  drawbacks,  which  will  be  discussed  in  detail  in  Chapter  IV. 

Previous  Simplified  Physical  Solutions 
Edelman's  Method 

It  appears  that  the  first  attempt  at  quantifying  the  amount  of 
dune  erosion  and  the  resulting  recession  of  the  dune  field  was  given 
by  Edelman  (1958).  He  assumed  that  the  poststorm  nearshore  profile 
could  be  approximated  by  a  straight  line  with  a  slope  of  about  1:50 
for  beaches  typical  to  those  found  in  the  Netherlands.  Under  the 
assumption  of  onshore-offshore  sediment  transport  where  the  eroded 
volume  equals  the  deposited  volume,  Edelman's  method  could  be  used  to 
predict  erosion  for  an  actual  beach  profile  as  illustrated  in  Figure 
la,  or  for  "ideal  conditions"  as  shown  in  Figure  lb. 

For  the  actual  beach  profile,  recession,  R,  is  obtained  by  shift- 
ing the  equilibrium  storm  profile  until  erosion  area  equals  deposition 
area.  For  the  "ideal  condition"  the  dune  recession  is  calculated  by 
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Erosion 


-Deposition 
Equilibrium  Storm  Profile 


la:  ACTUAL  PROFILE  APPLICATION. 


_StorrTi_Jide  _Le_vel j_ 


Erosion      i;:  f' 


■Deposition 
Equilibrium  Storm  Profile 

lb:   IDEALIZED  PROFILE  APPLICATION. 


FIGURE  1:   EDELilAN'S  METfiOD  APPLIED  TO  ACTUAL  AND  IDEALIZED 
BEACH  PROFILES. 


-10- 


where 

m  -   beach  slope, 

3  -  D  -  S  +  M, 

y  =  1.3H^  -  S  +  F, 

6  =  1. 3H,  u  -  ]A   , 

F  =  vertical  distance  from  mean  sea  level  to  the  base  of 
the  dune,  and 

S  =  storm  surge  level . 
The  other  assumptions  are  given  as:   (a)  the  width  of  the  surf 
zone  is  defined  by  breaking  wave  height  H,  =  1.5(S)  and  breaking  depth 
h,  =  1.3H,  ,  where  S  is  the  storm  surge  level;  (b)  the  deposition  is 
evenly  spread  over  the  surf  zone;  (c)  the  slope  of  the  poststorm  pro- 
file is  the  same  as  the  original  slope;  and  (d)  the  eroded  dune  toe 
is  slightly  lower  than  the  highest  surge  level  reached.  Both  (c)  and 
(d)  above  are  based  upon  field  observations. 

The  main  weaknesses  of  Edelman's  method  lie  in  the  assumption 
that  the  nearshore  beach  profile  is  a  straight  line  with  the  same 
slope  as  the  prestorm  profile,  and  in  the  assumptions  for  the  limiting 
depth  h,  =  1.3H,  and  the  surf  zone  width.   In  addition,  the  wave 
period  and  the  sediment  grain-size  have  not  been  taken  into  account, 
and  the  assumption  is  made  that  an  equilibrium  is  reached. 

Edelman  (1972)  later  conceded  that  the  beach  profiles  have  a 
more  intricate  shape  than  he  assumed,  and  that  storms  never  lasted 
long  enough  to  produce  the  equilibrium  condition. 
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The  Edelman  method  was  modified  by  Vallianos  (1974)  and  used  to 
design  dune  cross  sections.  One  modification  included  the  recognition 
of  a  flatter  beach  slope  after  the  storm  than  that  of  the  prestorm 
profile. 

A  modified  version  of  Edelman 's  method  has  also  been  used  by  the 
University  of  Florida's  Department  of  Coastal  and  Oceanographic 
Engineering  to  estimate  erosion  distance  for  use  in  establishing 
coastal  construction  setback  lines  for  the  State  of  Florida.  Al- 
though Edelman 's  method  was  used  successfully  by  the  Department  to 
test  an  eroded  profile  in  St.  Johns  County,  Florida,  during  a  three 
day  northeast  storm  in  early  February,  1973,  it  predicts  too  much 
erosion  for  Hurricane  Eloise.  This  points  to  the  influence  of  storm 
duration  since  Eloise  was  a  fast  moving  storm  with  a  peak  surge  last- 
ing only  about  half  an  hour  to  an  hour. 

Dean's  Method 

Besides  Edelman 's  assumptions,  Dean  (1976)  considered  uniform 
energy  dissipation  and  uniform  shear  stress  across  the  surf  zone  to 
establish  an  equation  for  the  equilibrium  beach  profile  given  as 

D 

where 

X  =  the  horizontal  coordinate  directed  onshore  from  the 

breaker  line, 
w  =  the  width  of  the  surf  zone. 
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h  =  the  water  depth, 
h,  =  the  breaking  depth, 
and  a  depends  upon  the  assumed  energy  dissipation  mechanism. 

By  assuming  different  profile  shapes  before  and  after  the  storm, 
and  considering  the  limiting  depth  and  wave  assumptions  of  Edelman, 
Dean  gave  the  following  dimensionless  relationship  for  beach  recession: 

h^ 

in  which 

R'  =  R 


S'  -   S/B  where  B  is  dune  height  above  prestorm  mean  sea  level, 
S  =  surge  level , 

\  -  y^ 
and  the  subscripts  1  and  2  refer  to  prestorm  and  storm  conditions, 
respectively.  Here  again  a  balance  between  erosion  volume  and  deposi- 
tion volume  is  assumed.  Dean  did  a  comparison  between  his  method  and 
Edelman's,  and  while  similar  trends  were  observed,  the  two  methods 
differed  significantly.  He  also  pointed  out  that  a  more  valid  and 
verified  method  is  needed  which  would  include  the  presence  of  an  off- 
shore bar,  the  transient  effects,  and  the  fact  that  surges  will  never 
last  long  enough  to  produce  equilibrium  conditions. 

Previous  Experimental  Research 

While  much  experimental  research  has  been  done  regarding  the 
general  aspects  of  beach  response  to  wave  conditions,  those  which  have 
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addressed  the  specific  problem  of  dune  erosion  during  storm  conditions 
are  relatively  few. 

Qualitative  Model  Tests 

Van  der  Meulen  and  Gourlay  (1968)  performed  two-dimensional  small- 
scale  model  tests  to  examine  the  influence  of  dune  height,  wave  steep- 
ness, and  wave  period  using  a  sediment  grain-size  of  drQ  =  0.22  mm. 
Due  to  their  uncertainty  regarding  the  correctness  of  the  model  law 
implemented  for  these  tests,  the  conclusions  drawn  were  given  as  quali- 
tative observations  which  should  be  considered  for  the  selection  of 
future  scale-model  studies. 

These  observations  include  the  following: 

1)  The  recession  of  a  sand  dune  increases  as  the  height  of  the 
dune  decreases; 

2)  For  given  wave  conditions,  storm  surge  water  level,  and  sedi- 
ment, the  equilibrium  profiles  seem  to  be  independent  of  the 
dune  shape; 

3)  The  dune  erosion  is  a  function  of  wave  steepness,  wave  length, 
sediment  size,  storm  surge  level,  dune  shape,  and  initial 
beach  profile; 

4)  The  dune  foot  recession  is  greater  for  wind  waves  than  for 
regular  waves  with  the  same  energy  and  modal  wave  period; 

5)  Bar  formations  are   less  predominant  for  wind  wave  experiments; 
5)  Dune  erosion  is  very   rapid  in  the  early  stages  of  the  tests 

(a  direct  result  of  using  an  instantaneous  surge  level  rise). 
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All  of  these  observations  are  basically  correct,  and  later  chap- 
ters will  attempt  to  physically  explain  and  quantify  each  point  given 
above. 

Model  Tests  to  Determine  Scaling  Relationships 

Van  de  Graaff  (1977)  conducted  a  series  of  dune  erosion  experi- 
ments with  the  twofold  purpose  of  developing  a  design  criterion  for 
the  Dutch  coastline  and  developing  relationships  for  scaling  between 
model  and  prototype  as  well  as  between  models  of  different  dimensions. 
All  the  tests  were  performed  using  an  Instantaneous  surge  level  rise 
and  a  single  prototype  storm  condition  represented  by  a  Pierson-Moskowitz 
energy  density  spectrum. 

An  empirical  correlation  of  the  data  resulted  in  the  scaling 
relations  of 


and 


where 


N,  ^1 


A     p 


N^  =  morphological  time  scale, 
N  =  horizontal  length  scale,  and 
N  =  vertical  length  scale, 

valid  between  the  ranges  of  26  <  N  <  150  when  the  same  bottom  material 

p 

in  the  prototype  is  used  in  the  scale  model. 

More  extensive  tests  by  Vellinga  (1978)  using  the  same  prototype 
design  conditions  as  Van  de  Graaff  concluded  that  the  best  scale-model 
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relationship  was  achieved  when  the  dimensionless  fall  velocity  parameter 
H/a)T  was  held  constant  between  prototype  and  model.  When  this  could 
not  be  met,  the  profile  could  be  distorted  using  the  relationship 

N.  =  ^^ (4) 

where  a  =  0.5  for  fine  sand  and  N  -   grain  fall-velocity  scale.  The 
morphological  time  scaling  in  this  case  was  given  as 

K   =  N  ^/^  ,  (5) 

which  is  the  same  as  the  hydrodynamic  time  scale. 

Unfortunately,  both  of  these  experimental  series  concentrated  on 
the  development  of  the  scaling  relationships  for  one  prototype  condi- 
tion, and  then  relied  upon  correlations  to  determine  the  distortion 
without  any  physical  arguments  to  support  them.  Hence,  it  is  quite 
possible  that  the  given  relationships  are  valid  only  for  the  tested 
prototype  condition.  Furthermore,  the  use  of  a  fairly  wide-banded 
spectrum,  such  as  the  Pierson-Moskowitz  in  shallow  water,  could  be 
criticized  since  shallow  water  spectra  are  typically  more  narrow-banded. 
This  subject  will  be  touched  upon  again  in  Chapter  IV.  Finally,  the 
testing  of  only  one  prototype  condition  gives  no  information  as  to 
the  effects  of  varying  the  storm  parameters.   In  addition,  the  in- 
stantaneous surge  level  rise  negates  any  type  of  time-dependency 
analysis  which  could  be  carried  out  upon  their  data.   It  will  be  shown 
in  the  present  study  that  the  storm  surge  must  rise  over  a  finite 
period  in  order  to  more  accurately  predict  prototype  response  from 
model  experiments. 
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Ma  (1979)  carried  out  a  series  of  dune  erosion  experiments  at  the 
University  of  Florida  as  a  preliminary  to  the  present  research.  The 
main  thrust  of  the  work  was  directed  towards  identifying  the  scaling 
effects  present  in  small-scale  models.  The  tests  were  designed  holding 
the  dimensionless  fall  velocity  parameter  constant  between  three  dif- 
ferent vertical  scales.  The  prestorm  profile  used  in  the  model  was 
selected  to  represent  an  average  beach-dune  profile  existing  on  the 
coast  of  the  Florida  Panhandle;  and  the  model  series  included  two  dif- 
ferent wave  heights,  wave  periods,  and  storm  surge  levels.  The 
scale-model  law  developed  by  Van  de  Graaff  was  utilized,  and  three 
different  sediment  grain  sizes  were  tested. 

Upon  scaling  the  results  for  both  fine  and  coarse  grains  up  to 
prototype,  Ma  found  good  agreement  between  the  differently  scaled 
models  for  the  same  prototype  storm  conditions.  Unfortunately  the 
comparison  between  model  and  actual  prototype  erosion  caused  by  Hur- 
ricane Eloise  was  not  too  successful  since  the  test  program  did  not 
include  a  time-dependent  surge  level  increase.  However,  this  research 
did  lend  some  credibility  to  the  idea  of  preserving  the  parameter 
H/a)T,  and  it  pointed  out  some  of  the  difficulties  to  be  expected  in 
future  modeling  efforts.   In  addition,  some  qualitative  insight  was 
gained  with  regard  to  the  dune  erosion  process,  and  some  of  the  prob- 
lems with  the  wave  flume  were  uncovered  and  rectified. 

Predictive  Techniques  Based  on  Model  Tests 

An  exhaustive  effort  by  Swart  (1974a,  1974b)  examined  the  results 
of  many  model  tests  carried  out  at  the  Delft  Hydraulics  Laboratory 
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over  the  years.  He  proposed  that  the  beach  profile  be  characterized 
by  three  distinct  zones,  each  with  its  own  transport  mechanism.  These 
zones  are  the  backshore  above  the  wave  run-up  limit  in  which  "dry" 
transport  takes  place,  the  developing  profile  where  a  combination  of 
bed-load  and  suspended-load  transport  occurs,  and  the  transition  area 
where  only  bed-load  transport  takes  place. 

Swart  wisely  avoided  the  internal  mechanisms  of  sediment  concen- 
trations and  bottom  velocities  by  making  a  schematization  of  the 
external  properties  of  the  profile  development  in  order  to  predict 
rates  of  offshore  sediment  transport.  The  final  empirical  relation- 
ships given  by  Swart  provided  a  means  of  determining  both  the  equili- 
brium profile  characteristics  and  the  time-dependent  sediment  transport 
rates. 

Swart  subsequently  realized  that  the  computational  method  he  de- 
scribed was  too  complicated  for  normal  use,  prompting  him  to  modify 
the  technique  in  order  to  ease  the  computations  without  significantly 
affecting  the  results  (Swart,  1975). 

However,  even  with  these  simplifications,  the  method  remains 
fairly  complex  and  has  several  drawbacks  which  preclude  it  from  use 
in  determining  dune  erosion  during  storms: 

1)  The  model  law  used  in  the  testing  has  never  been  fully  veri- 
fied as  representing  actual  prototype  behavior; 

2)  The  tests  used  primarily  a  fixed  water  level  which  does  not 
simulate  storm  surge  rise; 

3)  The  empirically  derived  expressions  offer  no  physical  explana- 
tion for  the  interaction  of  the  parameters  such  as  wave 
height  and  period; 
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4)  The  time-dependent  rates  for  erosion,  when  scaled  to  the 
prototype,  are  assumed  to  be  correct,  but  never  verified; 

5)  Many  of  the  initial  profiles  did  not  represent  near-equilibrium 
beach  conditions; 

5)  Barred  profiles  are  not  predicted  by  the  method  but  are 

usually  observed  in  nature  during  the  offshore  sediment  trans- 
port phenomenon. 

Discussion 

As  can  be  seen  in  the  above  literature  review,  the  problem  of 
coastal  dune  erosion  during  storms  has  been  approached  by  two  different 
methods.  On  one  hand  are  the  simplified  physical  models  offered  by 
Edelman  and  Dean,  and  on  the  other  are  small-scale  model  testing  pro- 
grams seeking  to  provide  some  insight  into  the  process. 

The  physical  models  suffer  due  to  the  lack  of  understanding  of 
the  actual  complex  physics  of  surf  zone  and  the  resultant  effects  upon 
the  coastal  dunes.  There  has  not  even  been  sufficient  field  data  of 
poststorm  erosion  to  indicate  that  these  methods  provide  even  an  order 
of  magnitude  estimate  of  the  dune  erosion. 

The  model  testing  programs,  to  present,  have  been  concerned 
mainly  with  the  establishment  of  suitable  model  laws,  to  be  used  in 
small-scale  testing,  which  would  provide  accurate  results  when  scaled 
up  to  prototype.  In  addition,  a  time-variable  surge  level  increase 
has  not  been  used,  which  later  will  be  shown  to  be  a  pertinent  factor. 

Thus,  it  is  seen  that  this  is  an  area  requiring  immediate  examina- 
tion in  order  to  both  identify  the  important  parameters  and  their 
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relationship  to  the  process,  and  to  provide  an  accurate  and  reliable 
means  of  calculating  the  dune  erosion  as  a  result  of  selected  storm 
conditions. 

Proposed  Method  of  Problem  Solution 

Drawing  upon  the  previous  experience  discussed  above,  the  follow- 
ing steps  are  proposed  in  order  to  logically  examine  and  formulate  a 
solution  method  for  determining  dune  erosion  during  a  storm: 

1)  Derive  a  suitable  small-scale  movable-bed  model  law  based 
upon  sound  physics  and  past  experiments  which  will  provide  a 
basis  for  the  operation  of  a  laboratory  model,  simulating 
storm  conditions  and  the  resulting  dune  erosion. 

2)  Verify  the  model  law  by  the  approximate  reproduction  of  a 
prototype  event  for  which  field  data  are  available. 

3)  Conduct  an  extensive  series  of  model  experiments  in  which  the 
important  storm  parameters  are   varied  within  their  expected 
ranges.  The  series  will  attempt  to  closely  reproduce  actual 
storm  conditions  by  including  such  factors  as  a  varying  storm 
surge  level  and  an  initial  equilibrium  profile. 

4)  The  test  results  will  be  initially  examined  overall  in  an 
effort  to  identify  the  effects  of  the  important  parameters  in 
a  qualitative  sense,  and  to  gain  some  insight  into  the  entire 
process.  From  this  overview  it  may  be  possible  to  recognize 
the  primary  physical  mechanisms  which  influence  erosion. 

5)  From  the  above  step,  decisions  can  be  made  as  to  which 
aspects  of  the  phenomenon  can  be  parameterized  in  such  a 
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way  that  some  physical  justification  for  the  parameters  could 
be  offered. 

6)  Correlations  will  be  sought  for  the  dominant  features  present 
in  the  erosion  process  in  terms  of  the  chosen  parameters,  and 
physical  explanations  will  be  given  when  possible.  While  the 
physical  explanations  will,  for  the  most  part,  be  qualitative 
in  nature,  it  is  expected  that  a  good  amount  of  insight  will 
be  gained  to  further  add  to  the  present  poor  understanding  of 
the  process. 

7)  Finally,  a  proposed  engineering  solution  will  be  constructed, 
based  upon  the  analysis,  with  the  hope  that  it  will  provide 
an  accurate  method  until  the  time  the  process  is  fully  under- 
stood and  analyzed  in  physical  terms.  Along  with  the  proposed 
solution  will  be  a  summary  of  the  restrictions  imposed  by  the 
assumptions  and  an  outline  of  the  range  of  applications  for 
which  the  solution  can  be  deemed  reasonably  accurate. 


CHAPTER  III 
WAVE  TANK  FACILITY 

Air-Sea  Wave  Tank 

The  Coastal  and  Oceanographic  Engineering  wind  and  wave  facility 
was  dedicated  in  1957  and  used  primarily  as  a  wave  tank.  Under  a  grant 
awarded  in  1957,  the  Air-Sea  tank  was  updated  and  improved,  giving  it 
the  capability  of  being  used  for  advanced  research  into  air-sea  inter- 
actions. These  improvements,  along  with  a  detailed  description  of  the 
facility,  are  given  by  Shemdin  (1969). 

The  Wave  Channel 


Figure  2  gives  a  general  layout  of  the  wave  facility  which  is 
6  feet  wide  and  has  a  total  length  of  150  feet.  The  mechanical  wave 
generator  section  occupies  11  feet,  while  the  wave  absorbing  beaches 
occupy  19  feet  of  the  total  length.  The  remaining  120  feet  are  divided 
into  two  bays,  each  34  inches  wide  and  6'4"  deep.  The  maximum 
water  depth  for  wave  generation  is  3  feet.  The  dune  erosion  tests 
were  conducted  in  the  eastern  bay  since  the  entire  length  of  the  test 
section  has  full  height  glass  observation  windows. 

While  the  tank  could  be  filled  from  a  well  at  the  rate  of  500 
gallons  per  minute,  the  water  used  for  this  series  of  experiments  came 
through  a  2-inch  city  water  line  equipped  with  a  calibrated  valve. 
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This  provided  a  controlled  means  of  increasing  the  surge  level  during 
the  experiments,  and  the  water  was  of  greater  clarity  than  the  well 
water,  providing  for  better  observation  of  the  sediment  movement. 

Mechanical  Wave  Generator 

Wave  generation  in  the  Air-Sea  tank  is  achieved  through  an  elec- 
tronically controlled,  hydraul ical ly  driven  wave  paddle  measuring  6 
feet  wide  and  4  feet  high.   Identical  waves  are  generated  in  both  bays. 

The  wave  generator  bulkhead  is  mounted  on  a  carriage  and  is 
driven  by  two  hydraulic  rams  governed  by  hydraulic  servo-valves.  The 
system  provides  independent  control  of  the  top  and  bottom  rams  in  such 
a  way  that  the  bulkhead  can  move  either  as  a  piston  on  the  carriage  or 
as  a  paddle.  Any  combination  of  piston  and  paddle  motion  is  possible. 
For  the  dune  erosion  tests,  the  bulkhead  motion  consisted  of  a  piston 
motion  combined  with  a  paddle  motion  about  the  bottom  of  the  bulkhead. 
This  gave  the  upper  water  particles  a  greater  displacement  than  the 
lower  particles,  which  more  closely  represents  the  vertical  velocity 
distribution  found  in  water  waves. 

Sinusoidal  wave  motion  was  produced  by  sending  an  analog  sine-wave 
signal  from  a  function  generator  to  the  servo-valves  with  the  wave 
amplitude  being  controlled  by  varying  the  input  amplifier  gain. 

Prior  to  the  experiments,  new  servo-valves  and  solid-state  elec- 
tronic controls  were  installed  in  order  to  provide  better  frequency 
response  and  system  reliability.  Throughout  the  entire  testing  pro- 
gram, the  new  system  performed  without  any  problems. 
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Computer 

As  is  the  case  with  most  modern  wave  facilities,  a  digital  computer 
has  been  employed  to  aid  in  the  gathering  and  processing  of  experi- 
mental data.  This  not  only  speeded  up  the  process,  it  also  provided 
a  more  accurate  means  of  data  acquisition,  with  a  smaller  chance  of 
error.  To  guard  against  errors,  methods  were  developed  to  immediately 
check  and  confirm  the  accuracy  of  the  computer-made  measurements. 

Two  further  benefits  of  using  the  computer  are  that  the  speed  of 
profile  measuring  and  recording  allowed  for  a  greater  number  of  experi- 
ments to  be  conducted  within  the  time  constraints  of  the  project,  and 
that  data  reduction  and  analysis  programs  could  be  run  immediately 
following  each  experiment  on  an  instantaneous  "turn-around"  basis. 

In  the  early  stages  of  the  project  the  computer  was  used  to 
generate  random  time-series  to  a  given  power  spectrum,  to  operate  the 
wave  paddle  to  produce  these  random  waves  in  the  tank,  to  monitor  and 
record  wave  data  on  command,  and  finally  to  analyze  the  recorded  wave 
data.  During  the  actual  test  series,  when  regular  waves  were  used, 
the  computer  was  used  only  to  measure  and  record  the  beach  profiles 
and  to  analyze  the  data  for  erosion  volume.  Listings  of  these  programs 
are  given  in  Appendix  A. 

Hardware 

The  computer  used  in  the  experimental  program  was  the  DECLAB-11/03 
with  a  CPU  memory  of  28  K  bytes  and  programmable  in  the  Fortran  IV 
language,  complete  with  diagnostics.  The  system  includes  1)  a  real- 
time clock,  2)  an  analog-to-digital  converter  capable  of  handling  15 
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single-ended  or  8  differential  analog  signal  inputs,  3)  a  digital-to- 
analog  converter  capable  of  outputting  4  different  analog  signals  on 
command,  and  4)  an  RXOl  dual  floppy  disk  drive  for  program  and  data 
storage. 

Two  terminals  were  interfaced  to  the  computer.  A  Tektronix  4006-1 
Graphics  video  terminal  was  used  for  writing  and  editing  programs, 
while  a  TI  Silent  700  hard-copy  terminal  was  remotely  located  nearer 
to  the  wave  tank  for  executing  profiler  programs  and  data  analysis. 

Software 

The  system  software  included  the  standard  components  of:  1)  an 
editing  program  for  writing  or  editing  Fortran  programs,  2)  a  Fortran 
compiler,  3)  a  linking  program  which  creates  the  final  machine  language 
program  including  any  hardware  input/output  commands,  4)  a  peripheral 
exchange  program  for  transferring  programs  or  data  to  different  de- 
vices such  as  tape  or  hard  disk,  5)  an  extensive  Fortran  Scientific 
Subroutine  package,  and  5)  support  routines  for  the  laboratory  peri- 
pherals. 

The  computer  realized  only  a  small  percentage  of  its  potential 
during  the  dune  erosion  experiments,  and  it  is  expected  that  many  more 
applications  will  be  found  in  future  studies  involving  movable-bed 
models. 

Wave  Gauge 

The  wave  gauge  used  for  monitoring  the  waves  in  the  tank  was  a 
capacitance-type  gauge  designed  and  built  by  the  Coastal  and 
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Oceanographic  Engineering  Laboratory  at  the  University  of  Florida. 
Figure  3  shows  the  gauge  in  position. 

The  capacitance  wave  gauge  uses  the  water  as  one  plate  of  a  vary- 
ing capacitor;  hence  the  capacitance  of  the  sensing  portion  of  the 
instrument  changes  as  the  water  level  fluctuates  during  wave  motion. 
The  signal  is  sent  to  the  control  room  in  the  form  of  a  varying  fre- 
quency where  it  is  then  converted  to  an  analog  signal  and  recorded  as 
a  continuous  record  on  a  strip  chart  recorder.  A  detailed  calibration 
at  the  beginning  of  the  experiment  confirmed  the  linearity  of  the 
gauge  (see  Appendix  B),  and  the  resulting  conversion  factor  for  volts 
to  centimeters  was  confirmed  prior  to  every  experimental  run.  The 
gauge  is  electronically  isolated  to  eliminate  any  interaction  between 
other  electrical  signals  present  in  the  laboratory. 

Rail  Mounted  Movable  Carriage 

A  horizontally  moving  instrument  carriage  was  installed  in  the 
eastern  bay  of  the  wave  tank  to  provide  a  variable-speed  platform  on 
which  sensing  instruments  could  be  mounted  for  the  monitoring  and  re- 
cording of  experimental  data.   In  this  research  program  the  carriage 
was  used  solely  for  the  purpose  of  transporting  the  beach  profiling 
instrument  over  the  length  of  the  profile  in  a  controlled  fashion. 
Figure  4  shows  the  carriage  with  profiler  in  place  from  several  angles. 

The  rails  for  the  carriage  have  been  installed  with  a  horizontal 
tolerance  of  ±0.001  inches,  providing  an  excellent  platform  from  which 
to  measure  beach  profiles.  The  carriage  drive  train  is  powered  by  a 
IHP  electric  motor  capable  of  moving  the  carriage  at  speeds  between 
0  to  20  feet  per  second. 
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FIGURE  3:  WAVE  GAUGE, 
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TOP  VIEW 


END  VIEW 
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SIDE  VIEW 


FIGURE  4:      INSTRUMENT  CARRIAGE. 
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A  hand-held  remote  control  unit  allows  the  drive  motor  to  be 
started  or  stopped,  activates  the  automatically  operated  bottom- 
following  probe  on  the  profiling  instrument,  and  also  allows  the 
manual  vertical  positioning  of  the  probe. 

The  horizontal  movement  of  the  carriage  is  monitored  by  a  "follower 
wheel"  which  has  a  ring  of  magnets  equally  spaced  around  it.  A  sensor 
is  activated  each  time  a  magnet  passes  by,  and  this  produces  a  voltage 
variation  similar  to  a  square  wave  between  zero  and  five  volts  as  the 
carriage  proceeds  down  the  rail.  The  computer  can  "count"  the  number 
of  downsteps  in  the  voltage  signal,  and  by  applying  a  conversion 
factor,  it  can  determine  relative  horizontal  carriage  displacement. 
The  approximate  conversion  factor  was  0.42  inches  of  horizontal  travel 
per  voltage  downstep. 

The  horizontal  displacement  can  be  referenced  to  any  initial 
or  final  carriage  location  by  inputting  the  appropriate  value  read 
from  a  measuring  tape  mounted  on  the  upper  portion  of  the  tank  side- 
wall. 

The  signals  from  the  horizontal  position  indicator  and  from  the 
profiling  instrument  were  relayed  to  the  computer  through  a  cable 
connected  between  the  carriage  and  a  junction  box  mounted  on  the  side 
of  the  tank. 

Profiling  Instrument 

The  MK-V  Electronic  Profile  Indicator  was  developed  at  the  Delft 
Hydraulic  Laboratory  in  Holland  for  the  purpose  of  continuously 
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measuring  bed  levels  in  hydraulic  movable-bed  models.  The  instrument 
and  its  mounting  on  the  carriage  is  illustrated  in  Figure  5. 

Probe  Description 

The  instrument  probe  is  placed  vertically  in  the  water  where  a 
feedback  servo-mechanism  maintains  the  tip  of  the  probe  at  a  constant 
distance  (adjustable  between  0.5  mm  to  2.6  mm)  above  the  sediment 
bedform.  When  the  instrument  is  placed  on  the  carriage  and  moving 
horizontally,  the  probe  continuously  follows  the  bottom  profile  con- 
figuration. 

Monitoring  of  a  ten-turn  potentiometer  connected  to  the  vertical 
shaft  holding  the  probe  produces  a  continuous  record  of  the  bottom 
variations  in  terms  of  voltage.  This  voltage  is  easily  converted  to 
vertical  distance  through  the  very  nearly  linear  potentiometer  cali- 
bration constant  determined  to  be  8.25  inches/volt.  This  calibration 
is  shown  in  Appendix  B.  The  vertical  probe  stroke  is  about  36  inches, 
and  the  instrument  is  placed  on  the  carriage  so  that  it  can  reach 
within  1  inch  of  the  bottom  of  the  wave  tank. 

Carriage  Mounting 

A  sliding  track  on  the  carriage  allows  the  profiler  to  be  placed 
in  any  desired  lateral  position  across  the  width  of  the  tank.  For  the 
dune  erosion  experiments  the  profile  was  measured  down  the  tank  center- 
line  with  the  thought  that  this  position  gave  the  most  representative 
profile  and  was  least  affected  by  the  side  walls. 
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PROFILER  IN 
OPERATION 


PROFILER  ON  CARRIAGE 


VERTICAL  PROBE 


FIGURE  5:  BOTTOM  PROFILING  INSTRUMENT. 
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Profjler  Operation 

The  profiler  is  normally  operated  in  conjunction  with  a  computer 
program  which  automatically  records  the  vertical  position  every   time 
the  horizontal  carriage  position  indicator  senses  a  voltage  drop. 
This  translates  to  a  set  of  profile  coordinates  every  0.42  inches. 

Since  it  was  felt  that  this  much  data  were  not  necessary,  the 
program  was  operated  in  a  different  mode  which  allowed  the  operator 
to  take  profile  data  only  where  it  was  deemed  necessary. 

The  recording  of  the  above  water  portion  of  the  profile  was  com- 
plicated by  the  fact  that  the  bottom-seeking  servo-mechanism  is  non- 
functioning when  the  probe  is  out  of  the  water.  However,  since  the 
probe  potentiometer  still  produces  a  voltage  signal,  the  computer 
program  was  written  in  a  way  as  to  permit  the  taking  of  a  data  point 
when  an  external  trigger  was  pressed. 

The  procedure  followed  in  recording  the  above-water  profile  was 
to  place  the  carriage  in  the  desired  horizontal  position,  lower  the 
probe  manually  by  remote  control  to  the  beach  surface,  press  the 
trigger,  and  then  move  the  cart,  using  the  remote  control,  to  a  new 
position  where  the  procedure  was  repeated.  Once  the  profile  was  below 
the  water  level,  the  bottom-seeking  feature  was  activated,  and  the 
points  were  recorded  where  desired  by  stopping  the  horizontal  move- 
ment of  the  carriage  and  pressing  the  trigger  of  the  hand  control. 
The  computer  kept  track  of  the  horizontal  positions  where  data  were 
taken  and  paired  them  with  their  respective  beach  elevations.  The  data 
were  then  stored  on  floppy  disks  along  with  the  necessary  initial 
references  and  conversion  factors  needed  to  convert  the  digitized 
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voltage  to  linear  ranges  and  elevations,  both  referenced  to  the  initial 
intersection  point  of  beach  and  water  level. 

The  operator  was  careful  to  record  all  the  data  necessary  to 
accurately  describe  the  profile  features  and  variations.  The  typical 
profile  measuring  consisted  on  the  average  of  35  to  40  profile  data 
pairs  and  covered  a  horizontal  range  of  14  feet  and  a  vertical  range 
of  12  inches.   In  all  cases  the  measured  profiles  extended  seaward 
past  the  region  of  any  net  sediment  transport. 

While  recording  the  underwater  profile,  the  carriage  moved  at 
about  8  feet/minute,  and,  after  experience  was  gained,  the  operator 
found  that  the  measuring  and  recording  of  the  entire  profile  could  be 
completed  in  under  5  minutes  from  the  time  the  wave  action  was  halted. 


CHAPTER  IV 
MODEL  LAW 

Small-scale  physical  model  testing  of  both  natural  phenomena  and 
man-made  structures  has  long  been  accepted  as  a  useful  engineering  aid 
for  the  analysis  and  prediction  of  the  prototype  behavior.  The  first 
known  hydraulic  scale-model  experiments  were  conducted  by  an  English 
engineer  during  the  period  1752-53  to  determine  the  performance  of 
water  wheels  and  windmills  (Hudson  et  al .  ,  1979).  The  earliest  known 
tests  of  a  river  model  using  a  movable-bed  were  conducted  in  France  in 
1875  and  in  England  by  Reynolds  in  1885.  Since  that  time  great  ad- 
vances have  been  made  in  the  development  of  the  model  laws  which  govern 
hydraulic  similitude,  the  instrumentation  used  in  the  models,  and 
modeling  techniques  in  general.  Perhaps  more  importantly,  many  years 
of  experience  and  basic  research  by  many  fine  hydraulic  laboratories 
throughout  the  world  have  made  the  simulation  of  more  complex  phenomena 
possible,  with  more  confidence  given  to  the  test  results. 

Many  hydraulic  problems  are  fairly  simple  in  nature  and  can  be 
adequately  solved  by  analytical  methods  to  a  reasonable  degree  of 
accuracy.  However,  there  are  many  more  cases  where  the  problem  being 
considered  is  far  too  complex  to  be  handled  analytically.  Often  our 
understanding  of  the  physics  of  the  problem  is  very  limited  or  even 
nonexistant.   In  these  cases  simple  approximations  are  insufficient 
and  other  means  of  obtaining  engineering  solutions  are  necessary.  This 


-34- 


-35- 


is  when  small-scale  modeling  of  a  hydraulic  phenomenon  can  prove 
beneficial  by  aiding  in  the  development  of  a  physical  solution  or  by 
helping  us  to  arrive  at  parametric  relationships  which  can  aid  in  our 
understanding  of  the  hydraulic  process  under  consideration.  This  is 
because  scale  models  of  hydraulic  phenomena  are  essentially  a  means  of 
replacing  the  analytical  integration  of  the  differential  equations 
governing  the  process,  including  initial  and  boundary  conditions.  Often 
model  testing  provides  enough  insight  into  the  physics  of  the  process 
so  that  refinements  can  be  made  to  simple  mathematical  models. 

Physical  models  allow  the  investigator  to  observe  the  process  in 
action  and  distinguish  certain  characteristics  of  the  flow  patterns 
which  might  have  been  neglected  in  the  analytical  approach.   In  addi- 
tion, the  model  provides  a  tight  control  of  the  important  parameters, 
a  control  which  is  not  available  in  the  prototype.  Examples  include 
extreme  events  where  accurate  prototype  data  will  be  lacking,  such  as 
hurricanes,  dam  bursts,  and  flooding.  Using  the  model  to  predict  the 
results  of  these  extreme  events  can  provide  invaluable  information 
which  more  than  justifies  the  high  cost  of  building  and  operating  the 
model . 

Of  all  the  hydraulic  engineering  models  which  can  be  performed, 
movable-bed  scale-model  investigations  of  coastal  erosion  and  coastal 
sediment  transport  phenomena  are  probably  the  most  difficult.   In 
fact,  so  many  different  model  laws  have  been  proposed,  modeling  of 
this  type  should  be  considered  more  of  an  art  than  a  science!  How- 
ever, by  carefully  identifying  the  major  forces  involved,  it  should  be 
possible  to  derive  a  model  law  which  can  be  verified  and  which  will 
provide  reasonable,  quantitative  results. 
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Previous  Scale-Model  Relationships 

Numerous  papers  have  been  written  proposing  similitude  relation- 
ships for  movable-bed  coastal  models.  For  a  rather  complete  bibli- 
ography of  the  subject,  see  page  308  of  "Coastal  Hydraulic  Models" 
(Hudson  et  al . ,  1979).  These  relationships  range  from  those  which 
were  derived  solely  by  theoretical  considerations  to  those  which  were 
established  on  a  strictly  empirical  foundation.   In  each  case  the 
relationship  is  assumed  valid  for  a  given  set  of  specific  flow  charac- 
teristics. For  example,  a  model  law  which  attempts  to  scale  the 
dominant  forces  involved  in  incipient  motion  of  sediment  in  a  tidal 
inlet  over  the  tidal  cycle  will  not  be  valid  for  sediment  transport 
in  the  surf  zone,  where  a  completely  different  set  of  forces  dominate 
the  regime.  For  this  reason  it  is  necessary  to  first  have  an  under- 
standing of  the  dominant  dynamic  forces  at  work  for  a  given  coastal 
hydraulic  problem,  and  then  to  examine  the  proposed  scale-model  rela- 
tionships to  determine  if  these  are  the  forces  being  scaled. 

Perhaps  the  most  thorough  investigation  using  known  relationships 
of  beach  processes  to  determine  the  proper  scaling  law  for  coastal 
movable-bed  models  was  that  of  Fan  and  LeMehaute  (1969).  Their  result 
was  a  table  containing  eight  proposed  scaling  relationships,  each 
being  a  combination  of  three  or  more  derived  similitude  conditions 
based  on  the  known  beach  process  relationships.  One  important  point 
to  note  is  that  all  of  the  proposed  model  laws  satisfied  the  condition 
of  the  time  scale  being  equal  to  the  square  root  of  the  vertical  scale. 
Using  available  data,  a  tentative  model  law  was  chosen  to  be  verified 
by  extensive  experiments. 
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The  experimental  program  was  carried  out  the  following  year  with 
the  results  reported  in  1971  (Noda,  1971).  Noda  conducted  the  experi- 
ments to  determine  the  validity  of  the  proposed  model  law,  and  he  then 
proceeded  to  derive  a  completely  empirical  model  law  based  on  simili- 
tude of  equilibrium  beach  profiles  in  the  breaker  zone.  Unfortunately, 
while  good  reproduction  was  given  in  the  surf  zone,  the  corresponding 
dune  and  beach  erosion  was  not  reproduced.  Hence,  this  model  law  is 
not  really  applicable  to  the  modeling  of  dune  erosion. 

The  literature  appears  to  contain  no  reported  studies  regarding 
model  similitude  relationships  which  were  physically  derived  specifi- 
cally for  the  case  of  dune  erosion  during  storms.  However,  two  recent 
studies  conducted  by  the  Delft  Hydraulic  Laboratory  in  the  Netherlands 
provided  a  totally  empirical  model  law  for  dune  erosion  (Van  de  Graaff, 
1977;  Vellinga,  1978). 

Without  going  into  details,  these  studies  derived  a  model  law  by 
empirical  correlations  of  tests  done  at  different  scale  dimensions 
when  compared  to  a  single  prototype  condition.  Several  questions  can 
be  raised  as  to  the  validity  of  the  derived  relationships: 

1)  Since  very  little  physical  reasoning  has  been  applied  in  the 
model  law  derivation,  it  is  more  than  likely  that  a  prototype 
condition  other  than  that  tested  could  not  be  reproduced. 

2)  Tests  were  conducted  with  a  fixed  surge  level  instead  of  a 
variable  surge  similar  to  what  occurs  in  the  prototype. 
Hence,  there  is  no  solid  evidence  that  the  morphological  time 
scaling  is  correct. 

3)  The  hydraulic  time  scale  was  derived  by  Froude  scaling  and  is 
equal  to  the  square  root  of  the  vertical  scale,  the  same  as 
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Noda's  derivation.   It  will  be  shown  shortly  that  this  is  not 
necessarily  correct  when  modeling  dune  erosion. 

4)  The  Pierson-Moskowitz  energy  density  spectrum  was  used  as  the 
wave  climate  for  the  tests.  This  spectrum  is  representative 
of  deep-water  wave  fields,  while  the  experiments  were  carried 
out  in  water  depths  equivalent  to  70  feet  in  the  prototype. 
At  these  water  depths,  energy  density  spectra  are  decidedly 
narrower  than  the  P-M  spectrum.   In  addition,  prototype  spectra 
contain  wave  grouping  which  is  not  understood,  but  definitely  is 
important.  These  model  tests  did  not  attempt  to  reproduce 

this  phenomenon.  This  point  will  be  discussed  in  the  follow- 
ing chapter  in  further  detail. 

5)  Attempts  to  apply  this  model  law  in  the  initial  verification 
stages  of  the  current  research  failed  to  even  come  close  to 
reproducing  prototype  erosion.   In  fact,  accretion  of  the 
beach  above  mean  sea  level  occurred  for  both  regular  and 
random  wave  conditions.  This  also  will  be  discussed  in  the 
following  chapter. 

These  factors,  combined  with  the  apparent  lack  of  any  other  suit- 
able model  law  for  dune  erosion,  have  led  to  the  following  development 
of  a  new  scale-model  relationship. 

Proposed  Model  Law  for  Dune  Erosion 

The  requirements  for  similarity  between  hydraulic  scale-models 
and  their  prototypes  are  found  by  the  application  of  several  relation- 
ships generally  known  as  the  laws  of  hydraulic  similitude.  These  laws, 
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which  are  based  on  the  principles  of  fluid  mechanics,  define  the 
requirements  necessary  to  ensure  correspondence  between  model  and 
prototype  (Hudson  et  al.,  1979). 

Complete  similarity  between  model  and  prototype  requires  that  the 
system  in  question  be  geometrically,  kinematically,  and  dynamically 
similar.  Geometric  similarity  implies  that  the  ratios  of  all  linear 
dimensions  between  model  and  prototype  are  equal,  kinematic  similarity 
is  similarity  of  motion,  and  dynamic  similarity  between  two  geometric- 
ally and  kinematically  similar  systems  requires  that  the  ratios  of  all 
forces  in  the  two  systems  be  the  same.  From  Newton's  Second  Law  the 
dynamic  similitude  is  achieved  when  the  ratio  of  inertial  forces 
between  model  and  prototype  equals  the  vector  sums  of  the  active 
forces,  which  are  recognized  as  gravity,  viscous,  elastic,  surface 
tension,  and  pressure  in  the  coastal  regime.  An  additional  requirement 
is  that  the  ratios  of  each  and  e^ery   force  must  be  equal.  Or,  in 
equation  form, 

!!li£    !!ai£    l^L,     <!A,     (!A     (U^ 

where  subscripts  p  and  m  are  for  prototype  and  model,  respectively. 
Five  of  these  force  ratios  are  taken  as  independent,  with  one  (usually 
pressure)  being  determined  after  establishment  of  the  others. 

Since  it  is  considered  impossible  to  satisfy  equation  (6)  except 
with  a  full-scale  model,  it  is  necessary  to  examine  the  flow  situation 
being  modeled  to  determine  which  forces  contribute  little  or  nothing 
to  the  phenomenon.  These  forces  can  then  be  safely  neglected  with  the 
goal  of  reducing  the  flow  to  an  interplay  of  two  major  forces  from 


-40- 


which  the  pertinent  similitude  criterion  may  be  theoretically  de- 
veloped (Rouse,  1950). 

For  models  of  wave  action  and  ensuing  sediment  transport  the 
elastic  forces  and  the  surface  tension  forces  are  sufficiently  small 
that  they  can  be  neglected,  provided  that  the  water  wave  length  in  the 
model  is  greater  than  about  10  centimeters.  Since  inertial  forces  are 
always  present  in  fluid  flow,  the  condition  for  dynamic  similitude 
reduces  to  equating  the  ratio  of  inertial  forces  to  the  ratio  of 
either  gravity  forces  or  viscous  forces. 

For  the  particular  case  of  dune  erosion  the  main  area  of  interest 
is  not  the  offshore  zone,  but  the  surf  zone  and  beachface.  Here  the 
waves  rush  up  the  beach  and  then  return  down  the  slope,  eroding  and/or 
depositing  sediment  in  their  wake.  During  this  process,  particularly 
during  severe  storm  conditions,  the  fluid  particle  velocities  near 
the  bed  are  well  in  excess  of  the  critical  velocity  for  incipient 
motion,  and  sediment  is  in  a  state  of  nearly  constant  motion.  Thus, 
over  the  wave  cycle,  the  fluid  process  can  be  idealized  as  unsteady, 
unidirectional,  open  channel  flow  up  a  slope  followed  by  flow  down  the 
slope.  Of  course  this  simplification  can  not  be  used  in  an  analytical 
approach  due  to  the  complexities  involved.  However,  through  this 
visualization,  it  is  easy  to  recognize  that  the  two  major  forces  act- 
ing on  a  sand  grain  are  the  inertia  forces,  due  to  the  turbulent  flow 
fluctuations  near  the  bed,  and  the  nearly  horizontal  component  of 
gravity  acting  parallel  to  the  beach  slope.  The  viscous  forces  are 
small  compared  to  the  forces  due  to  the  turbulent  fluctuations  and 
thus  can  be  neglected  in  this  instance. 
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Before  deriving  the  scale-model  relationships,  it  is  important  to 
point  out  that  the  horizontal  and  vertical  length  scales  will  be  dif- 
ferent, providing  a  distortion  of  the  model.  This  is  necessary  because 
the  reduction  of  size  of  the  sand  grains,  as  required  to  obtain  geo- 
metric similarity,  would  result  in  sand  particles  so  small  that  co- 
hesive forces  not  present  in  the  prototype  would  be  present  in  the 
model.  When  sediment  material  of  nearly  the  same  size  and  specific 
weight  as  found  in  the  prototype  is  used  in  the  model,  distortion  of 
the  model  will  occur  naturally  in  the  mathematical  derivation.  Some- 
times choosing  a  lighter  material  to  use  as  the  sediment  in  the  model 
will  result  in  an  undistorted  model,  but  the  expenses  and  difficulties 
which  arise  are  seldom  worth  the  effort. 

Derivation  of  Dynamic  Similarity  Condition 

For  convenience  it  is  customary  to  introduce  the  notation 


|u  _  value  of  parameter  "a"  in  prototype 
a   value  of  parameter  "a"  in  model 


to  represent  the  scale  ratio  of  a  given  parameter  between  the  proto- 
type and  model.  Using  the  symbols: 

L  =  horizontal  length, 

B  =  horizontal  width, 

D  =  vertical  depth, 

T  =  time,  and 

F  =  force; 
the  fundamental  model  scale  ratios  can  be  defined  as: 
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L    B 
Horizontal  Length  Scale:  N,  =  -j-^  =  h^  (7) 

A    L     D 

m    m 

D 
Vertical  Length  Scale:  N  =  tt^  (8) 


-       ^ 


T 
'T  "  T 


Time  Scale:  N^  =  y^  (9) 


m 

F 
Force  Scale:  K  =  p^  ,  (10) 

m 

where  again  the  subscripts  p  and  m  are  for  prototype  and  model,  respec- 
tively. From  these  four  scales,  all  other  model  scales  can  be  derived. 

Following  the  development  of  Christensen  and  Snyder  (1975),  the 
force  due  to  gravity  in  the  nearly  horizontal  direction  of  the  prin- 
cipal flow  may  be  written  as 


F  =  pg(volume) (sin3) 

where  p  is  the  fluid  density,  g  is  the  gravitational  acceleration,  and 
3  is  the  beach  slope.  For  small  beach  slopes,  sinB  ^j  3  ^i  -p,  so  the 
force  scale  for  gravity  can  be  written  as 


Dp 
P  qjL  B  D J(i:^ 


^  ,  p^p-  p  p  p--p^  ^  plp]N  n2  .  (11) 

g  L)     ^p  -^  ^g  -*  A  u  ^ 

P  g  (L  B  D  ){j^)  ^   "^ 

m^m  m  m  m^  ^L  ' 
m 

The  inertial  force  is  best  represented  as  a  horizontal,  or  nearly 
horizontal,  area  multiplied  by  the  shear  stress  acting  over  this 
area,  or 


^•nertial  ^  ^^^^  ^  ^^^^^  ^^^^^^^ 
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For  the  turbulent  flow  experienced  next  to  the  bed,  which  is  in  the 
rough  range,  the  shear  stress  depends  on  the  rate  of  momentum  transfer 
and  can  be  expressed  as  the  Reynold's  shear  stress  which  is  proportional 
to  the  fluid  density  and  the  time  mean  value  of  the  product  of  a 
vertical  velocity  fluctuation  and  the  velocity  fluctuation  in  the 
direction  of  the  time  mean  flow.  Consequently,  the  inertial  force  is 


F.  =  pu' v' (area)  , 


and  the  inertial  force  scale  ratio  can  be  written  as 


L_  D 


K,  =     ;,P   pP '    (f)^-  (12) 

p^(^)(^)(L^Bj         T 


m   m 


For  dynamic  similitude  requiring  K.^g^^.^-j  =  K  ay^ty'  ^"'^  noting 
that  g  =  g^,  equations  (11)  and  (12)  yield  the  time  scale,  i.e., 

A  p    ^2 


Nj  =  ^  .  (13) 


(N 


M 


Equation  (13)  is  essentially  the  same  as  a  similarity  of  the  Froude 
number  between  prototype  and  model,  when  the  Froude  number  is  based  on 
a  vertical  length  and  a  horizontal  velocity,  i.e.. 


c*  _    horizontal  ,,., 

(g  D   ^-  J^ 
.  vertical 
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Physically,  this  can  be  interpreted  as  a  measure  of  near-horizontal 
displacement  of  a  sand  grain  being  held  up  just  above  the  bed  by 
turbulent  fluctuations.  The  grain  is  moved  horizontally  by  a  velocity 
as  it  falls  back  to  the  bed  vertically. 

It  should  also  be  noted  that  in  the  case  of  an  undistorted  model 
where  N,  =  N  ,  equation  (13)  reduces  to 


y 


Ny  =  (N  )^^^  '  (15; 


which  is  the  time  scale  for  wave  motion  used  in  the  previous  modeling 
attempts.   It  also  arises  from  scaling  of  the  Froude  number.  However, 
the  distortion  required  due  to  the  use  of  beach-size  sand  in  the  model 
ns  that  equation  (13)  must  be  used  for  dynamic  similarity. 


mea 


Sediment  Transport  Similarity  Criterion 

Besides  having  dynamic  similarity  in  the  model,  it  is  necessary  to 
find  some  method  of  determining  the  required  distortion  of  the  model 
arising  from  the  use  of  natural-sized  beach  sand.  The  lack  of  physical 
understanding  with  regard  to  this  question  has  led  many  investigators 
to  propose  a  variety  of  parameters  to  be  scaled  in  the  model  which 
result  in  a  distortion  relationship  (Kemp  and  Plinston,  1968;  Noda, 
1971).  Others  have  used  an  empirical  approach  (Van  de  Graaff,  1977; 
Vellinga,  1978),  as  already  mentioned. 

Currently,  the  most  promising  parameter  used  for  the  prediction 
of  equilibrium  beach  slopes  is  the  dimensionless  fall  velocity,  as 
presented  by  Dean  (1973).  The  physical  significance  of  this  parameter, 
herein  designated  by 
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P  =  4  '  (16) 

where 

H  =  wave  height, 

T  =  wave  period,  and 

u  =  fall  velocity  of  the  sediment, 
is  whether  a  sediment  particle  thrown  into  suspension  by  the  passage 
of  a  wave  will  settle  to  the  bed  during  the  time  that  the  water  par- 
ticle motion  is  shoreward  or  seaward,  resulting  in  onshore  or  offshore 
movement  of  the  particle.  This  parameter  has  proven  to  be  a  good 
predictor  of  onshore-offshore  sediment  transport  (Dean,  1973). 

More  recent  investigations  have  shown  that  the  dimensionless  fall 
velocity  parameter  is  also  a  good  predictor  of  several  other  surf  zone 
features.  Dalrymple  and  Thompson  (1976)  demonstrated  that  this 
parameter  could  be  used  successfully  to  predict  the  foreshore  slope 
when  H  was  given  as  the  deepwater  wave  height.  While  significant 
scatter  was  present  in  the  data  representation,  they  reported  a  much 

better  correlation  than  when  the  same  data  were  plotted  versus  H  /L 

^  o'  0 

and  a  third  parameter  pertaining  to  grain-size.  They  concluded  that 
the  parameter  H  /;^T  should  be  preserved  between  the  model  and  the 
prototype  in  order  to  reproduce  the  same  equilibrium  profile. 

Noda  (1978)  investigated  both  full-scale  and  small-scale  model 
results  for  profile  similarity  and  found  that  a  much  closer  similarity 
could  be  obtained  when  the  H/'jjT  parameter  was  conserved  than  when  wave 
steepness,  H  /L  ,  was  held  constant.  He  also  offered  an  empirical 
relationship  for  the  selection  of  model  grain-sizes  and  concluded  that 
movable-bed  coastal  models  could  be  distorted,  but  the  validity  still 
needed  to  be  confirmed. 
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Savllle*  (1980)  compared  prototype  model  tests  using  sand  against 
small-scale  tests  employing  coal  as  the  sediment.  The  tests  were 
designed  so  that  the  fall  velocity  of  the  sediment  would  be  correctly 
scaled.  His  preliminary  findings  indicated  that  profile  similarity 
was  best  in  the  surf  zone  and  on  the  beachface,  where  settling  velocity 
might  be  expected  to  be  a  major  parameter  affecting  the  modeling. 
The  comparison  seaward  of  the  surf  zone  is  not  as  good,  possibly  be- 
cause a  shear  stress  modeling  relationship  would  be  more  appropriate 
in  this  region. 

As  pointed  out  in  Chapter  II,  comparisons  by  Ma  (1979)  indicate 
that  scaling  between  models  of  different  dimensions  and  sediment  is 
best  achieved  by  preservation  of  the  dimensionless  fall  velocity 
parameter.   In  this  way  there  appears  to  be  a  profile  similarity,  and 
consequently,  erosion  volume  similarity. 

Based  on  this  growing  amount  of  evidence  it  becomes  increasingly 
clear  that  the  parameter  H/uT  should  be  the  same  in  the  model  as  in  the 
prototype  for  profile  similarity.  This  will  also  allow  the  sediment 
grain-size  and  specific  weight  to  be  incorporated  into  the  model  law 
as  a  single  variable,  uj.  This  requirement  becomes 


J_  =  ^IL  (17) 


or 


p  p    mm 


N,  =  /.  (18; 


"Personal  communication  with  Saville  indicates  that  this  research  has 
not  been  completed  as  of  June,  1981,  and  thus  was  not  included  in  the 
17th  Coastal  Conference  Proceedings. 
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Dune  Erosion  Model  Law 

Since  the  time  scale  for  wave  motion  is  the  same  as  the  time  scale 
for  the  resulting  turbulent  velocity  fluctuations,  equation  (18)  can 
be  equated  to  the  dynamic  similarity  criterion  of  equation  (13), 
yielding 


N,     N 

(19) 


^   =  ^ 


Rearranging  gives  the  final  scale-model  relationship  for  the  model 
distortion  as  a  function  sediment  fall  velocity,  i.e., 


3/2 


N,  =  -^-  ,  (20) 


which  combined  with  equation  (13),  given  again  below. 


\ 


7   ,  1/2  ' 


(13) 


provides  the  complete  requirements  for  scale-model  testing  of  dune 
erosion  during  storms  using  a  movable-bed  model.  Actually,  equations 
(13)  and  (20)  should  fulfill  the  modeling  requirements  for  the  deter- 
mination of  any  profile  alterations  due  to  wave  action  in  the  surf 
zone,  not  just  those  due  to  storm  conditions. 

It  is  interesting  to  note  that  equation  (20)  is  identical  to  the 
empirically  derived  distortion  given  by  Vellinga  (1978)  for  fine  sand 
(see  equation  (4),  Chapter  II),  with  the  only  difference  between  the 
model  laws  being  the  time  scale  (equation  (5)  and  compared  to  equation 
(13)). 


Morphological  Time  Scale 

Making  the  morphological  time  scale  the  same  as  the  hydraulic 
time  scale  will  conserve  the  number  of  incoming  waves  per  unit  time, 
thus  conserving  the  incoming  wave  energy  per  unit  time,  and  this  seems 
most  plausible,  in  view  of  the  scaling  of  the  particle  fall  velocity 
and  the  length  scale  distortion.  The  time  scaling  of  the  surge  dura- 
tion should  also  be  to  the  same  scale. 

Unfortunately  not  enough  is  known  about  beach  process  reaction 
times  to  determine  the  morphological  time  scaling  unequivocally.  How- 
ever, results  from  Saville  (1980)  indicate  that  better  time-dependent 
profile  comparisons  are  given  if  the  time  scale  is  distorted  since  the 
small-scale  model  beach  in  Saville' s  experiments  deformed  faster  than 
would  be  expected  when  scaled  up  to  prototype  using  the  time-scale 
relationship 

Nt  -  N  ^/^  (15) 

I     y 

It  is  easily  seen  that  the  proposed  time  scaling  given  by  equa- 
tion (13)  can  be  expressed  as 


Nt  =  qN  ^/^  ,  (21) 


where 

n  =  tA  =   model  distortion. 

Hence,  the  time  scale  has  the  same  distortion  as  the  lengths  in  the 
model.  If  the  model  distortion  is  increased,  the  beach  process  will 
occur  in  a  shorter  time  in  the  model. 
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While  the  above  reasoning  appears  to  be  qualitatively  correct,  a 
further  confirmation  of  the  morphological  time  scale  for  beach  pro- 
cesses is  needed.  However,  the  results  of  the  model  verification, 
given  in  Chapter  V,  indicate  that  the  assumed  time  scaling  is  quite 
reasonable. 

Discussion  of  the  Model  Law 

The  scale-model  relationships  just  derived  have  resulted  from  a 
combination  of  basic  physics  and  an  observed  parametric  relationship 
which  results  in  a  model  distortion  as  a  function  of  the  sediment  grain- 
size  scale  ratio.  The  two  equations  given,  equation  (13)  and  equation 
(20),  are  expressed  in  terms  of  four  variables,  N, ,  N  ,  N^.,  and  N  . 
This  allows  the  experimenter  the  freedom  of  selecting  two  of  the 
scaling  ratios  to  suit  the  model  facility. 

The  main  difference  between  this  model  law  and  previous  attempts 
is  the  derived  time  scaling  relationship.  The  other  model  laws  have 
expressed,  for  the  most  part,  the  time  scaling  as 


N^  -  N 


1/2 


'T    u 

which  usually  arises  as  a  result  of  trying  to  preserve  the  wave 
steepness  parameter  H  /L^.  However,  the  wide  range  of  values  of  wave 
steepness  found  for  transition  between  summer  and  winter  profiles  by 
numerous  investigators  (Johnson,  1949;  Watts,  1954;  Rector,  1954; 
Saville,  1957a;  etc.)  indicates  that  wave  steepness  may  not  be  such 
a  good  parameter  to  preserve. 
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In  contrast,  this  proposed  model  law  preserves  the  dimensionless 
fall  velocity  parameter,  which  has  been  shown  to  be  a  better  indicator 
of  beach  processes.  This  results  in  a  time  scale  given  by 


where 


\i  =  model   distortion 


Likewise,   the  wave  steepness   is   scaled  as 


N  H 

li    0 


P        2Tr        j  2Tr     T 


N        H  ^ 

_U 0  ) 

N  2   LJ 

T  m 


or 


'IN  IN  ^ 


L 


(22) 


One  drawback  that  might  arise  from  this  distortion  is  that  the  reflec- 
tion of  the  incipient  waves  may  become  significantly  greater  in  the 
model  than  in  the  prototype,  and  care  must  be  taken  to  minimize  this 
effect  by  selection  of  scales  which  give  a  small  distortion.  For  ex- 
ample, reflection  is  not  appreciable  for  beach  slopes  milder  than 
about  1:20,  so  a  prototype  beach  having  a  slope  of  1:40  could  be 
modeled  with  a  distortion  up  to  fi  =  2  without  much  effect  due  to  wave 
reflection. 

It  is  possible  to  have  an  undistorted  model,  using  these  proposed 
relationships,  by  the  proper  selection  of  sediment  for  the  model.  This 
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represents  the  ideal  condition,  as  long  as  the  resulting  model  sediment 
size  is  still  outside  of  the  cohesive  sediment  range.  Using  equations 
(13)  and  (20)  with  N,  =  N  ,  the  requirement  becomes 

N  =  N  ^/^  =  N,  ,  (23) 

with  both  the  parameters  H^/L^  and  H/wT  being  conserved.  However,  this 
condition  is  more  often  than  not  impossible  to  satisfy  due  to  the 
fairly  large  length  scale  ratio,  N  ,  required  to  model  typical  sandy 
beaches. 

It  is  worth  noting  that  Battjes'  (1974)  surf  similarity  parameter, 
given  as 

5  -   ^^  ,  (24) 

is  preserved  by  the  proposed  model  law.  Representing  tang  =  -^  as 
before,  and  L  =  -§—  T  ,  then 

0     27T 


D   g  1/2  T 

Si  (%  P  .,       1/2 

P  =  _2 '^p    _  I   y 

^m   D   g  1/2  T       ^A    ' 

L  ^2u^    ,,  1/2 
m       H  ' 
m 

when  (q  /q  )  =  1. 

Substitution  of  equation  (13)  for  N^  yields 
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Thus,  similarities  which  Battjes  noted  were  related  to  the  simi- 
larity parameter  should  be  successfully  scaled  using  the  proposed 
model  law. 

An  interesting  aspect  of  the  power  relationship  in  the  scale-model 
law  between  the  horizontal  and  vertical  scales  is  that  it  is  the  same 
as  was  derived  by  Dean  for  equilibrium  beach  profiles  (Dean,  1977) 
and  further  examined  by  Hughes  (1978).  Dean  proposed  an  equation  for 
the  equilibrium  beach  profile,  given  as 

h  =  Ax^/^  ,  (26) 

where 

h  =  depth  below  mean  sea  level, 

X  =  horizontal  distance  seaward  from  the  intersection  of  the 
beach  and  the  water  level,  and 

A  =  function  of  sediment  characteristics, 
developed  from  the  theoretical  concept  of  uniform  wave  energy  dissi- 
pation per  unit  volume  of  water  in  the  surf  zone.  An  examination  of 
a  large  number  of  profiles  by  both  Dean  and  Hughes  lent  credibility 
to  this  proposed  equation.  Rewriting  equation  (26)  as 

^\-=  f(x,h)  ^  1 

and  requiring  similarity  of  this  parameter  between  prototype  and  model 
gives 

Ax  2/3   A,  2/3 
p  p    _  m  m 

p       m 
or 
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A    x2/3   h 

A^     2/3   h  •  ^^^^ 

m   x^  '     m 
m 

Replacing  x  /x^  with  the  horizontal  length  scale  N, ,  and  h  /h  with 

P  11  ^         ,\      p  m 

the  vertical  length  scale  N  ,  and  rearranging  gives 

^3/2 

"A 
^A  > 


^\         A  3/2  '  (^^^ 


m 

which  results  in  the  same  distortion  power  as  was  found  in  the  model 
law  equation  (20).  In  essence,  the  proposed  model  law  is  also  preserv- 
ing the  equilibrium  beach  profile  as  given  by  equation  (26).   It  is 
quite  possible  that  the  similarity  is  only  coincidental,  but  later 
analysis  in  Chapter  VIII  seems  to  indicate  that  this  is  not  the  case. 

Finally,  as  a  sampler  of  what  follows  in  the  next  chapter,  it 
seems  appropriate  to  mention  that  attempts  to  reproduce  a  prototype 

event  in  the  wave  tank  resulted  in  complete  failure  when  previous  model 

1/2 
laws  based  on  the  time  scaling  N,  =  N  '  were  tried,  while  a  quite 

\  \i 

reasonable  verification  was  achieved  using  the  proposed  relationships. 
This  fact,  plus  the  seemingly  logical  development  of  the  model  law, 
makes  it  appear  that  the  proposed  model  relationships  given  by  equa- 
tions (13)  and  (20)  provide  a  viable  means  of  investigating  dune 
erosion  using  a  small-scale  movable-bed  model. 


CHAPTER  V 
MODEL  VERIFICATION 

Complete  verification  of  a  movable-bed  scale-model  is  perhaps  the 
most  difficult  task  in  the  whole  realm  of  modeling.  The  main  principle 
behind  verification  is  that  of  being  able  to  reproduce  in  the  model  the 
results  of  a  prototype  event  by  the  scaling  of  the  known  parameters  of 
the  event  in  the  model. 

Major  difficulties  arise  in  two  areas: 

1)  In  most  cases,  complete  data  on  all  the  important  parametric 
values  during  the  prototype  event  are  lacking  or  unreliable. 

2)  Small-scale  testing  can  introduce  secondary  effects  caused  by 
the  model  facility  itself.  For  example,  many  erosion  experi- 
ments are  done  in  two-dimensional  wave  tanks,  whereas,  the 
process  may,  in  fact,  be  a  three-dimensional  phenomenon,  com- 
plete with  rip-currents  and  longshore  sediment  transport. 

The  first  of  these  difficulties,  the  lack  of  complete  data,  will 
probably  never  be  resolved  for  the  case  of  dune  erosion  during  storms, 
due  simply  to  the  fact  that  prior  knowledge  of  a  storm  landfall  is 
not  known  in  time  to  install  all  the  necessary  instrumentation  for 
recording  wave  climate,  surge  level  rise,  profile  changes,  etc.  Even 
if  all  instrumentation  were  in  place,  the  severity  of  a  hurricane- 
related  storm  has  the  habit  of  destroying  or  misplacing  the  equipment. 
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The  second  difficulty,  scaling  effects,  can  often  be  resolved  by 
the  careful  examination  of  these  effects  and  the  magnitude  of  change 
they  cause  in  the  model.  In  the  case  of  dune  erosion,  the  assumption 
is  made  that  the  process  is  strictly  that  of  onshore-offshore  sediment 
movement.  While  perhaps  not  totally  correct,  the  general  feeling  among 
most  investigators  is  that  onshore-offshore  motion  is  the  primary 
mechanism  at  work  during  storms,  and  thus,*the  process  can  be  success- 
fully modeled  in  two  dimensions  with  the  hope  that  the  other  effects 
are  small . 

Model  verification  proceeds  loosely  in  the  following  order: 

1)  Select  a  prototype  event  with  as  much  available  data  as 
possible. 

2)  Select  the  model  scales  to  give  the  best  performance  and 
range  of  the  facility.  For  instance,  a  wave  tank  must  be 
able  to  produce  good  wave  forms  over  the  frequency  and  wave 
height  ranges  determined  to  represent  the  upper  and  lower 
extremes  expected  for  the  phenomenon. 

3)  Attempt  to  reproduce  the  prototype  results  by  running  the 
model  with  the  given  input  conditions  scaled  as  accurately 
as  possible. 

4)  After  reproducing  the  desired  end  condition  to  reasonable 
satisfaction,  repeat  the  experiment  to  show  that  the  same 
result  can  be  obtained  for  the  same  conditions. 

5)  Investigate  the  sensitivity  of  the  model  by  examining  the 
effects  of  perturbing  the  input  parameters. 

6)  Select  a  different  set  of  scale  ratios,  as  determined  by  the 
model  law,  and  repeat  the  verification. 
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Verification  of  the  scale  model  is  a  time-consuming  but  essential 
step  in  an  experimental  program.  With  a  verified  model  the  experi- 
menter can  proceed  with  a  carefully  planned  series  of  model  tests 
designed  to  shed  light  on  the  importance  of  each  independent  parameter 
with  reasonable  confidence  that  the  results  accurately  depict  what 
would  happen  in  the  prototype  under  the  same  conditions. 

Selection  of  Prototype  Conditions 

When  Hurricane  Eloise  struck  the  Florida  Panhandle  in  September, 
1975,  its  storm  surge  and  accompanying  destructive  wave  action  eroded 
large  sections  of  the  relatively  undeveloped  natural  beach-dune  system 
of  Walton  County  to  the  right  of  the  storm's  landfall,  in  addition  to 
the  havoc  it  played  with  the  highly  developed  Panama  City  area 
(Chiu,  1977). 

Under  the  Florida  Coastal  Construction  Setback  Line  Program 
(Purpura,  1972),  the  coastal  areas  primarily  affected  by  Hurricane 
Eloise  were  surveyed  about  two  years  before,  with  beach-dune  profiles 
taken  approximately  every   1000  feet  and  offshore  soundings  every  3,000 
feet  (Coastal  and  Oceanographic  Engineering  Laboratory,  1974;  Sensa- 
baugh  et  al . ,  1977).  Immediately  after  the  passage  of  Eloise,  survey 
teams  from  the  Florida  Department  of  Natural  Resources  (DNR)  resur- 
veyed  the  profiles  in  the  areas  most  affected  by  the  storm,  and  these 
195  sets  of  beach-dune  profiles  were  made  available  by  DNR  for  selec- 
tion of  a  suitable  profile  to  be  modeled. 
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Profile  Selection 

Since  the  beach  and  dunes  of  Walton  County  came  closest  to  repre- 
senting the  natural,  undisturbed  shoreline  being  considered  in  this 
study,  these  profiles  were  examined,  and  the  decision  was  made  to  use 
the  beach-dune  profile  designated  as  R-41.  Selection  was  made  based 
on  a  number  of  factors: 

1)  A  25  feet  high  dune,  which  means  no  wave  overtopping  occurred. 

2)  A  fairly  simple  dune  geometry  making  remolding  of  the  dune  in 
the  wave  tank  considerably  easier. 

3)  Profile  R-41  is  located  almost  exactly  on  the  track  of  the 
hurricane's  eye,  which  allows  use  of  existing  methods  for 
calculation  of  peak  storm  surge. 

4)  The  other  Walton  County  profiles  to  the  right  of  the  storm 
landfall  seem  to  exhibit  pretty  much  the  same  erosion  charac- 
teristics as  R-41,  thus  indicating  the  decision  was  fairly 
arbitrary  with  regards  to  poststorm  profile  configuration. 

5)  Sand  samples  taken  in  the  immediate  vicinity  of  profile  R-41 
were  available  from  DNR  for  analysis. 

The  profile  location  is  shown  in  Figure  6,  and  the  profile  selected 
for  verification  tests  is  given  in  Figure  7  with  the  poststorm  profile 
also  shown.  Since  there  was  no  offshore  profile  taken  at  this  range, 
the  offshore  portion  to  be  modeled  was  obtained  by  averaging  profiles 
R-39  and  R-42  where  such  data  were  recorded. 
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FIGURE  5:   PROFILE  R-41  LOCATION  MAP. 
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Profile  Sediment  Analysis 

Sand  samples  taken  in  the  immediate  vicinity  of  profile  R-41  were 
obtained  from  DNR  and  analyzed,  in  order  to  arrive  at  a  representative 
grain  size  for  the  profile.  The  sieve  analysis  used  standard  sieves 
numbered  10,  18,  25,  35,  45,  60,  80,  120,  170,  and  230,  which  gave  a 
linear  progression  in  ^   units.*  The  results  of  the  analyses  are  shown 
in  Table  1  below. 

TABLE  1 
WALTON  COUNTY  SAND  ANALYSES 


Parameter 

Dune 

Mean  High 

Mean  Low 

Profile 

in   (mm) 

Water 

Water 

So 

0.226 

0.215 

0.330 

R-39 

^ 

0.2199 

0.2092 

0.3202 

So 

0.260 

0.280 

0.270 

R-42 

d 
e 

0.2517 

0.2705 

0.2627 

So 

0.260 

R-45 

'e 

0.2511 

In  the  above  table,  d^Q  is  given  as  the  mean  grain-size  diameter 
taken  from  the  distribution  curve  at  the  point  where  50  percent  of  the 
sample  is  finer  than  that  size.  The  parameter  d  has  been  defined  by 


•log2d^^:definition  with  diameter,  d^^,  given  in  millimeters. 
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Christensen  (1969)  to  be  the  effective  grain-size  of  a  nonuniform  natural 
sediment.  It  represents  the  grain-size  of  a  uniform  spherical  sediment 
that  behaves  in  the  same  way  as  the  natural  nonuniform  sediment  from 
which  it  was  derived,  and  as  such,  provides  a  convenient  means  for  ex- 
pressing well-sorted  sediment  distributions.  Calculation  of  d  is  done 
using  equation  (29)  when  the  grain-size  distribution  is  a  straight 
line. 

e  _  2  Cu  In  Cu  /„„, 

d        2      ■  ^^^' 

%0  Cu^  -  1 

dgQ 
In  equation  (29),  Cu  =  -j — ,  Hazen's  uniformity  coefficient,  and  the 

°10 
values  for  d,^,  dr^,  and  d^Q  are  obtained  from  the  sediment  distribu- 
tion as  illustrated  in  Figure  8.  Since  the  prototype  grain-size  dis- 
tributions for  Walton  County  were  nearly  straight  lines  (as  can  be 
seen  in  Figure  9),  d  was  calculated  using  equation  (29).  As  shown 
in  Table  1,  well-sorted,  narrow  grain-size  distributions  have  an  effec- 
tive grain-size  nearly  equal  to  the  mean.  The  average  value  of  d  for 
all  seven  samples  is  0.255  mm,  while  the  average  of  the  three  samples 
of  profile  R-42  is  0.262  mm.  Since  R-42  is  only  1000  feet  from  the 
selected  profile  R-41,  it  was  decided  to  use  the  value  of  d  obtained 
from  the  R-42  samples.  Actually  the  difference  between  that  value  and 
the  overall  average  is  quite  insignificant. 


Surge  Hydrograph  Selection 

Since  no  comprehensive  storm  surge  data  exist  for  Hurricane 
Eloise,  it  was  necessary  to  make  a  close  approximation  for  the  surge 
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hydrograph  using  a  combination  of  existing  prediction  methods  and 
logical  arguments. 

A  first  estimate  of  the  peak  surge  level  can  be  obtained  by  look- 
ing at  the  eroded  dune  profile  shown  on  Figure  7.  It  is  seen  that  the 
dune  toe  is  located  at  a  prototype  elevation  of  about  10  feet  above 
mean  sea  level.  It  is  obvious  that  any  surge  level  above  10  feet 
would  have  eroded  the  dune  more  than  it  did.  Allowing  for  wave  run-up 
would  put  the  peak  surge  level  at  somewhere  between  7  and  9  feet. 

Using  the  nomogram  method  for  prediction  of  peak  surge  for  a  storm 
moving  onshore  given  in  the  Shore  Protection  Manual  (Coastal  Engineer- 
ing Research  Center,  1975),  it  was  possible  to  calculate  a  first 
estimate  of  the  surge  on  the  storm  track,  which  passed  over  the  pro- 
file in  question.  Using  the  storm  values  of 

Vp  -   20  knots--forward  speed, 
P  =  -1.72  inches  H  --pressure  differential, 
R    =  15  naut.  miles--radius  of  maximum  winds, 

obtained  from  the  National  Weather  Service,  the  peak  surge  was  found 
to  be  approximately  7.7  feet. 

Although  the  method  cannot  be  considered  accurate,  it  is  at  least 
in  the  neighborhood  predicted  by  common  sense. 

Finally,  a  computer  prediction  for  peak  surge  using  the  NCAA  SPLASH 
data  is  shown  in  Figure  V-5  of  Pidgeon  and  Pidgeon  (1977),  as  calcu- 
lated for  Hurricane  Eloise.  From  this,  the  surge  elevation  of  about 
7.5  feet  is  given  at  the  location  of  profile  R-41. 

The  above  reasoning,  combined  with  observations  during  the  veri- 
fication runs  as  to  the  effects  of  varying  the  peak  surge,  has  led  to 
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the  final  estimate  of  slightly  over  8  feet  for  the  value  of  peak  surge 
experienced  at  this  location  during  Hurricane  Eloise. 

The  time  history  of  the  surge  level  rise  was  approximated  as  a 
linear  increase  from  mean  sea  level  to  peak  surge  over  a  time  span  of 
12  hours,  a  constant  peak  surge  for  one  hour,  and  a  linear  decrease 
over  6  hours.  This  was  determined  by  examining  the  recorded  surge  time 
histories  associated  with  storms  of  similar  strengths  as  Eloise 
(Harris,  1963). 

While  the  assumed  surge  hydrograph  may  not  be  exactly  what  occur- 
red, it  is  felt  that  it  is  a  very  reasonable  estimate,  and  hopefully  the 
beach  response  in  the  model  will  bear  this  out.  In  fact,  preliminary 
tests  showed  that  the  majority  of  erosion  occurred  around  the  peak 
surge  value;  thus  the  assumed  duration  is  not  quite  as  important  as 
the  peak  elevation.  Later  observations  from  the  experimental  test 
series  confirmed  this. 

Wave  Climate  Selection 

While  no  nearshore  wave  data  were  taken  precisely  at  the  location 
of  R-41,  it  was  possible  to  estimate  the  probable  wave  climate  using 
data  recorded  at  other  locations  during  Hurricane  Eloise. 

A  wave-rider  buoy  operated  by  the  National  Oceanic  and  Atmospheric 
Administration  recorded  wave  data  in  the  middle  of  the  Gulf  of  Mexico 
as  the  hurricane  eye  passed  within  10  miles  (Withee  and  Johnson,  1975). 
At  the  peak  of  the  storm  at  that  location,  it  was  found  that  the 
dominant  spectral  wave  period  was  11  seconds,  and  it  remained  constant 
over  a  five  hour  period.  Peak  significant  wave  height  was  29  feet.  As 
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the  storm  made  landfall,  wave  data  were  also  being  recorded  on  the 
Naval  Coastal  Systems  Laboratory  tower  located  11  miles  offshore  of 
Panama  City  in  105  feet  of  water.  Here  the  modal  wave  period  was  also 
11  seconds  for  the  two  hours  before  power  failure  which  occurred  within 
an  hour  of  the  closest  approach  of  the  storm's  eye  (Pidgeon  and 
Pidgeon,  1977).  At  that  time  the  maximum  significant  wave  height 
was  found  to  be  14  feet. 

From  these  two  sets  of  data,  it  can  be  seen  that  the  dominant  wave 
period  appears  to  be  conserved  while  the  storm  moved  into  shallow 
water.  Thus  a  dominant  wave  period  of  11  seconds  can  confidently  be 
nominated  as  a  prototype  parameter  for  profile  R-41,  located  about 
40  miles  to  the  west  of  the  Panama  City  tower. 

Selection  of  significant  wave  height  is  a  little  less  precise  due 
to  the  effects  of  energy  dissipation  experienced  by  the  wave  field 
during  shoaling.  In  addition,  the  significant  wave  height  should  be 
greater  to  the  right  of  the  track  of  the  storm  with  the  maximum  at 
about  the  radius  of  maximum  winds  (20  miles  in  this  case). 

Some  computer  numerical  models  do  exist  which  attempt  to  predict 
energy  losses  of  spectra  due  to  shoaling  (Hsiao,  1978),  but  the  ex- 
pense and  difficulty  in  running  them  make  their  use,  in  this  case, 
unjustified.*  For  this  reason,  the  significant  wave  height  (SWH)  in 
shallow  water  (50  feet)  at  the  R-41  profile  site  was  estimated  by 
assuming  a  14  feet  SWH  existed  11  miles  offshore,  the  same  as  Panama 
City,  and  that  losses  incurred  in  shoaling  reduced  the  significant  wave 
height  to  around  12  feet.  This  value  represents  the  peak  value  at  the 


''Dr.  O.H.  Shemdin,  Department  of  Coastal  and  Oceanographic  Engineering, 
University  of  Florida,  personal  communication. 
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height  of  the  storm  and  associated  surge.  At  lower  surge  levels,  the 
SWH  value  was,  of  course,  less. 

While  the  chosen  value  for  wave  period  is  probably  very   close  to 
what  actually  occurred  during  the  storm,  not  as  much  confidence  can  be 
given  to  the  selected  SWH.  However,  the  value  of  14  feet  represents  a 
maximum,  while  values  of  under  10  feet  seem  slightly  low  for  a  category 
3  hurricane  at  landfall.  However,  it  was  shown  during  the  experimental 
phase  that  wave  height  perturbations  do  not  significantly  alter  the 
amount  of  dune  erosion. 

Selection  of  Model  Scales 


In  Chapter  IV  the  scale-model  relationships  were  derived,  result- 
ing in  two  equations  in  four  unknowns:  the  vertical  length  scale  ratio 
(N  ),  the  horizontal  length  scale  ratio  (N  ),  the  sediment  fall  velocity 

U  A 

scale  ratio  (N  ),  and  the  time  scale  ratio  (Ny).  Having  the  freedom 
of  selecting  two  of  the  ratios,  it  was  found  that  choosing  N  and  N 
was  the  most  convenient.  Originally  it  was  hoped  that  the  vertical 
scaling  could  be  1:16  (or  N  =  16),  but  this  placed  a  limitation  on 
the  range  of  wave  heights  which  could  be  generated  in  the  tank  before 
extreme  nonlinearities  began  to  dominate  the  wave  motion.  As  a  con- 
sequence, the  vertical  scale  ratio  was  chosen  to  be  N  -   25. 

The  sand  used  to  mold  the  profile  in  the  wave  tank  was  first 
sifted  to  remove  all  sizes  greater  than  0.3  mm.  This  was  done  to  re- 
duce the  effective  grain-size  diameter,  which  in  turn  increases  N  with 

CO 

the  result  of  decreasing  the  model  distortion.  A  second  reason  for 
removing  the  larger  grains  was  to  try  to  prevent  "armoring"  of  the 
foreshore  beach  slope. 
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Two  samples  were  sieve  analyzed,  producing  nearly  identical  dis- 
tributions. A  composite  using  both  sample  results  was  plotted,  and 
the  values  of  d^Q  =  0.152  mm  and  d  =  0.147  mm  were  found.  The  dis- 
tribution of  grain-sizes  for  the  model  sand  compared  very   favorably 
with  the  prototype  distributions,  one  almost  appearing  to  be  a  linear 
offset  of  the  other.  Figure  9  shows  both  the  model  sand  and  the 
approximate  prototype  sand  grain-size  distributions. 

The  values  of  particle  fall  velocity  for  both  the  prototype  and 
model  effective  grain-size  diameters  at  25°C  are  given  in  the  Shore 
Protection  Manual  as : 


o„  =  4.0  cm/sec  for  d  )„  =  0.262  mm 
p  e'^p 


J  =  1.65  cm/sec  for  d  )  =  0.147  mm  , 
m  em 


which  results  in 


N  =  S.  =   2.424 


From  equation  (13)  and  (20)  of  Chapter  IV,  the  scale  ratios  for  hori- 
zontal length  and  time  were  calculated  along  with  the  model  distortion, 
fi  =  Tj— .  These  values  are  listed  in  Table  2  below,  and  were  used 
throughout  the  model  verification  and  the  experimental  test  series. 


TABLE  2 
MODEL  SCALE  RATIOS 


Vertical  Length  Scale 

N 

25 

Horizontal  Length  Scale 

h 

51.55 

Time  Scale 

"i 

10.31 

Model  Distortion 

n 

2.06 
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Requirements  for  Verification 

There  are  three  major  characteristics  of  dune  erosion  which  must 
be  reproduced  in  the  model  test  in  order  to  satisfactorily  obtain 
verification: 

1)  For  the  correctly  scaled  input  parameters  determined  for  the 
prototype,  the  total  volume  of  dune  material  eroded  above  mean 
sea  level  in  the  model  must  approximate  the  same  volume  eroded 
in  the  prototype  when  the  proper  scale  factors  are  applied. 

2)  This  erosion  must  occur  over  the  time  span  of  the  surge  dura- 
tion as  scaled  in  the  model. 

3)  The  poststorm  foreshore  beach  slope  must  duplicate  that  which 
existed  in  the  prototype  when  the  distortion  factor  is  applied. 
By  doing  this,  the  derived  distortion  of  the  model  is  verified. 

Referring  to  Figure  7,  the  above  three  points  essentially  mean 
reproducing  the  poststorm  profile  in  the  region  above  mean  sea  level. 
It  is  not  expected  that  the  berm  feature  between  ranges  -10  inches 
to  +5  inches  will  be  reproduced  since  it  was  probably  formed  by  low- 
steepness  swell  wave  conditions  which  occurred  between  the  end  of  the 
storm  and  the  measurement  of  the  profile  several  days  later.  The  re- 
gion above  this  elevation  was  not  influenced  after  the  surge  elevation 
decreased.  Reproduction  of  the  profile  in  the  surf  zone  is  not 
necessary  because  the  poststorm  wave  climate  also  caused  changes  in 
this  region. 
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Early  Attempts  at  Model  Verification 

The  first  attempts  to  verify  a  model  law  were  carried  out  using 
the  scaling  relationships  as  proposed  by  the  Delft  Hydraulics  Laboratory 
(Vellinga,  1978).  The  first  run  was  made  using  irregular  waves  with  a 
vertical  scale  of  N  =  16.  It  was  seen  that  this  scale  was  too  large 
to  give  proper  waves  which  necessitated  reducing  the  vertical  scale  to 
N  =  25.  Several  tests  were  conducted  at  this  scaling  range  using 
both  regular  and  irregular  waves.   In  all  cases,  sediment  was  moved 
from  offshore  and  deposited  on  the  beach  as  the  surge  level  was  raised. 

At  this  point  it  became  obvious  that  preserving  the  wave  steepness  in 

1/2 
the  model,  which  results  in  the  time  scaling  of  N-r-  =  N  '      =5,  was 

^     T    M 

not  a  proper  technique.  What  perhaps  worked  for  others,  when  an  in- 
stantaneous surge  level  was  applied,  did  not  work  when  the  water  level 
rise  occurred  like  it  does  in  nature.  In  the  latter  case,  the  waves 
had  time  to  transport  sediment  landwards  and  thus  altered  the  beach 
configuration,  which  in  turn  affects  the  waves. 

At  this  point  in  the  verification  program,  a  series  of  undocu- 
mented runs  were  made  for  the  purpose  of  developing  experience  and 
intuition  into  the  effects  of  varying  the  parameters  and  of  gaining 
knowledge  of  the  wave  tank  characteristics,  an  essential  part  of  any 
experimental  endeavor.  Concurrently,  the  scale-model  relationships 
presented  in  Chapter  IV  were  developed,  and  the  resulting  distorted 
profile  was  remolded  in  the  tank. 

Following  the  establishment  of  the  new  model  law,  over  16  attempts 
were  made  at  model  verification  using  an  irregular  wave  field.  None 
of  these  runs  produced  results  which  satisfactorily  fulfilled  the  re- 
quirements for  verification.  An  example  is  shown  in  Figure  10. 
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As  can  be  seen,  the  amount  of  dune  recession  is  good,  and  the 
resultant  beach  slope  is  nearly  right;  but  the  recurring  problem  was 
that  the  erosion  was  not  deep  enough,  and  the  sand  was  not  being 
transported  offshore  far  enough. 

Each  of  the  input  parameters  was  varied  within  reasonable  limits 
to  see  if  the  estimation  of  the  storm  characteristics  was  inaccurate. 
None  of  these  variations  produced  the  desired  result,  although  much 
was  learned  in  the  process.  After  looking  at  possible  causes  for  this 
failure,  from  wave  reflection  to  paddle  response  to  spectral  represen- 
tations and  so  forth,  it  was  finally  decided  that  the  problem  was  in 
trying  to  use  irregular  waves  in  the  model.  The  next  run  after  this 
conclusion  was  conducted  using  regular  waves  which  gave  an  almost  per- 
fect verification  of  the  model.  Since  it  was  only  a  trial,  control 
was  not  very   strict,  and  consequently  this  run  was  not  documented. 
Knowing  that  verification  of  the  model  was  possible  with  regular  waves 
makes  the  explanation  of  the  failures  using  irregular  waves  much  easier. 
These  arguments  appear  in  the  following  section. 

Irregular  Versus  Regular  Waves 

Ideally  the  operation  of  a  small-scale  wave  model  should  include 
a  random-type  wave  field  as  would  exist  in  nature.  In  deep  water  it  is 
possible  to  invoke  the  assumption  of  a  Gaussian  process  in  order  to 
derive  a  suitable  power  spectrum  representation  of  the  wave  heights 
and  frequencies.  However,  as  these  waves  propagate  into  shallow  water, 
they  become  altered  by  such  effects  as  shoaling,  refraction,  and  re- 
flection, thus  turning  the  wave  field  into  a  very  non-Gaussian  phenomenon 
as  individual  waves  begin  to  interact  with  each  other  in  ways  that  are 
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still  not  understood.  One  of  the  main  characteristics  of  shallow-water 
spectra  is  the  presence  of  phase  grouping,  sometimes  known  as  surf 
beat.  This  groupiness  is  a  function  of  both  space  and  time,  with  well 
grouped  waves  at  one  location  becoming  less  well  defined  as  they  pro- 
pagate in  time.  This  must  certainly  have  an  effect  upon  the  erosion 
of  the  beach  and  dune. 

The  narrow-band  spectra  generated  for  use  in  the  model  verifica- 
tion were  developed  under  the  assumption  of  a  Gaussian  process  since 
research  into  shallow  water  waves  is  still  in  its  infancy.  And  while 
the  random  time  series  conformed  to  the  statistical  representation  given 
by  the  desired  wave  spectrum,  there  are  an  infinite  number  of  different 
time  series  which  will  also  conform,  none  of  which  depict  the  wave 
groupiness  present  in  shallow-water  random  wave  fields.  This  is  the 
major  argument  against  using  spectral  wave  representations  in  a  small- 
scale  model  until  more  is  known  about  the  resonance  effect  of  nearshore 
wave  spectra,  and  it  probably  was  the  cause  of  the  failure  to  verify 
the  model  when  spectral  wave  distributions  were  employed.  This  view 
is  also  supported  in  the  reviewer's  comments  presented  in  a  report  by 
Jain  and  Kennedy  (1979),  where  the  statement  is  made  that  small-scale 
movable-bed  modeling  should  be  restricted  to  regular  waves  until  more 
is  known  about  nearshore  random  wave  phenomena.  Further  confirmation 
of  the  importance  of  wave  grouping  was  given  by  Johnson  et  al .  (1978) 
who  found  that  substantially  more  damage  occurred  to  rubble  mound  struc- 
tures for  grouped  waves  than  for  nongrouped  waves  conforming  to  the 
same  power  spectrum  distribution. 

To  further  complicate  matters,  it  has  been  shown  that  the  genera- 
tion of  random  waves  in  model  facilities  has  the  tendency  of  generating 
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long  "parasitic"  waves  because  the  boundary  conditions  for  second 
order  wave  effects  are  not  satisfied  at  the  wave  paddle  (Hansen  and 
Hebsgaard,  1978).  This  fact  was  observed  in  the  verification  process 
on  a  strip  chart  recording  of  a  random  wave  train,  which  illustrated 
a  parasitic  long  wave  with  a  period  between  15  and  20  seconds.  These 
long  waves  are  reflected  and  continue  to  build  in  amplitude.  Methods 
exist  which  will  eliminate  these  long  waves  (Hansen  et  al.,  1980; 
Hudspeth  et  al.,  1978),  but  are  primarily  theoretical  analyses  which 
are  difficult  to  apply  with  the  present  equipment  in  the  Air-Sea  wave 
flume. 

For  the  above  reasons,  plus  the  fact  that  verification  could  be 
reasonably  achieved  using  regular  waves,  it  was  decided  that  the  model 
verification  and  ensuing  experimental  test  series  be  conducted  using 
regular  wave  trains.  Perhaps  future  developments  will  allow  the  use 
of  random  waves  in  small-scale  erosion  tests. 

Verification  Runs 

Two  documented  verification  runs  were  performed  using  regular 
waves  and  reproducing  nearly  the  same  storm  conditions.  The  second 
run  was  required  in  order  to  demonstrate  that  the  results  of  the  first 
run  could  be  reproduced,  and  because  a  slight  wave  period  variation 
was  observed  near  the  end  of  the  first  run.  Both  tests  are  detailed 
below. 

First  Verification  Run 

The  parameters  for  the  first  run  are  given  in  Table  3  with  both 
model  and  prototype  values  given. 
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TABLE  3 
FIRST  VERIFICATION  RUN 


Wave 

period 

(sec) 

Maximum 

wave 

height 

Peak 

surge 

level 

Time  to 

reach  peak 

surge 

Time  at 
peak 
surge 

Time  of 

surge 
decrease 

Prototype 
Value 

11 

8.5  ft. 

8.3  ft. 

12.4  hr. 

1.9  hr. 

5.2  hr. 

Model 
Value    ' 

1.07 

4.1  in. 

4.0  in. 

72  min. 

11  min. 

30  min. 

The  surge  increase  was  approximately  linear  over  the  time  span, 
and  the  wave  period  was  held  constant  throughout.  The  prototype  wave 
height  of  8.5  feet  gave  the  same  wave  energy  density  as  a  narrow-band 
spectrum  with  a  significant  wave  height  of  about  12  feet;  and  while 
the  surge  level  was  on  the  increase,  the  wave  height  was  gradually  in- 
creased until  the  peak  value  was  reached  at  38  minutes  into  the  run, 
or  about  when  the  surge  was  half  its  peak  value.  This  value  was  main- 
tained until  the  surge  level  began  to  decrease  at  which  time  the  wave 
height  was  again  reduced.  This  process  hopefully  typifies  what  happens 
in  nature  over  the  duration  of  a  storm. 

Figure  11  gives  the  measured  profile  taken  down  the  wave  flume 
centerline  at  the  end  of  the  verification  run.  Profiles  measured  near 
the  glass  wall  and  the  back  wall  indicated  that  there  was  a  slight 
lateral  variation  in  the  eroded  profile.  However,  the  erosion  which 
occurred  compared  very  favorably  with  the  prototype,  with  the  exception 
of  the  extra  dune  recession.  A  standard  posttest  check  of  the  wave 
frequency  revealed  that  the  function  generator  knob  had  been  bumped. 


-76- 


0) 

■- 

o 

•  — 

O 

a. 

2 

Q. 

c 
o 

E 
o 

b 

O 

o 
u 

V) 

<u 

^ 

o 

0) 

Q. 

Q. 

> 

O  lO  O  ^O 

(S3H3NI)  N0llVA3n3  IBOOIN 


-77- 


causing  an  increased  wave  period  (equivalent  to  13  seconds  in  the 
prototype)  during  the  latter  portions  of  the  peak  surge  duration.  This 
was  assumed  to  be  the  cause  of  the  extra  recession,  a  fact  later  con- 
firmed by  the  experimental  test  series. 

Second  Verification  Run 

Table  4  gives  the  storm  parameters  for  the  second  verification  run 
which  proceeded  in  much  the  same  manner  as  the  first  run. 


SECOND 

TABLE  4 
VERIFICATION  RUN 

1 

:i  Wave 
■!  period 
|i  (sec) 

Maximum 

wave 

height 

Peak 

surge 

level 

Time  to 

reach  peak 

surge 

Time  at 
peak 
surge 

Time  of 

surge 

decrease 

Prototype 
Value 

11 

8.5  ft. 

8.3  ft. 

12  hr. 

1.4  hr. 

5.2  hr. 

Model     ' 
Value 

1.07 

4.1  in. 

4.0  in. 

70  min. 

8  min. 

30  min. 

As  before,  the  surge  rise  was  approximately  linear  with  the  wave 
period  being  held  constant  throughout.  The  wave  height  was  gradually 
increased  to  its  peak  value  which  occurred  at  about  48  minutes  into 
the  run,  and  the  wave  height  was  decreased  as  the  surge  level  dropped, 
signifying  the  passage  of  the  storm. 

This  time  there  was  little  discernible  lateral  variation  in  the 
eroded  profile;  and  the  eroded  quantity  and  dune  recession  gave  an 
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almost  unbelievable  reproduction  of  the  prototype  event  as  can  be  seen 
in  Figure  12,  which  gives  the  posttest  profile  measured  on  the  flume 
centerline.  Near  the  glass,  the  profile  was  slightly  elevated,  but  not 
by  very  much.  Figure  13a  shows  the  profile  at  the  end  of  the  peak 
surge  duration,  while  Figure  13b  was  taken  at  the  end  of  the  surge 
decrease. 

General  Observations 

Several  observations  which  provide  some  basic  insight  into  the 
erosion  process  were  made  during  the  verification: 

1)  As  the  surge  level  began  to  rise,  there  was  little  change  in 
the  profile  until  the  latter  stages  of  the  increase,  when  more 
of  the  dune  was  exposed  to  the  erosive  wave  action.  This 
points  to  the  importance  of  the  surge  level  in  the  erosion 
process. 

2)  The  erosion  that  had  occurred  by  the  time  the  peak  surge  had 
been  reached  represented  between  80  to  90  percent  of  the  final 
erosion  quantity.  Thus  it  is  seen  that  the  most  damage  is 
done  during  the  surge  rise  and  the  time  the  surge  remains  at 
its  peak  value. 

3)  The  effect  of  lowering  the  surge  level  with  decreased  wave 
energy  had  only  a  minor  effect  on  the  profile,  the  most 
noticeable  being  the  smoothing  out  of  the  profile  as  the  water 
level  dropped.  This  was  caused  by  the  deposition  of  sediment 
over  the  beachface  and  the  transition  from  plunging  breakers 
back  to  spilling  breakers  during  the  surge  decrease. 
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13a:  PROFILE  AT  END  OF  PEAK  SURGE  DURATION. 


13b:  PROFILE  AT  END  OF  SURGE  LEVEL  DECREASE. 


FIGURE  13:  ERODED  BEACH-DUNE,  RUN  2. 
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Discussion 

The  verification  of  the  proposed  model  law  seems  to  have  been 
successful  based  upon  the  above  results.  The  estimation  of  the  storm 
parameters  was  fairly  reasonable  in  view  of  the  limited  data  avail- 
able, and  the  resultant  eroded  profiles  obtained  during  the  model 
testing  closely  resemble  the  actual  prototype  profile  in  terms  of  re- 
cession, eroded  volume,  and  beachface  slope.  The  offshore  portion  was 
not  as  good,  due  to  the  fact  that  this  area  was  altered  by  the  mild 
wave  climate  which  prevailed  between  the  end  of  the  storm  and  the  time 
the  prototype  profile  was  surveyed.  The  subsequent  experimental  series 
demonstrated  that  the  erosion  is  most  effected  by  surge  level  which 
was  probably  the  most  accurate  of  the  estimated  parameters.  Thus  the 
model  would  not  be  too  sensitive  to  incorrect  estimates  in  wave  height, 
as  long  as  they  were  within  reason. 

Of  course  there  is  always  a  need  for  further  verification  since 
the  modeling  laws  have  only  been  verified  on  one  prototype  event,  but 
this  must  wait  until  more  prototype  data  become  available.  Until  such 
time,  the  derived  small-scale  movable-bed  modeling  relations  presented 
in  Chapter  IV  represent  suitable  modeling  laws  for  use  in  small-scale 
investigations  into  dune  erosion  during  severe  storms.  The  relation- 
ships are  based  on  physics  and  observations  and  have  been  verified  in 
the  laboratory  by  the  reproduction  of  storm  conditions  to  a  degree  of 
realism  never  before  attempted,  including  a  variable  surge  hydrograph 
and  a  varying  wave  height. 


CHAPTER  VI 
EXPERIMENTAL  TEST  SERIES 

The  experimental  test  series  to  investigate  the  role  of  the  primary 
factors  involved  in  dune  erosion  by  severe  storms  was  designed  in  such 
a  way  as  to  cover  the  range  of  the  storm  parameters  expected  to  occur 
during  a  hurricane  making  landfall  on  an  undeveloped  portion  of  the 
Florida  coastline.  The  order  of  the  test  runs  was  usually  arranged  so 
that  each  successive  run  reproduced  conditions  more  severe  than  the 
previous  run.  In  this  way  the  regrading  of  the  beach  and  dune  was 
optimized,  and  progressive  patterns  could  be  recognized  early  in  the 
test  series  in  order  to  make  changes  in  the  parameters  if  necessary. 
This  chapter  details  the  experimental  test  series,  the  parameters 
tested,  and  the  procedures  followed. 

Profile  Selection 

The  selected  profile  for  the  model  series  was  chosen  to  be  repre- 
sentative of  an  average  Florida  beach.  From  an  analysis  by  Hughes 
(1978),  it  was  found  that  a  curve  fit  of  the  profile  form 

h  -   Ax2/3 

on  over  400  actual  beach  profiles  from  the  Florida  each  coast  and 
panhandle  indicated  a  modal  value  of  A  =  0.15  when  both  depth  (h)  and 
range  (x)  are  given  in  feet.  Since  these  profiles  probably  exhibited 
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the  characteristics  of  near-equilibrium  beach  profiles,  it  was  decided 
to  adopt  this  equation  for  the  below  mean  sea  level  portion  of  the 
profile.  This  represents  a  better  choice  than  a  plane-sloped  beach 
since  it  is  desirable  to  begin  the  model  tests  with  at  least  near- 
equilibrium  conditions.  In  this  way,  the  resultant  changes  would  more 
closely  duplicate  nature. 

Modeling  the  profile  equation  using  the  previously  derived  rela- 
tionship obtained  by  comparing  equations  (20)  and  (28),  i.e., 


A„    w„  2/3 
m    m 


where 


CO  =  grain-size  fall  velocity, 
p  =  prototype  value,  and 
m  =  model  value, 
resulted  in  the  equation 

where  both  h^  and  x„  have  units  of  inches.  Hughes  also  observed  that 

mm  ^ 

beyond  a  certain  depth  the  offshore  profile  became  distinctly  steeper, 
most  probably  a  relic  feature  from  the  rising  sea  level  over  the  past 
4000  years  (or  from  previous  severe  storms).  An  inspection  of  prototype 
profiles  indicated  that  this  portion  could  be  well  represented  by  a 
straight  line  with  a  prototype  slope  of  1:40  extending  seaward  from 
about  the  13  foot  depth. 

The  prestorm  beachface  was  represented  by  a  straight  line  having 
a  prototype  slope  of  1:10  with  a  dune  toe  elevation  of  7  feet  above 
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mean  sea  level.  These  values  were  chosen  after  discussions  with  an 
individual  familiar  with  Florida  beaches.*  While  slightly  steeper  than 
the  average  beachface  calculated  for  Bay  and  Walton  Counties  in  the 
Florida  panhandle  (Chiu,  1977),  it  was  felt  that  the  values  used  were 
a  good  representation  of  Florida  beaches  on  the  whole. 

The  dune  was  chosen  with  an  initial  height  of  20  feet  above  mean 
sea  level  in  the  prototype  with  a  dune  face  slope  of  1:2.  This  slope 
is  approximately  the  angle  of  repose  for  seciments  in  the  range  of 
grain  sizes  found  in  coastal  dunes,  and  this  slope  was  retained  when 
the  dune  height  was  increased  later  in  the  test  series. 

The  same  sediment  size  as  was  used  during  the  verification  tests 
was  used  for  the  model  series.  Hence  the  prototype  profile  repre- 
sents one  with  an  effective  grain-size  of  d  =  0.27  mm.  For  this  reason 
all  the  scaling  relationship  values  given  for  the  model  verification 
in  Table  2  still  apply  for  the  test  series.  Figure  14  gives  the  ini- 
tial profile  used  for  the  experimental  series  along  with  the  model 
dimensions  and  slopes  as  determined  by  the  scale-model  relationships. 

Storm  Parameters 


The  storm  parameters  which  were  varied  during  the  test  series 
were  as  follows: 

1)  Wave  period  (10,  12,  and  14  seconds  in  the  prototype); 

2)  Wave  height,  when  measured  before  shoaling  began  in  the  tank 
(5,  8,  10  feet  in  the  prototype); 

3)  Peak  surge  level  above  mean  sea  level  (8,  11,  14  feet  in  prototype] 

*Dr.  T.Y.  Chiu,  Department  of  Coastal  and  Oceanographic  Engineering, 
University  of  Florida,  personal  communication. 
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4)  Instantaneous  surge  level  rise  as  opposed  to  finite-time 
surge  increases; 

5)  Surge  level  rise  duration  (4,  7,  15  hours  in  the  prototype); 
and 

6)  Prestorm  dune  height  (20,  24,  28  feet  in  the  prototype). 
Regular  waves  were  used  throughout  the  entire  experimental  program 

which  began  November  12,  1980,  and  ended  January  20,  1981.  The  38 
experimental  runs,  involving  a  time-dependent  surge  level  increase, 
approximated  the  surge  as  a  linear  increase  over  time.  Table  5  lists 
both  the  model  and  prototype  values  for  the  parameters  involved  in 
the  test  series. 

The  experimental  test  runs  can  be  divided  roughly  into  five  dis- 
tinct series,  as  noted  on  Table  5.  Series  1  compared  three  different 
wave  heights  at  three  different  surge  levels  for  a  given  wave  period 
and  dune  height.   In  addition,  the  first  8  runs  of  the  series  investi- 
gated the  effects  of  different  time  spans  in  bringing  the  surge  level 
from  mean  sea  level  to  its  peak  value.  The  first  run  used  an  in- 
stantaneous peak  surge  level  in  order  to  see  if  a  finite  surge  is 
even  required  to  obtain  the  final  equilibrium.  Runs  36  and  37  of  series 
4  repeated  this  test  using  different  surge  levels  and  initial  dune  - 
heights. 

Series  2  repeated  series  1  at  a  different  wave  period,  the  only 
difference  being  that  the  effect  of  surge  level  rise  duration  had 
been  established  so  further  tests  on  this  variation  were  not  required. 

Test  series  3  was  an  abbreviated  version  of  series  1  and  2  using 
yet  another  wave  period,  but  maintaining  the  same  initial  dune  height 
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elevation.  In  this  series  only  two  different  wave  heights  at  3  dif- 
ferent surge  levels  were  compared. 

Series  4  and  5  repeated  some  of  the  earlier  tests  using  dif- 
ferent values  of  dune  height. 

Time-variation  plots  of  the  experimental  profiles  for  each  run 
are  shown  in  Appendix  C. 

General  Observations 

A  few  general  observations  are  worth  noting: 

1)  As  the  surge  level  began  to  rise,  very   little  change  in  the 
original  profile,  outside  of  ripple  formation,  occurred  until 
the  surge  level  reached  a  height  where  the  wave  uprush  reached 
the  toe  of  the  original  dune.  This  was  typically  between 

15  and  20  minutes  into  the  run,  depending  upon  the  rate  of  rise 
being  used.  After  this,  dune  erosion  occurred  quite  rapidly. 

2)  In  test  series  1,  as  the  surge  approached  its  peak  value, 
some  unevenness  in  the  profile  across  the  width  of  the  tank 
was  observed.  However,  this  was  usually  rectified  by  the 
wave  action  before  the  peak  surge  value  was  reached.  The 
cause  of  this  is  the  fact  that  the  wave  frequency  was  close 
to  the  fundamental  mode  of  sloshing  in  a  basin,  and  the  wave 
crests  could  be  seen  exhibiting  a  rocking  motion  as  they 
propagated  down  the  tank  toward  the  beach.  On  the  gentle 
prestorm  slope,  this  resulted  in  uneven  uprush;  but  as  the 
sediment  moved  offshore  and  built  the  bar  system,  the  slope 
became  steeper  which  acted  to  decrease  this  effect.  The 
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profile  then  quickly  became  uniform  across  the  width  of  the 
wave  flume. 

It  is  strongly  felt  that  this  phenomenon  had  little  in- 
fluence on  the  results.  The  other  series  sometimes  exhibited 
this  unevenness,  but  to  a  very  small  degree. 

3)  Wave  reflection  was  most  evident  in  test  series  1,  and  when 
the  bar  formation  became  prominent,  the  monitored  wave  reflec- 
tion sometimes  approached  the  value  of  15  percent.  When  this 
occurred,  the  test  was  temporarily  halted  and  the  water  was 
allowed  to  settle  before  the  test  was  resumed.  The  other 
series  at  higher  values  of  wave  period  showed  little  sign  of 
reflection,  but  subsequent  data  analysis  indicated  that  re- 
flection was  present  in  all  runs.   It  will  be  seen  that  this 
problem  is  overcome  in  the  data  analysis  by  considering  only 
the  properties  of  the  transmitted  wave  instead  of  the  deep 
water  characteristics. 

4)  Dune  erosion  occurred  very  rapidly  as  the  surge  level  in- 
creased. Typically  70  to  90  percent  of  the  final  "equilibrium 
erosion"  occurred  by  the  time  the  surge  reached  its  peak  value. 
However,  the  process  was  not  one  of  passing  through  discrete 
steps  of  equilibrium.  Chapter  IX  presents  a  more  complete 
analysis  of  these  observations. 

5)  The  model  waves  initially  were  breaking  in  a  spilling  manner 
until  the  surge  level  increased  to  the  point  that  the  profile 
was  steep  enough  to  cause  plunging  breakers.  The  steepening 
of  the  profile  was  caused  by  a  bar  formation  offshore. 
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6)  It  appeared,  at  least  in  the  two-dimensional  case,  that  all 

the  sediment  that  went  into  building  of  the  offshore  bar  during 
the  experiments  was  derived  solely  from  the  coastal  dune. 
This,  however,  may  not  be  the  case  for  offshore  bar  forma- 
tion in  the  absence  of  a  storm  surge. 

7)  As  the  erosion  was  taking  place,  it  was  not  uncommon  for  the 
dune  to  be  substantially  undercut  before  the  weight  of  the 
overhanging  dune  brought  it  plunging  down  into  the  water. 
This  was  caused  by  the  surface  tension  of  the  wet  sand  in  the 
model.  Since  surface  tension  was  not  scaled  in  the  model, 
this  effect  would  not  be  present  in  nature  where  the  weight 
of  the  sand  would  cause  the  dune  to  shear  off  before  much 
undercutting  could  occur.   In  order  to  provide  more  reasonable 
results,  it  was  at  times  necessary  to  tap  the  top  of  the  dune 
to  collapse  it  when  the  undercutting  had  become  excessive. 
Surface  tension  did,  however,  work  to  the  model's  advantage 
since  the  nearly  vertical  sheared  dune  face  closely  resembled 
those  found  in  nature  where  the  roots  of  vegetation  aid  in 
maintaining  a  nearly  vertical  dune  face  during  storm  erosion. 

Typical  Test  Procedures 

Procedures  are  important  in  any  experimental  research  in  order  to 
provide  consistent  results  and  to  help  eliminate  human  error.  When- 
ever possible,  procedures  and  standard  practices  should  be  established 
in  advance  of  the  testing  program.  The  following  is  a  brief  outline  of 
some  of  the  procedures  and  techniques  used  during  the  experimental  phase 
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of  this  project.  By  listing  them  here,  future  investigators  might  be 
able  to  identify  the  pitfalls  and  develop  better  methods  for  their 
work. 

Regrading  of  the  Prestorm  Profile 

After  every  run  it  was  necessary  to  regrade  the  beach  to  the  ori- 
ginal prestorm  profile.  This  profile  was  marked  on  the  outside  of  the 
flume's  glass  wall  with  an  indelible  marking  pen.  It  was  found  that 
the  easiest  method  was  to  drain  the  water  from  the  tank  and  then  to 
immediately  regrade  the  profile  while  the  sand  was  still  wet.  This 
gave  better  sand  compaction  and  a  smoother  profile.  The  grading  was 
always  begun  offshore  and  proceded  toward  the  beach.  Sighting  the  pro- 
file marking  through  the  glass,  the  profile  was  graded  using  a  timber 
2x4  cut  to  slightly  less  than  the  tank  width.  The  profile  was  kept 
uniform  by  checking  with  a  "bubble-type"  level  across  the  width  of  the 
tank.  Since  it  was  found  that  the  grading  proceeded  faster  than  the 
tank  could  be  filled  with  water,  the  inlet  water  valve  was  open  while 
the  grading  was  being  done.  The  process  of  draining,  regrading,  and 
filling  the  tank  could  be  accomplished  within  U  to  2  hours. 

Pretest 

Before  the  start  of  each  test,  the  experimental  values  for  the  run 
were  recorded  in  a  journal,  the  initial  water  level  was  checked,  the 
peak  surge  point  gauge  was  set,  the  wave  gauge  calibration  was  checked, 
the  electronics  were  wanned  up,  the  wave  frequency  was  set  or  checked, 
and  the  stopwatch  was  wound  and  set. 
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Beqinninq 

After  the  hydraulic  system  was  brought  up  to  pressure,  small  waves 
were  generated,  and  the  inlet  water  valve  was  opened  to  its  predeter- 
mined position  to  provide  the  nearly  linear  surge  level  increase  over 
the  stipulated  time  span.  At  this  point,  the  timing  of  the  experiment 
began.  After  a  few  minutes  the  wave  amplitude  was  increased  to  that 
stipulated  for  the  run.  While  a  steadily  increasing  amplitude  and 
period  over  the  surge  rise  would  more  closely  duplicate  nature,  it  was 
felt  that  this  introduced  too  many  factors  into  the  present  analyses. 
During  the  surge  rise,  the  wave  amplitude  was  frequently  checked  and 
adjusted  to  compensate  for  the  increasing  hydrostatic  pressure  on  the 
wave  paddle.  Any  irregularities  were  noted  in  the  journal. 

Peak  Surge 

As  the  time  for  peak  surge  approached,  the  wave  height  was  closely 
watched,  and  the  water  level  was  monitored  on  a  preset  point  gauge 
located  in  the  stilling  basin  to  the  rear  of  the  tank.  When  peak  surge 
was  reached,  the  water  valve  was  closed,  the  wave  generation  was  termi- 
nated, the  time  was  noted,  and  the  profile  was  measured  and  recorded. 
The  computer  immediately  listed  all  values  for  spot  checking  for 
accuracy. 

Equil ibrium 

At  this  fixed  peak  surge  level,  the  same  wave  condition  was  then 
run  until  near-equilibrium  appeared  to  be  reached.  This  time  varied 
between  30  and  45  minutes,  with  profiles  being  recorded  every  15  minutes, 
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Equilibrium  was  assumed  when  the  dune  toe  position  had  not  shifted  and 
the  nearshore  beachface  did  not  appear  to  be  changing  over  a  period  of 
at  least  5  minutes.  As  a  check  of  these  criteria,  one  test  was  allowed 
to  continue  running  for  several  hours  with  very  little  noticeable 
change. 

Posttest 

Immediately  following  the  measurement  of  the  final  profile,  the 
drainage  valves  were  opened,  and  the  tank  was  drained.  Simultaneously, 
the  profiles  were  plotted  and  compared  with  previous  results  to  con- 
firm trends  recognized  earlier.  After  the  drainage  valves  were  closed, 
the  water  valve  was  opened,  the  profile  regraded,  and  the  water  level 
brought  up  to  its  initial  level. 

Using  this  method,  it  was  possible  to  conduct  two  tests  in  one 
working  day;  however,  other  factors  sometimes  prevented  this. 

The  series  of  photographs  shown  in  Figure  15  depict  the  dune  ero- 
sion that  occurred  during  run  41.  The  times  given  are  those  from  the 
start  of  the  surge  level  increase. 
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CHAPTER  VII 
QUALITATIVE  ANALYSIS  OF  EQUILIBRIUM  PROFILES 

It  is  instructive  to  first  examine  the  resulting  experimental 
profiles,  recorded  when  near-equilibrium  conditions  were  reached,  for 
general  trends  before  attempting  an  analytical  parameterization.  This 
provides  an  insight  into  the  relative  importance  of  each  of  the  tested 
parameters,  as  well  as  providing  a  view  of  the  overall  general  mechan- 
ism involved  in  the  process.   In  this  comparison,  only  the  last  profile 
recorded  for  each  run  will  be  examined  in  an  effort  to  visualize  the 
final  steady-state  condition.  The  time-dependent  trends  will  be 
discussed  in  Chapter  IX.  For  some  aspects  of  the  comparison,  it  is 
possible  to  involve  a  fairly  large  number  of  profiles.   In  these 
instances,  representative  examples  are  given  with  the  remainder  shown 
in  Appendix  D. 

Surge  Level  Rise  Duration 

Paramount  to  the  entire  experimental  series  was  the  establishment 
of  the  effects  of  the  surge  level  rise  duration.  Definite  knowledge 
on  this  subject  had  to  be  obtained  before  the  test  series  could  pro- 
gress into  the  investigation  of  other  parameters. 

Firstly,  a  comparison  was  made  between  the  equilibrium  profile 
obtained  for  given  wave  conditions  when  the  peak  surge  level  was  reached 
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instantaneously  and  a  profile  formed  under  the  same  wave  conditions  but 
with  a  finite-time  surge  level  rise. 

Figure  16  compares  run  1,  where  the  peak  surge  was  fixed  before 
the  waves  began,  and  run  2  where  the  surge  level  was  raised  to  its  peak 
value  over  45  minutes  and  then  allowed  to  run  until  near-equilibrium 
was  established.  As  can  be  seen,  the  variable  surge  level  run  experi- 
enced more  dune  recession,  apparently  because  sediment  can  be  moved 
further  offshore  while  the  water  rises  and  thus  allows  more  sediment 
to  be  eroded  from  the  dune.   In  the  case  of  the  instantaneous  surge, 
the  seaward  slope  of  the  offshore  bar  is  steeper,  indicating  that  the 
depth  limit  for  significant  sediment  transport  prohibits  sand  from 
being  moved  as  far  offshore.  Later  tests  conducted  using  a  different 
dune  height  and  two  additional  surge  levels  confirmed  this  trend  as 
shown  in  Figure  17.  Details  of  the  runs  are  given  in  Table  5  of 
Chapter  VI.  Thus,  the  first  important  conclusion  to  be  drawn  is  that 
a  time-dependent  surge  level  rise  allows  more  erosion  to  take  place 
than  could  occur  when  the  surge  rise  is  instantaneous.  This  comes 
about  quite  simply  because  the  limiting  depth  for  sediment  transport 
covers  a  wider  range  of  the  profile  when  the  water  level  is  gradually 
increased.  Since  time-dependent  surge  increases  occur  in  nature,  then 
they  must  also  be  reproduced  in  the  model  experiments. 

The  next  question  to  be  answered  about  surge  level  rise  regards 
the  time  span  over  which  the  rise  takes  place.  Runs  numbered  2  through 
7  were  conducted  with  this  question  in  mind.  Figure  18  compares  runs 
2,  3,  and  4  which  were  conducted  using  the  same  wave  condition  but  with 
surge  level  rise  times  of  45,  85,  and  25  minutes,  respectively.  As  can 
be  seen,  the  profiles  are  very   similar  at  peak  surge  but  show  some 
variation  30  minutes  later. 
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Since  these  runs  proved  inconclusive,  runs  5,  6,  and  7  were  con- 
ducted using  a  different  wave  height.  As  can  be  seen  from  Figure  19, 
the  profiles  show  good  similarity  with  the  single  exception  of  the  peak 
surge  profile  of  run  6  representing  a  rise  time  of  25  minutes,  which 
simulated  the  prototype  condition  of  an  8  feet  surge  rise  over  a  4  hour 
period. 

Experience  tells  us  that  a  surge  rise  as  rapid  as  given  by  run  6 
is  highly  unlikely  to  occur  in  nature.  However,  the  close  similarity 
shown  between  runs  5  and  7,  which  represent  prototype  surge  rise  times 
of  12  and  8  hours,  respectively,  point  to  the  fact  that  the  finite  surge 
increase  duration  is  not  an  important  consideration,  as  long  as  it  is 
within  reasonable  ranges  expected  in  the  prototype. 

From  this  evidence  it  is  concluded  that  the  analysis  of  the  dune 
erosion  need  not  be  concerned  with  the  duration  time  of  the  surge  level 
rise,  but  instead  only  with  the  peak  surge  value  and  the  length  of 
time  the  peak  value  is  maintained. 

For  the  remainder  of  the  tests,  the  surge  level  rise  was  accom- 
plished usually  between  50  and  60  minutes  (8.6  to  10.3  hours,  prototype). 
It  was  felt  that  this  was  fairly  typical  of  most  hurricanes  making 
landfall. 

The  Effect  of  Dune  Height 

Comparisons  of  the  final  recorded  profiles  from  test  series  2,  4, 
and  5  were  made  to  investigate  the  effects  of  different  dune  height 
elevations  on  the  erosion  process.  Figure  20  gives  the  comparison  in 
two  instances  where  wave  height,  wave  period,  and  peak  surge  level  were 
the  same.  The  remaining  comparisons  are  given  in  Appendix  D.l. 
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It  is  seen  that  while  the  nearshore  section  and  the  dune  recession 
show  significant  variation,  the  offshore  profiles  from  the  bar  trough 
seaward  are  almost  identical.  While  some  variation  is  present,  due  to 
experimental  effects,  the  strong  trend  is  toward  similar  offshore  bar 
features.  Also  in  e^ery   case  it  is  plainly  evident  that  the  higher 
dunes  have  less  recession  than  the  lower  dunes.  These  two  observations 
lead  to  the  interesting  hypothesis  that,  for  a  given  prestorm  beach- 
offshore  profile  and  given  storm  condition,  there  appears  to  be  a 
maximum  sediment  storage  capacity  in  the  offshore  bar  for  equilibrium 
conditions.   In  other  words,  the  storm  will  continue  to  move  sediment 
from  the  dunes  to  the  offshore  bar  until  this  "storage  capacity"  has 
been  reached  at  which  time  the  dune  recession  ceases.  Since  higher 
dunes  provide  more  sediment,  the  dune  recession  would  be  less  in  that 
case. 

This  concept  implies  that  there  exists  an  equilibrium  barred  storm 
profile  configuration  which  is  characteristic  of  wave  conditions,  pre- 
storm profile,  and  sediment  grain-size.  Once  this  profile  is  deter- 
mined, it  would  be  possible  to  shift  the  profile  horizontally  along 
the  prestorm  profile  (relative  to  peak  surge  level)  until  a  sediment 
balance  between  erosion  and  deposition  is  achieved,  thus  determining 
the  equilibrium  erosion  which  could  occur.  Obviously  a  higher  dune 
will  mean  less  dune  recession.  The  peak  surge  level  will  be  the  most 
important  factor  since  the  barred  storm  profile  is  taken  relative  to 
the  peak  surge.  The  analytical  development  presented  in  the  following 
chapter  is  based  upon  this  hypothesis. 
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The  Effect  of  Storm  Surge 

The  effect  of  the  peak  surge  level  upon  the  dune  erosion  at  near- 
equilibrium  conditions  is  seen  qualitatively  in  Figure  21  which  gives 
comparison  of  runs  of  the  same  wave  period,  wave  height,  and  dune 
height  but  with  different  peak  surge  levels.  Typically  it  is  seen 
that  a  surge  level  increase  of  3  feet  in  the  prototype  is  responsible 
for  dune  recessions  on  the  order  of  30  to  40  feet  in  the  prototype, 
for  this  particular  beach-dune  configuration.  Another  interesting 
aspect  is  that  the  portion  of  the  profiles  seaward  of  the  offshore  bar 
seem  to  coincide,  indicating  that  the  "storage  capacity"  of  the  off- 
shore bar  is  directly  related  to  the  surge  level,  as  it  should  be. 
However,  this  point  will  not  be  stressed  because  the  other  comparisons, 
shown  in  Appendix  D.2,  do  not  always  illustrate  this  trend  as  strongly. 
The  two  following  sections  show  that  while  wave  height  and  wave  period 
have  some  effect  upon  the  dune  recession,  neither  exhibits  the  large 
variations  caused  by  different  surge  levels. 

Therefore,  it  can  be  concluded  that  the  surge  level  acts  as  a 
vehicle  for  exposing  greater  portions  of  the  dune  to  the  incoming  wave 
energy,  and  that  the  surge  level  remains  the  single  most  important 
parameter  involved  in  the  process. 

The  Effect  of  Wave  Height 

The  comparisons  of  runs  involving  the  same  wave  period,  surge 
level,  and  dune  height,  but  different  wave  heights,  are  given  in 
Figure  22  and  Appendix  D.3.  An  increasing  wave  height  has  the  effect 
of  increasing  the  dune  recession  but  to  a  much  lesser  extent  than  an 
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increase  in  surge  level.  The  interesting  thing  to  note  is  that  in- 
creasing the  wave  height  lengthens  the  distance  between  the  eroded 
dune  toe  and  the  crest  of  the  offshore  bar,  as  well  as  increasing  the 
depth  of  the  bar  crest.  These  two  points  are  quantified  in  Chapter  VIII 
to  help  establish  some  important  characteristics  of  the  equilibrium 
barred  profile. 

The  Effect  of  Wave  Period 

Comparisons  between  runs  where  the  wave  period  was  varied,  while 
the  surge  level,  wave  height,  and  dune  elevation  were  held  constant, 
indicate  that  increasing  the  wave  period  increases  the  amount  of  dune 
recession,  but  once  again,  not  by  a  large  amount  (see  Figure  23  and 
Appendix  D.4).  This  effect  is  attributed  to  the  fact  that  the  longer 
period  waves  run  further  up  the  beach  and  are  thus  capable  of  eroding 
more  of  the  dune.  This  relationship  is  quantified  in  the  following 
chapter  and  used  in  determining  the  equilibrium  barred  profile.  The 
form  of  the  remaining  portion  of  the  profile,  especially  in  the  off- 
shore sector,  does  not  seem  to  be  particularly  affected  by  variations 
in  wave  period. 

Nearshore  Beach  Profile 

One  additional  profile  comparison  was  conducted  by  superimposing 
the  near-equilibrium  profiles,  formed  under  the  same  wave  condition  at 
different  surge  levels,  in  such  a  way  that  the  common  reference  point 
was  the  intersection  of  the  beachface  and  the  still  peak  surge  level. 
Representative  examples  are  shown  in  Figure  24  with  the  remainder 
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given  in  Appendix  D.5.  From  these  comparisons  it  is  seen  that  the 
nearshore  profile  shape  is  independent  of  the  surge  level  and  initial 
prestorm  profile,  while  the  profile  seaward  of  the  offshore  bar  is 
definitely  related  to  the  initial  profile  and  surge  level.   It  is  also 
encouraging  to  note  that  the  location  of  the  dune  toe  and  the  location 
of  the  trough  and  crest  of  the  bar  show  a  good  similarity  under  the 
same  wave  conditions.  This  is  an  indication  of  the  repeatability  of 
the  experimental  results. 

From  the  above  qualitative  look  into  the  effects  of  the  primary 
factors  involves  in  dune  erosion,  it  has  been  possible  to  draw  several 
conclusions  which  give  valuable  insight  into  the  erosion  of  dunes  and 
formation  of  barred  profiles.  Using  these  conclusions  as  guidelines, 
it  is  now  possible  to  begin  the  quantitative  analyses  using  the 
experimental  test  results.  However,  a  general  discussion  on  the  over- 
all topic  of  dune  erosion  is  warranted  in  order  to  place  everything  in 
its  proper  perspective. 

Dune  Erosion  in  General 


Besides  the  assumptions  listed  in  Chapter  I,  there  are  several 
additional  aspects  regarding  storms  and  the  resulting  erosion  of  the 
beach-dune  system  which  need  to  be  examined  in  the  hopes  of  recognizing 
some  of  the  limitations  or  failures  of  the  ensuing  analysis. 

Waves 

The  experimental  test  series  were  run  using  regular  waves  of  a 
"nonforced"  type,  i.e.,  gravity  was  the  only  force  determining  wave 
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motion.  It  is  well  known  that,  as  these  waves  shoal,  they  steepen  and 
break,  most  often  in  a  spilling  or  plunging  mode.  The  transition 
between  these  two  types  of  breaking  depends  upon  wave  steepness  and 
beach  slope  (Galvin,  1972)  with  plunging  occurring  for  steeper  waves 
and  beach  slopes.  Galvin  relates  these  two  steepnesses  through  a 
parameter  given  as 

"  -    °  (31) 


°   L^(tan3)2  ' 

where  tang  is  the  beach  slope.  Since  this  parameter  can  be  easily 
shown  as  being  conserved  by  the  model  law  given  in  Chapter  IV,  it  is 
reasonable  to  believe  that  the  model  waves  behaved  in  a  similar  fashion 
to  prototype  waves  when  it  came  to  breaking.  Thus  the  type  of  wave 
breaking  which  occurred  during  the  various  stages  of  the  experimental 
run  should  be  the  same  as  what  would  be  expected  in  nature.  For  ex- 
ample, the  waves  initially  were  spilling  breakers  due  to  the  mild 
prestorm  beach  slope,  but  changed  to  the  plunging  type  as  the  surge 
level  rose  and  the  offshore  bar  began  to  form,  resulting  in  a  steep- 
ening of  the  offshore  beach  slope. 

During  storms,  however,  the  waves  are  "forced"  by  strong  onshore- 
directed  winds.  So,  instead  of  breaking  strictly  under  the  force  of 
gravity,  the  waves  never  reach  their  maximum  steepness  because  the 
wind  blows  the  tops  off,  causing  the  waves  to  appear  to  be  of  the 
spilling  type. 

This  "forced"  condition  may  cause  the  final  erosion  obtained  in 
the  prototype  to  be  somewhat  less  than  indicated  by  the  small-scale 
model  since  plunging  waves  cause  more  turbulence  and  hence  move  more 
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sediment.   In  addition,  the  offshore  bar  features  are  less  predominant 
for  spilling  breakers. 

For  the  present  it  must  be  assumed  that  this  "forcing"  condition 
either  contributes  little  to  the  final  erosion,  or  that  it  is  offset 
by  another  natural  phenomenon  not  reproduced  in  the  laboratory  experi- 
ments, such  as  spectral  distribution. 

Spectral  Distribution 

The  reasons  that  regular  waves  were  chosen  over  a  narrow-banded 
spectrum  of  waves  for  use  during  the  model  tests  are  given  in  Chapter 
V,  the  main  reason  being  the  inability  to  reproduce  the  "surf  beat 
resonance"  phenomenon.   In  nature,  where  the  incoming  wave  field  ex- 
hibits a  variety  of  wave  heights  and  frequencies,  it  would  be  expected 
that  the  bar  feature  would  be  a  little  less  defined  than  with  regular 
waves  since  the  different  height  waves  break  at  different  depths,  thus 
acting  as  a  "smoothing"  agent.  The  presence  of  resonance,  however, 
helps  to  make  the  bar  a  little  more  well  defined  in  addition  to  causing 
the  wave  runup  to,  at  times,  progress  further  up  the  beach  than  in  the 
case  of  regular  waves.  This  additional  uprush  causes  more  erosion  and 
might  possibly  compensate  for  the  "forced"  wave  condition  discussed 
above. 

Until  more  is  understood  about  nearshore  wave  spectra  and  forced 
wave  conditions,  it  must  be  assumed  that  the  results  obtained  by  the 
experimental  series  are  close  to  reproducing  natural  events  with  the 
error  being  in  favor  of  overestimation  of  the  dune  recession. 
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Nonbarred  Profiles 

There  are  several  instances  where  the  formation  of  an  equilibrium 
barred  profile  would  not  occur: 

1)  When  conditions  are  such  that  the  waves  always  break  in  a 
spilling  fashion,  a  bar  formation  would  be  less  likely. 

2)  If  the  prestorm  beach  slope  is  "jery   mild  and  the  surge  level 
rise  is  not  too  severe,  the  lengthy  surf  zone  would  dissipate 
a  considerable  amount  of  energy,  leaving  less  energy  available 
to  move  sediment  and  build  up  the  bar.  While  some  erosion 
would  take  place  and  equilibrium  conditions  could  be  obtained, 
the  development  in  the  following  chapter  is  not  addressed  to 
this  particular  condition  since  it  does  not  represent  the 
severe  erosion  condition. 

3)  A  similar  condition  to  that  described  above  is  the  case  of 
fairly  low  elevation  dunes.  For  a  substantial  surge  level 
rise,  the  dune  recession  will  be  tremendous  as  the  storm  moves 
sediment  offshore.  It  may  be  that  the  lengthened  surf  zone 
could  reach  its  "storage  capacity"  without  the  formation  of 

an  offshore  bar.  The  verification  of  the  small-scale  model 
discussed  in  Chapter  V  is  an  example  of  such  a  case. 

Discussion 

In  view  of  the  above  points,  it  is  necessary  to  restrict  the 
definition  of  an  equilibrium  profile  by  stating  that  the  equilibrium 
storm  profile  is  one  in  which  all  the  incoming  wave  energy  is  completely 
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dissipated  over  the  length  of  the  surf  zone.  This  profile  may  or  may 
not  include  an  offshore  bar  feature. 

Concentration  for  the  remainder  of  the  analysis  will  be  given  to 
barred  storm  profiles  since  they  were  most  prevalent  in  the  experimental 
tests  and  represent  the  most  severe  prototype  conditions. 

Summary 

Before  becoming  entangled  in  the  quantitative  analysis  of  the  next 
chapter,  it  is  beneficial  to  briefly  review  the  conclusions  reached  in 
this  overview  of  the  final  near-equilibrium  test  profiles. 

1)  In  model  testing,  it  is  necessary  to  employ  a  time-dependent 
surge  level  increase  of  sufficient  duration,  in  order  to  ob- 
tain more  realistic  erosion  results. 

2)  There  appears  to  be  an  offshore  bar  "sediment  storage  capa- 
city" for  a  given  prestorm  profile  and  set  of  storm  condi- 
tions. This  is  in  the  form  of  an  equilibrium  barred  profile 
taken  relative  to  the  peak  surge  level. 

3)  Dune  erosion  can  be  found  by  shifting  this  barred  profile 
horizontally  until  the  erosion  area  equals  the  deposition 
area. 

4)  The  storm  surge  level  is  the  most  important  parameter  in  de- 
termining erosion,  while  wave  height  and  period  have  a  lesser 
effect. 

5)  The  profile  between  the  dune  toe  and  offshore  bar  is  indepen- 
dent of  surge  level,  and  thus,  independent  of  prestorm  profile 
configuration. 


CHAPTER  VIII 
EQUILIBRIUM  BARRED  STORM  PROFILE 

It  was  seen  in  the  preceding  chapter  that  an  equilibrium  barred 
profile  can  be  formed  under  storm  conditions  if  the  process  is  allowed 
to  continue  indefinitely.  This  profile  appears  to  exhibit  three 
distinct  regions:  the  nearshore  profile,  the  bar  feature,  and  the 
offshore  profile.  While  very   few  storms  are  of  sufficient  duration  to 
allow  establishment  of  the  final  equilibrium  state,  the  equilibrium 
condition  lends  itself  to  analytical  development  since  the  transient 
effects  may  be  discarded. 

Once  the  equilibrium  profile  for  a  given  set  of  storm  conditions 
has  been  established,  it  will  be  possible  to  then  determine  the  maximum 
expected  erosion  and  dune  recession  which  could  occur.  This  is  done 
by  simply  shifting  the  storm  profile  relative  to  the  prestorm  profile 
until  a  sediment  balance  is  achieved  between  the  eroded  quantity  and 
the  deposited  quantity.  With  the  "equilibrium  erosion"  solution  in 
hand,  it  will  then  be  possible  to  investigate  the  amount  of  erosion 
which  occurs  from  the  time  that  peak  surge  is  reached  until  the  time 
that  a  near-equilibrium  is  established.  The  quantity  of  erosion  which 
occurs  over  this  time  span  approximately  equals  20  percent  of  the 
total  final  equilibrium  volume. 

The  actual  physics  of  the  bar  formation  are  understood  only  in 
the  overall  general  description  of  the  phenomenon.  The  increased  wave 
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energy  and  water  level  rise  associated  with  storm  activity  act  to  upset 
the  beach  equilibrium  established  during  normal  weather  conditions 
because  the  surf  zone  cannot  dissipate  this  increased  energy  over  its 
length.  As  a  consequence,  a  large  volume  of  water  is  transported  onto 
the  beachface.  The  upflow  provides  a  quick-sand  effect  which  sets  the 
stage  for  rapid  removal  of  beach  material,  and  the  downflow  drags  sedi- 
ment seaward  from  the  surface  of  the  beachface.  The  large  volume  of 
water  reaching  the  beach  must  return  to  the  sea  in  order  to  maintain 
a  balance.  This  return  flow  is  sediment-laden  and  continues  to  carry 
its  load  offshore  until  the  flow  velocity  is  sufficiently  reduced  to 
the  point  that  the  load  is  dropped.  Through  this  mechanism  an  offshore 
bar  is  formed  by  the  deposition  of  sediment  removed  from  the  beachface 
and  dune.  This  bar  continues  to  grow  and  migrate  seaward  until  its 
crest  becomes  high  enough  to  steepen  and  break  the  majority  of  the 
incoming  storm  waves.  When  this  dissipation  mechanism  is  fully 
developed,  and  the  incoming  wave  energy  is  completely  dissipated  through 
breaking  and  turbulent  propagation  over  the  extended  length  of  the 
surf  zone,  the  equilibrium  barred  storm  profile  is  then  established. 
From  this  description,  the  importance  of  the  bar  in  providing  the 
ultimate  defense  against  beach-dune  erosion  can  be  seen,  making  the 
inclusion  of  the  bar  feature  in  the  analysis  imperative. 

Breaker  Depth  Versus  Breaker  Height 

When  dealing  with  waves  in  the  surf  zone,  it  is  an  almost  uni- 
versally accepted  practice  to  relate  the  water  depth  at  breaking  (h,  ) 
to  the  breaking  wave  height  (H,  )  through  the  simple  relationship 

\  =   1-28  H^   ,  (32) 
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as  derived  by  Munk  (1949)  using  a  modified  solitary  wave  theory.  Sub- 
sequent observations  of  both  lab  and  field  data  by  Iversen  (1952a, 
1952b),  Galvin  (1969),  and  Goda  (1970),  just  to  name  a  few,  have 
established  a  slight  dependency  of  the  ratio  h^/H,  on  beach  slope  and 
deep  water  wave  steepness  H  /L  .  However,  the  data  were  somewhat 
scattered,  and  it  is  felt  that  the  use  of  the  simpler  form  of  equation 
(32)  gives  satisfactory  results  as  illustrated  in  Figure  25  which  was 
reproduced  from  an  unpublished  report  by  James  H.  Balsillie  of  the 
Florida  Department  of  Natural  Resources,  Bureau  of  Beaches  and  Shores. 

Munk  also  derived  an  expression  for  the  breaking  wave  height  in 
terms  of  deep  water  wave  characteristics  which  is  given  by 


Hk       H„  -1/3 

H^-0.30(l^)     ,  (33) 

0  0 


where 


H  is  deepwater  wave  height  and 
L  =  -§—  T  =  deepwater  wave  length. 

0    2lT  ^ 

Van  Dorn  (1978)  gives  a  slightly  different  expression  for  H, 
given  as 

Hk      H^  -1/3 

H^=0.4(-p^)     ,  (34) 

0         0 

which  varies  from  Munk's  theoretical  derivation  only  by  the  constant, 
and  which  seems  to  fit  the  available  data  a  bit  better. 

For  the  model  experiments,  the  values  of  H,  obtained  from  equa- 
tions (32),  (33),  and  (34)  are  shown  in  Table  6  along  with  the  value 
of  H  ,  calculated  using  linear  wave  theory  from  the  intermediate  depth 
where  the  wave  height  was  measured  in  the  tank.  The  calculation  allows 
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FIGURE  25:   BREAKER  HEIGHT  VERSUS  BAR  DEPTH. 
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TABLE  6 
BREAKING  WAVE  HEIGHTS 


Run 

# 

Period 
(sec.) 

Surge 
Level 

(in.) 

Wave 
Height 

"o 
(in.) 

Breaking  Wave  Height--H,  (in.) 

Eq.  (32) 

Eq.  (33) 

Eq.  (34) 

2 
8 
9 
10 
11 
12 
13 
14 

0.970 

3.84 

2.97 
3.95 

1.96 
1.96 

2.39 
2.9 

3.19 
3.87 

5.28 

2.95 
3.93 
4.92 

1.58 
2.29 
2.89 

2.39 
2.89 
3.35 

3.18 
3.85 
4.47 

5.72 

2.94 
3.19 
4.89 

1.73 
2.16 
2.64 

2.38 
2.88 
3.35 

3.17 
3.84 
4.46 

15 
16 
17 

1.163 

3.84 

3.08 
4.11 
5.13 

1.96 
1.96 
1.51 

2.77 
3.35 
3.89 

3.69 
4.47 
5.19 

18  1 

19  ; 

20  t 

5.28 

3.06 
4.08 
5.09 

1.69 
2.12 
2.82 

2.76 
3.34 
3.88 

3.68 
4.45 
5.17 

21  1 

22  1 

23  ! 

6.72 

3.04 
4.05 
5.06 

1.62 
2.25 
2.73 

2.75 
3.33 
3.86 

3.66 
4.44 
5.15 

24  1 
25 

26  ' 
27 

28  ! 

29  ; 

1.357 

3.84 

3.14 
4.19 

1.77 
2.21 

3.11 
3.77 

4.15 
5.03 

5.28 

3.13 
4,18 

1.18 
2.26 

3.11 
3.77 

4.14 
5.02 

6.72 

3.12 
4.16 

1.65 
2.13 

3.10 
3.76 

4.13 
5.01 

30  ! 

31  1 

3.84 

3.08 
4.10 

1.86 
1.91 

2.77 
3.35 

3.69 

4.47 

32  i   1.163 

33  1 

5.28 

3.06 
4.08 

1.62 
2.20 

2.76 
3.34 

3.68 
4.45 

34  1 

35  1 

6.72 

3.04 
4.05 

1.82 
2.30 

2.75 
3.33 

3.66 
4.44 

38  1 

39  1 

1.163 

5.28 

3.06 
4.08 

1.62 
1.89 

2.76 
3.34 

3.68 
4.45 

40 

41  li 

6.72 

3.04 
4.05 

1.90 
2.36 

2.75 
3.33 

3.56 
4.44 
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for  surge  level  increase  and  applies  to  the  final  near-equilibrium 
beach  profiles.  The  value  of  h|^  was  taken  to  be  the  minimum  offshore 
bar  depth  with  respect  to  still  water  surge  level. 

Of  these  three  relationships,  the  values  calculated  using 
^b  ^  h./1.28  are  probably  the  most  accurate  in  view  of  the  large  amount 
of  supporting  data  for  this  relationship,  and  the  fact  that  the  values 
are  determined  at  the  break  point  independently  of  the  shoaling  charac- 
teristics of  the  deep  water  wave. 

The  large  discrepancy  between  the  predicted  values  of  H,  using 
either  equation  (33)  or  (34)  from  the  probable  value  of  H.  found  from 
equation  (32)  undoubtedly  is  a  result  of  the  laboratory  scale  effects, 
more  specifically,  wave  reflection. 

The  distortion  inherent  in  the  modeling  law  gives  offshore  slopes 
twice  as  steep  as  those  found  in  nature.  As  a  consequence,  a  greater 
proportion  of  the  incoming  wave  energy  ended  up  being  reflected  in  the 
wave  tank,  resulting  in  smaller  wave  heights  at  breaking. 

Additionally,  equations  (33)  and  (34)  are  not  preserved  in  the 

model  law,  so  there  exists  some  doubt  as  to  whether  prototype  or  model 

values  of  H  and  L  should  be  used  to  calculate  H,  . 
0     0  b 

Since  the  main  concern  in  this  effort  deals  with  the  area  extend- 
ing from  the  dune  toe  to  just  offshore  of  breaking,  it  was  decided  to 
use  values  of  H,  calculated  from  equation  (32). 

Location  of  Bar  Crest 

One  of  the  more  important  factors  in  constructing  an  equilibrium 
barred  profile  is  the  determination  of  the  bar  location  with  respect 
to  the  dune  toe,  representing  the  limit  of  wave  uprush. 
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Two  previous  papers  touching  upon  this  point  were  those  of  Kemp 
(1960)  and  Sunamura  and  Horikawa  (1974). 

Kemp's  experiments  were  conducted  in  a  three-dimensional  wave  tank 
measuring  15  feet  by  9  feet,  with  a  water  depth  of  6  inches.  For  sand 
ranging  between  0.45  mm  and  2.0  mm,  Kemp  gave  an  empirical  relation- 
ship for  the  distance  between  the  dune  toe  and  the  bar  crest,  X,  ,  as 

be 

24H^^-^/^ 
^bc  "  jl/2   '  ^^^^ 

where  H,  and  X,   are  given  in  feet  and  the  grain  size,  d,  is  in  milli- 
meters. 

Some  doubt  can  be  cast  upon  the  validity  of  this  expression  due 
to  the  very   small  size  of  the  wave  facility  employed  and  the  large 
grain-sizes  used.  Data  collected  in  the  current  dune  erosion  experi- 
ments did  not  match  when  plotted  along  with  Kemp's  results.  However, 
the  difference  is  not  too  great,  indicating  that  the  basic  form  of 
equation  (35)  may  be  nearly  correct  (as  will  be  shown  shortly). 

The  more  recent  work  of  Sunamura  and  Horikawa  (1974)  detailed 
experiments  in  a  25  meter  long  wave  tank  using  a  water  depth  of  1.5 
meters.  Sand  grain-sizes  employed  were  0.2  mm  and  0.7  mm,  and  the 
authors  gave  the  expression  for  bar  crest  location  in  non-dimensional 
form  as 


be  0   n  „,o/  b 


2^  =  0.018(-^)   .  (36) 

^b 

Their  experiments  began  with  a  plane-sloping  beach  which  was  eroded 
by  the  ensuing  wave  action. 
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The  main  objection  to  equation  (36)  is  that  H,  is  given  as  the 
breaker  height  at  the  onset  of  the  erosion  when  the  profile  is  still 
a  plane  slope.  Since  nature  does  not  form  barred  profiles  from  plane- 
sloped  beaches,  equation  (36)  has  no  real  application  outside  of  model 
testing.  To  be  more  practical,  H^^  must  be  the  breaker  height  when  the 
final  equilibrium  profile  is  reached. 

In  addition,  neither  of  the  two  above  mentioned  papers  gives  any 
physical  interpretation  of  the  predictive  equation  for  X, 


^bc- 


Development  of  the  Equation  for  Bar  Crest  Distance,  X, 


The  important  characteristics  which  determine  the  equilibrium 
position  of  the  offshore  bar  formed  by  plunging  breakers  obviously  are 
the  wave  period,  the  breaking  wave  height,  and  the  sediment  grain-size. 
Formulation  of  an  expression  using  the  breaking  wave  characteristics 
rather  than  the  deepwater  wave  parameters  eliminates  any  need  for  con- 
sideration of  the  shoaling  of  the  waves,  the  surge  level,  or  the  pre- 
storm  profile.  Since  it  is  desirable  to  nondimensionalize  the  breaker 
distance,  X,  ,  the  most  logical  choice  is  the  breaker  wave  length,  L,  , 
since  it  incorporates  both  wave  period  and  breaker  wave  height,  and 
the  ratio  X,  /L,  is  preserved  between  model  and  prototype.  Experi- 
mental results  indicate  that  this  "dimensionless  breaker  distance"  is 
not  constant,  but  instead  it  varies  directly  with  H,  and  inversely 
with  wave  period,  T.  This  corresponds  with  recent  results  obtained  by 
Noda  (1978)  and  Saville  (1980),  who  point  out  that  profile  similarity 
can  be  obtained  when  the  dimensionless  fall  velocity  parameter  is  con- 
served, thus  implying  that  H^^/wT  is  indeed  an  important  parameter  for 
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determi nation  of  the  breaker  distance.  It  should  then  be  expected 
that 


where 


X,   =  breaker  distance  (distance  from  bar  crest  to  uprush 
limit), 

L,  =  breaking  wave  length, 

H,  =  breaking  wave  height, 

T  =  wave  period,  and 

0)  =  sediment  fall  velocity. 

Using  the  expression  given  by  Van  Dorn  (1978)  for  the  breaking  wave 

length,  expressed  as 

Lb  -  C^T  =  T[(2g)(0.77Hb)]^/2  ,  (37) 

values  of  \^/^^   were  plotted  in  Figure  26  against  H,/uT  for  the  29 
experiments  where  the  offshore  bar  feature  was  clearly  discernible. 
While  the  data  contain  considerable  scatter,  a  linear  relationship 
appears  to  exist,  which  was  found  to  be 

H  X 

-Y  =   2.691  T^  -  0.068  (38) 

by  the  method  of  least-squares.  The  coefficient  of  determination  for 

2 
this  curve  fit  was  r  =  0.684,  which  is  not  as  good  as  hoped  for,  but 

still  indicates  a  strong  trend. 

The  scatter  of  the  data  can  be  attributed  to  the  difficulty  in 

selecting  the  actual  position  of  the  bar  crest,  which  is  a  gently 
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sloping  crown,  and  also  to  the  possibility  that  a  full  equilibrium 
condition  had  not  been  reached  at  the  time  of  final  data  taking.  None- 
theless, it  is  felt  that  this  relationship  is  basically  correct. 

In  equation  (38)  it  is  seen  that  the  linear  offset  of  0.068  is 
quite  small,  and  that  the  data  can  be  represented  almost  as  well  by 
a  least-squares  fit  through  the  origin  given  as 

H,       X, 

^=2.53^.  (39) 

b 

It  should  be  noted  that  attempts  to  locate  prototype  data  to  fur- 
ther verify  equation  (39)  proved  fruitless.  Data  obtained  from  nature 
typically  do  not  represent  severe  storm  conditions  since  milder  wave 
climates  have  altered  the  bar  features  before  poststorm  measurements 
could  be  taken.   In  addition,  grain-size  and  wave  period  data  are 
generally  lacking  in  the  prototype.  Large  scale  model  tests,  such  as 
those  by  Rector  (1954)  and  Watts  (1954),  were  conducted  using  an 
initially  plane-sloped  beach,  and  the  resulting  barred  profiles  did 
not  exhibit  a  reasonably  well-defined  location  for  the  bar  crest  and 
the  uprush  limit.  Thus,  they  could  not  be  used  either. 

Physical  Interpretation 

A  physical  interpretation  of  equation  (39)  can  be  obtained  by 
rearranging  it  into  the  form 

.X,^  =  0.38-p. 
By  replacing  L,  with 
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Lb  =  ^T 


and  multiplying  both  sides  by  H,  ,  the  equation  becomes 

""b^bc  ^  ^■^^%\   • 

By  noting  that  the  average  wave  energy  flux  (or  wave  power)  per  unit 
width  of  surf  zone  perpendicular  to  the  beach  at  breaking  is  given  by 
linear  wave  theory  as 

P  =  i  %\   '  (40) 

equation  (39)  finally  becomes 

In  equation  (41)  the  product  H^^X^  represents  a  water  volume  per  unit 
width  in  the  surf  zone  which  is  directly  related  to  the  total  volume 
of  water  per  unit  width  contained  in  the  surf  zone.  Thus,  equation 
(41)  indicates  that  the  incoming  breaking  wave  energy  flux  per  unit 
width,  entering  the  surf  zone  across  a  vertical  plane  located  at  the 
bar  crest,  is  directly  proportional  to  the  product  of  the  surf  zone 
volume  and  the  sediment  fall  velocity. 

Considering  that  for  equilibrium  conditions  this  energy  must  be 
totally  dissipated,  primarily  through  turbulence  associated  with  the 
wave  breaking  and  its  propagation  as  a  bore  up  the  beach,  it  can  be 
logically  seen  that  an  increase  in  energy  flux  must  also  increase  the 
breaker  distance  in  order  to  provide  a  larger  surf  zone  volume  to 
dissipate  this  increased  energy  input. 
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It  can  also  be  seen  that  for  a  given  energy  flux,  increasing  the 
grain-size  of  the  sediment  increases  the  fall  velocity,  which  in  turn 
decreases  the  breaker  distance.  This  is  because  the  larger  grains 
take  more  energy  to  be  turbulently  tossed  about,  and  they  settle  more 
rapidly,  resulting  in  greater  energy  dissipation.  As  a  consequence, 
beaches  with  coarser  grains  will  have  a  shorter  breaker  distance,  and 
hence  they  will  be  steeper.  This  fact  has  long  been  known  from  both 
model  and  prototype  observations. 

For  nonequilibrium  conditions,  such  as  were  observed  in  the  current 
dune  erosion  experiments,  the  breaker  distance,  X.  ,  steadily  increased 
until  the  time  when  near-equilibrium  was  established.  Since  more 
energy  flux  per  unit  width  was  entering  the  surf  zone  than  could  be 
dissipated  over  the  shorter  breaker  distance,  there  was  an  "excess  of 
energy"  which  was  expended  through  erosion  of  the  dune  and  the  offshore 
movement  of  sediment.  Eventually  the  equilibrium  breaker  distance  was 
established,  and  erosion  of  the  dune  ceased. 

A  further  interesting  aspect  of  equation  (39)  can  be  seen  by 
rearranging  the  equation  and  replacing  L,  with  equation  (37),  which 
yields 

\c :r-^  •  («> 

Equation  (42)  closely  resembles  the  empirical  equation  proposed  by 
Kemp  (equation  (35)),  with  H,  being  raised  to  the  same  power.  The 
main  difference  lies  in  Kemp's  choice  of  grain-size  diameter  instead 

of  grain-size  fall  velocity.  However,  it  should  be  pointed  out  that 

1/2 
in  the  turbulent  flow  range,  to  is  proportional  to  d  '  ,  so  Kemp's 

equation  is  essentially  of  the  same  form  as  equation  (42). 
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It  should  also  be  pointed  out  that  wave  period  does  not  appear 
explicitly  in  equation  (42);  but  Hj^  is  certainly  a  function  of  wave 
period,  as  seen  in  equations  (33)  and  (34)  which  relate  H,  to  L  . 

A  final  feature  of  equation  (39)  is  the  fact  that  it  is  preserved 
when  scaling  from  the  model  to  the  prototype,  due  to  the  fact  that 
both  Xj^^/Lj^  and  H, /.jT  are  preserved  in  the  model  law.   In  other  words, 
equation  (39)  is  valid  as  it  stands  for  either  the  model  or  the 
prototype. 

In  conclusion,  the  relationship  given  by  equation  (39)  appears  to 

be  a  reasonable  estimate  of  the  breaker  distance,  X,  ,  between  the 

be 

offshore  bar  crest  and  the  limit  of  uprush  on  the  beach.  The  physical 
interpretation  seems  logically  sound,  and  it  is  felt  that  further 
developments  and  measurements  will  bear  this  out.  Finally  it  should 
be  reiterated  that  this  analysis  applies  only  to  plunging  breakers  in 
storm  conditions. 

Location  of  Bar  Trough 

The  location  of  the  bar  trough  was  also  assumed  to  be  related  to 

the  breaking  wave  length  since  the  plunging  distance  of  a  breaking  wave 

determines  the  location  of  the  trough  and  this  distance  must  be  a 

function  of  wave  length.  Denoting  the  distance  between  the  bar  trough 
and  the  dune  toe  as  X,   and  assuming 

—  =  f  (— ) 

29  data  points   representing  near-equilibrium  conditions  were  obtained 
from  the  model   experiments  and  plotted  as   shown  on   Figure  27. 
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Once  again  the  scatter  is  disturbing,  but  it  appears  that  a  linear 
relationship  does  exist.  The  least-squares  best  fit  is  given  by 

H        X 

-Y  =   3.166  j^  +  0.004  ,  (43) 

with  the  coefficient  of  determination  being  r  =  0.722. 

Since  the  linear  offset  is  very  small,  it  is  more  convenient  to 
force  the  best-fit  through  the  origin  yielding 

-^=3.17^^.  (44) 

The  scatter  can  once  again  be  attributed,  somewhat,  to  the  difficulty 
in  selection  of  the  actual  position  of  the  bar  trough. 

As  was  the  case  for  equation  (39),  equation  (44)  is  preserved  in 
the  model  law,  thus  it  can  be  used  for  either  model  or  prototype  pre- 
dictions. 

Dividing  equation  (39)  by  equation  (44)  gives  the  interesting 
result  that 

^bc 

y^=  1.21  .  (45) 

%t 

Or  as  X,   increases  so  does  X,  ,  ,  proportionately.  Figure  28  shows  the 
variation  in  X,  /X,  .  as  recorded  from  the  model  experiments.  The 
average  of  all  points,  not  surprisingly,  is  1.21  with  a  root-mean- 
square  (RMS)  error  of  0.039. 

Building  upon  the  physical  interpretation  given  for  the  breaker 
distance,  X,  ,  it  can  be  seen  that  equation  (44)  and  the  resulting 
equation  (45)  imply  that  an  increase  in  the  wave  energy  flux  per  unit 
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width  entering  the  surf  zone  increases  the  breaker  distance,  X^  , 

be' 

which  in  turn  increases  X,  ,  in  such  a  way  that  the  ratio  of  X,  /X,  , 
remains  constant.  This  can  only  be  done  if  the  breaking  wave  plunging 
distance  increases  with  the  wave  energy  flux  per  unit  width,  as  would 
logically  be  expected.  This  aspect  is  enlarged  upon  in  the  following 
section. 

Plunge  Distance  Hypothesis 

Galvin  (1969)  developed  an  elementary  expression  for  the  plunging 
wave  plunge  distance,  X  ,  given  as 

Xp  ^  2H^  .  (46) 

This  development  assumes  that  at  breaking  the  water  particles  at  the 
crest  are  moving  forward  with  the  speed  of  the  wave,  C,  ,  and  the  par- 
ticles are  moving  vertically  subject  only  to  the  acceleration  of 
gravity.  However,  Galvin' s  experimental  results,  conducted  using  a 
rigid,  plane  beach  slope,  gave  plunge  distances  which  were  as  large 
as  4.5H,  ,  with  an  average  of  about  3H, .  Galvin' s  data  are  shown  plotted 
on  Figure  29a.  As  can  be  seen,  an  attempt  to  sketch  a  straight  line 
through  the  origin  fitting  the  data  meets  with  limited  success. 
Galvin  attributed  this  poor  fit  to  the  simple  assumptions  made  in  the 
"free-fall"  analysis.  As  he  correctly  pointed  out,  there  is  a  sig- 
nificant upward  vertical  particle  velocity  present  as  the  wave  breaks, 
which  further  complicates  matters. 

By  approximating  the  plunge  distance  as  the  distance  between  the 
bar  crest  and  the  bar  trough,  equation  (39)  and  equation  (44),  when 
subtracted,  result  in 
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FIGURE  29:   GALViri'S  BREAKING  WAVE  PLUNGE  DISTANCE  DATA. 
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which  is  the  nondimensional  equation  for  plunge  distance,  and  which  is 
also  "preserved"  between  prototype  and  model.  Equation  (47)  carries 
the  same  physical  interpretation  as  before.  The  plunge  distance  is 
proportional  to  the  wave  energy  flux  per  unit  width  crossing  the  bar 
crest.  Rearranging  equation  (47)  and  substituting  for  L,  from  equa- 
tion (37)  gives 

0.081g^/^H,  ^/^ 

\  ' ^^  •  CS) 

So  instead  of  X  being  proportional  to  H,  ,  it  is  seen  that  X  is  related 

3/2 
to  H,  '  .  Figure  29b  shows  Galvin's  results  replotted  as  a  function 

3/2 
of  H,  '  .   It  is  clearly  seen  that  the  values  more  closely  follow  the 

linear  relationship  as  suggested  by  equation  (48).  The  straight  line 

was  sketched  in  to  illustrate  this  trend. 

The  fact  that  Galvin's  data,  which  were  obtained  on  a  rigid  plane 

3/2 
slope  with  no  bar  formation,  exhibit  a  trend  toward  H,  '      adds  further 

credibility  to  the  physical  interpretations  given  thus  far. 

Pursuing  a  simplified  theoretical  development  in  the  same  fashion 
as  Galvin,  assume  that  the  trajectory  of  a  water  particle  at  breaking 
is  that  of  a  free-falling  body  under  the  force  of  gravity.  However, 
the  horizontal  velocity  of  the  wave  crest  is  not  C,  ,  the  breaking  wave 
celerity,  but  C,  ',  which  is  larger  than  C,  .  This  comes  about  by  vir- 
tue of  the  upward  vertical  accelerations  under  the  crest.  As  a  wave 
breaks  in  a  plunging  fashion,  these  vertical  accelerations  transfer 
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momentum  to  the  upper  portion  of  the  wave  crest,  acting  like  a  "thrust" 
in  the  horizontal  direction.  For  this  simplified  analysis,  it  is 
assumed  that  this  "thrust"  results  only  in  an  increase  of  the  hori- 
zontal velocity  at  the  crest. 

From  the  equations  of  uniformly  accelerated  motion,  and  in  accor- 
dance with  the  sketch  given  in  Figure  30,  the  equation  for  the  vertical 
distance  becomes 

H5  =  ^  gt^  ,  (49) 

while  the  horizontal  motion  is  described  by 

Xp  =  Cj^'t  ,  (50) 

where  t  -   time.  Solving  equation  (49)  for  t,  and  placing  the  result 
into  equation  (50),  yields 


2H,  1/2 
Xp  =  C^'[^]    .  (51) 

Solving  equation  (47)  for  X  ,  equating  the  result  to  equation  (51),  and 


solving  for  C,  '  gives 


1/2 
'b  "       0)  '2'       ^'b   "b 


=  M65  (^,  '   1/2   .  (53, 


1/2 
Replacing  H^    in  equation  (52)  with  the  relationship  obtained  from 

equation  (37) ,  i.e., 

H  1/2  =     '^b 
^     [2(.77)g]^/2  ' 

it  is  seen  that 


■137- 


Plunge  Trajectory 


Bar  Trough 


FIGURE   30:      PLUfJGE   DISTANCE   DEFIUITION   SKEl 


-138- 
C,  '  =  0.037  -2-  .  (53) 

DO) 


This  simply  states  that,  in  the  range  of  plunging  breakers,  the  "wave 
thrust  velocity"  at  breaking  varies  as  the  square  of  the  breaking  wave 
celerity. 

A  further  interpretation  of  equation  (53)  can  be  given  by  assuming 
that  C,  '  is  proportional  to  C,  through  the  simple  relationship  of 

where  K  is  greater  than  1  and  is  somehow  related  to  the  vertical 
water  particle  accelerations.  From  equations  (53)  and  (54)  it  is 
immediately  recognized  that 

0.037C, 
^c  =   .         '  (55) 

which  is  a  dimensionless  function  relating  K  to  the  breaking  wave 
celerity.  Thus,  as  the  wave  celerity  increases,  K  increases,  which 
implies  that  the  vertical  accelerations  responsible  for  the  horizontal 
"thrust  velocity"  also  increase.  The  presence  of  the  sediment  fall 
velocity  has  come  about  from  the  use  of  equation  (47),  and  it  could 
possibly  be  interpreted  as  a  measure  of  bottom  friction  or  as  a  damping 
agent  acting  upon  the  vertical  water  particle  accelerations. 

While  the  above  analysis  is  quite  simplistic  and  is  not  intended 
to  be  an  accurate  physical  description  of  plunging  breakers,  the 
hypothesis  put  forth  seems  to  provide  at  least  a  heuristic  model 
which  could  serve  as  a  basis  for  further  research. 
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Nearshore  Curve  Fit 

Examination  of  the  nearshore  portion  of  the  near-equilibrium 
barred  profiles  clearly  indicates  that  the  profile  is  not  a  straight 
line,  but  instead  a  curve  which  usually  exhibits  the  characteristics 
of  having  an  increasing  slope  as  the  water  depth  decreases.  While  it 
is  impractical  to  determine  every  minor  feature  of  the  nearshore  pro- 
file, caused  by  any  number  of  complex  phenomena,  it  is  possible  to 
derive  a  semi -empirical  relationship  for  this  portion  based  on  some 
rather  simplified  physics. 

Dean  (1977)  gave  three  derivations  for  equilibrium  beach  profiles 
in  the  nearshore  region  between  the  surf  zone  and  the  beach.  He 
assumed  that  the  ratio  of  breaker  height  to  water  depth  is  constant, 
landward  of  the  point  of  first  wave  break,  thus  restricting  the  formu- 
lations to  spilling  breakers.  The  first  model  assumed  the  beach  pro- 
file was  due  to  a  uniform  alongshore  shear  stress.  Since  dune  erosion 
is  considered  here  only  as  an  onshore-offshore  sediment  transport 
problem,  this  model  is  not  applicable.  Dean's  second  model  considered 
that  the  profile  was  a  result  of  uniform  energy  dissipation  per  unit 
surface  area  of  the  sea  bed,  with  the  energy  first  being  transferred 
into  turbulence  and  then,  through  viscous  action,  into  heat.  The 
final  result  was  an  equation  of  the  form 

h  =  Ax2/5  , 
where 

h  =  depth  below  mean  sea  level, 

X  =  horizontal  range  directed  seaward,  and 

A  =  a  function  primarily  of  grain-size. 
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h  =  Ax^/^  .  (56) 

Subsequent  evaluation  of  a  large  number  of  actual  beach  profiles  by 
both  Dean  and  Hughes  (1978)  found  that  equation  (56)  provided  the  best 
fit  to  the  prototype  data  in  a  substantial  majority  of  the  cases,  the 
implication  being  that  the  mechanism  of  beach  profile  formation  within 
the  surf  zone  is  best  described  as  a  uniform  energy  dissipation  per 
unit  volume,  in  the  simplest  of  approximations. 

It  was  seen  earlier  in  this  chapter  that  the  breaker  distance 
(or  surf  zone  length)  can  be  expressed  as  a  function  of  wave  power 
per  unit  width,  which  is  also  proportional  to  h,  . 

Multiplying  the  left-hand  side  of  equation  (41)  by  b/b,  where 
b  =  unit  width  of  surf  zone  perpendicular  to  the  beach,  leads  to  the 
conclusion  that  the  incoming  energy  flux  per  unit  width  is  proportional 
to  some  surf  zone  volume  per  unit  width.  Since  the  same  conclusion 
can  be  reached  for  the  bar  trough  distance,  it  can  be  reasoned  that 
the  energy  present  after  the  initial  plunging  wave  break  is  dissipated 
over  the  volume  of  the  remaining  surf  zone,  as  is  the  case  for  Dean's 
third  derivation. 

After  the  initial  plunge,  which  results  in  a  significant  energy 
loss,  the  remaining  wave  energy  moves  as  a  bore  toward  the  beach.  So, 
from  the  bar  trough  to  the  shoreline,  the  wave  can  be  considered  as  a 
spilling  breaker,  and  Dean's  profile  derivation  may  be  used.  For  com- 
pleteness. Dean's  derivation  is  briefly  outlined  in  Appendix  E.l. 


-141- 
Experimental  Results 

Using  a  direct  power-curve  fitting  technique  as  described  by 
Hughes  (1978),  the  "best-fit"  of  the  equation  h  =  hy}^'^   was  performed 
on  the  nearshore  beach  portion  for  33  of  the  experimental,  near- 
equilibrium  profiles.  The  curve-fitting  extended  from  the  dune  toe 
position  seaward  to  the  offshore  bar  trough,  and  the  method  employed 
the  technique  of  least-squares. 

Figure  31  presents  the  calculated  values  of  "A"  for  each  of  the 
runs,  which  show  only  slight  deviation  from  the  computed  average,  as 
given  by  the  straight  line.  In  fact,  the  root  mean  square  error  about 
the  mean  value  is  only  0.0386. 

The  RMS  errors,  calculated  for  each  profile  by  the  formula 

21  ^/2 


error   [_n  \    ^   measured  "  calculated''  J 


varied  between  0.142  inches  to  0.634  inches  with  an  average  RMS  error 
equal  to  0.336  inches.  This  corresponds  to  0.7  feet  in  the  prototype. 
Errors  of  this  magnitude  are  certainly  tolerable  and  are  much  better 
than  the  errors  found  when  the  profile  is  approximated  by  a  straight 
line  in  the  nearshore  region. 

Discussion 

From  the  results  of  the  curve-fitting,  it  is  plainly  evident  that 
"A"  is  neither  a  function  of  breaking  wave  height,  H,  ,  or  wave  period, 
T.  This  was  also  pointed  out  by  Hughes  when  he  examined  the  values  for 
"A"  taken  from  natural  beach  profiles  representing  three  different  wave 
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energy  conditions.  This  can  be  reasonably  explained  by  the  fact  that 
the  extra  surf  zone  volume,  necessary  to  dissipate  an  increased  incoming 
wave  energy  flux,  is  obtained  by  a  lengthening  of  the  surf  zone.  In 
addition,  as  will  be  seen  shortly,  increased  energy  flux  per  unit  width 
also  increases  the  depth  of  the  bar  crest  and  bar  trough  with  respect 
to  still  water  level  in  such  a  manner  that  the  ratio  of  the  depths 
appears  to  be  constant. 

Thus  for  increasing  wave  energy,  the  position  of  the  bar  trough 
moves  offshore  and  becomes  deeper  in  such  a  manner  that  the  curve 

h  =  Ax2/3 

can  be  extended  seaward  to  intersect  this  position  without  changing 
the  value  of  "A."  The  increased  energy  is  dissipated  by  the  longer 
plunge  distance  and  the  increased  surf  zone  volume. 

However,  it  can  be  stated  that  "A"  is  a  function  of  grain-size 
fall  velocity,  ,:j,  since  a  larger  value  of  co  will  dissipate  more  energy 
and  permit  a  shorter  surf  zone  length.  Because  the  depth  of  the  bar 
trough  will  remain  the  same,  it  will  be  necessary  for  "A"  to  increase 
as  some  function  of  uj.  Thus,  there  seems  to  be  some  sort  of  trade-off 
between  grain-size  and  surf  zone  volume  since  an  increased  co  decreases 
volume  by  decreasing  surf  zone  length,  but  it  also  increases  volume  by 
providing  a  steeper  profile  out  to  the  bar  trough.  It  is  qualitatively 
felt  that  the  net  volume  change  will  bring  about  a  decrease  in  total 
surf  zone  volume  for  increasing  grain-size,  with  the  larger  grain-sizes 
dissipating  more  energy. 

It  is  therefore  concluded  that  the  shape  factor.  A,  in  the  near- 

2/3 
shore  beach  profile  equation,  given  by  h  =  Ax   ,  is  independent  of  the 
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wave  climate  and  is  dependent  only  upon  the  grain-size,  more  specific- 
ally, the  grain-size  fall  velocity.  Furthermore,  "A"  probably  increases 
as  <^  '  as  was  indicated  in  the  derivation  of  Appendix  E.l.  Results 
by  Hughes  also  illustrate  this  trend,  but  only  in  a  general  sense  be- 
cause the  grain-size  range  covered  was  not  sufficient  to  quantify  the 
results. 

Offshore  Bar  Crest  Depth  Versus  Bar  Trough  Depth 

Several  studies  involving  both  model  experiments  and  prototype 
measurements  have  been  conducted  with  regard  to  offshore  bar  geometry. 

In  one  of  the  earliest  experimental  studies,  Keulegan  (1948) 
pointed  out  that  natural  bars  seemed  to  be  somewhat  longer  and  flatter 
than  those  produced  in  model  experiments.  As  already  demonstrated, 
this  is  a  direct  cause  of  the  need  to  use  model  sediments  similar  in 
size  to  those  found  in  nature,  making  it  necessary  to  use  a  distorted 
model  when  the  parameter  H/toT  cannot  be  scaled.  Keulegan,  however, 
did  point  out  that  the  ratio  of  bar  crest  depth  (hj^^)  to  bar  trough 
depth  (h.  ^)  was  similar  in  both  model  and  prototype,  and  it  appeared 
to  be  practically  independent  of  the  wave  characteristics  or  the  beach 
slope.  He  computed  an  average  value  of 

r^=   1.69 
^c 

for  his  model  experiments. 

Shepard  (1950)  examined  this  ratio  using  a  large  number  of  beach 
profiles  along  the  California  coast.  While  the  scatter  is  signifi- 
cant, the  data  do  indicate  a  trend,  given  by  Shepard  as 
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A  study  of  beach  profiles  from  three  locations  along  the  Texas  coast  by 
Herbich  (1970)  gave  values  of  h.  ^/h.   ranging  between  1.25  and  1.35. 

In  the  current  dune  erosion  experiments,  the  values  ranged  between 
1.15  and  1.97,  with  the  average  at 

^=   1-52  . 

The  29  values  obtained  are  shown  plotted  on  Figure  32,  along  with  the 
averages  given  by  the  other  investigators. 

While  the  relationship  seems  to  be  somewhere  between  1.2  and  1.7, 
it  is  felt  that  too  much  variation  exists  for  this  to  be  used  for  pre- 
diction of  bar  trough  depth.  It  is  interesting  to  note  that  the 
laboratory  tests  give  higher  values  for  h,  ./h,   than  those  found  in 
nature.  This  can  be  primarily  attributed  to  the  fact  that  the  labora- 
tory tests  were  conducted  using  monochromatic  waves,  while  the  wave 
trains  in  nature  contain  a  spectrum  of  waves  which  act  to  "smooth"  the 
bar  features  due  to  the  fact  that  the  different  height  waves  will  break 
in  different  locations  above  the  bar.  A  secondary  influence  is  the 
tidal  variations  found  in  the  prototype,  which  also  help  to  smooth  the 
bar  features. 

It  can  be  concluded  that  the  ratio  of  bar  trough  depth  to  bar 
crest  depth  is  not  necessarily  a  constant  because  of  the  considerable 
variation  shown  on  Figure  32.  This  ratio  is  more  than  likely  a  func- 
tion of  the  plunging  wave  and  of  the  huge  turbulent  vortex  it  generates, 
which  scours  out  the  bar  trough. 
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FIGURE  32:  BAR  TROUGH  DEPTH  VERSUS  BAR  CREST  DEPTH. 
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Wave  Runup 

Since  the  nearshore  profile  is  calculated  referenced  to  the  dune 
toe,  which  corresponds  to  the  point  of  maximum  wave  runup,  it  is  neces- 
sary to  be  able  to  predict  the  elevation  of  the  maximum  runup  above  the 
still  surge  level  in  order  to  construct  the  equilibrium  barred  profile. 
The  vertical  runup  elevation  was  selected  over  the  horizontal  runup 
distance  because  the  experiments  demonstrated  much  more  consistency 
than  when  the  horizontal  distance  was  chosen.  This  was  caused  by  the 
gentle  beach  slopes  which  necessitated  a  linear  interpolation  between 
profile  data  points  to  determine  the  intersection  of  still  surge  level 
and  the  beachface.  Thus,  small  errors  in  elevation  could  cause  sig- 
nificant errors  in  the  horizontal  position  of  this  intersection  point, 
and  consequently,  give  errant  values  for  the  horizontal  runup  distance. 

Most  previous  experimental  work  dealing  with  wave  runup  was  con- 
ducted using  plane,  rigid  beach  slopes.  A  simple,  reliable  expression 
for  runup  on  smooth  plane-slopes  was  given  by  Hunt  (1959)  and  can  be 
written  as 


V  _  tana  [cyy. 


where 


R  =  vertical  runup  distance  above  still  water  level, 
tana  =  beach  slope, 
H  =  wave  height,  and 
L  =  deep  water  wave  length. 
Unfortunately,  equation  (57)  is  not  valid  for  cases  other  than  plane- 
sloped  beaches.  As  has  been  seen,  actual  beaches  have  curving  profiles 
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which  get  progressively  steeper  nearer  to  the  beach.  In  addition,  the 
rigid  bottom  must  alter  the  wave  characteristics  from  those  present 
on  natural,  barred  profiles. 

A  method  for  calculating  wave  runup  on  composite  slopes  was  given 
by  Saville  (1957b).  While  perhaps  better  than  Hunt's  formula,  the 
method  requires  an  iteration  procedure,  and  unfortunately  it  still 
deals  only  with  rigid  beds. 

In  the  present  experiments,  the  values  of  vertical  wave  runup  with 
respect  to  still  surge  level  were  calculated  from  the  profile  measure- 
ments and  then  examined  for  any  significant  trend.  No  apparent  rela- 
tionship could  be  seen  as  a  result  of  variations  in  the  breaking  wave 
height,  which  is  related  to  energy  flux.  However,  a  definite  trend  was 
seen  to  exist  with  regard  to  wave  period.  As  wave  period  increased, 
the  vertical  runup  increased,  which  is  a  pretty  common  observation. 

It  also  seems  logical  that  R  should  somehow  be  directly  related 
to  beach  slope,  as  it  is  in  equation  (57). 

Rearranging  Hunt's  formula  into  the  form 


R^  -  (HLQ)^/2tana 


and  replacing  L  with 

=      0 


yields 


V         2iT 
For  gentle  beach   slopes. 


0  ZiT 


nH    1/2 

R     =  [?1]         T  tana  .  (58) 
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tana  =  f 


but  the  profile  has  been  shown  to  be  h  =  Ax  '    ,   from  which  it  is  easiV 
shown  that 


dh  ^  2  A^ 
dx  ~  3  ^1/2 


tana  .  (59; 


Placement  of  equation  (59)  into  equation  (58)  yields 

1/2  u  1/2 


\  -  !  t4f^  ^l^         Ta3/2  .  (60: 


But  within  the  surf  zone,  it  is  assumed  that  H  =  Kh,  where  K  =   0.78,  and 
equation  (70)  becomes 

R  =  0.235g^/^TA-^/^  ,  (61) 

which  is  independent  of  the  wave  height.  Figure  33  gives  the  runup 
calculated  when  placing  the  value  of  "A"  found  for  each  experiment  into 
equation  (61),  as  compared  to  the  actual  runup  value.  As  can  be  seen, 
equation  (61)  typically  gives  larger  values  for  vertical  runup,  which 
indicates  that  either  the  assumed  value  for  K  is  too  large  or,  more 
likely,  effects  such  as  bottom  friction  and  percolation  on  the  beach- 
face  aid  in  reducing,  slightly,  the  value  of  R  . 
Writing  equation  (61)  as 


R  =  TA^^  X  constant 

V 

2/3 
and  noting  that  A  is  proportional  to  ij  '  ,  as  shown  in  Appendix  E.l, 

the  final  expression  for  vertical  runup  on  a  natural  beach  slope  can 
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FIGURE  33:  CALCULATED  RUNUP  VERSUS  MEASURED  RUNUP. 
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be  given  as 


R 

"Y  =  constant  .  (62^ 


This  runup  parameter  has  the  features  of  being  dimensionless  and  also 
of  being  preserved  in  the  model  law.  In  other  words,  the  value  of  the 
constant  will  be  the  same  in  both  the  prototype  and  the  model. 

Physically,  equation  (62)  states  that  the  vertical  runup  distance 
is  influenced  only  by  the  wave  period  for  a  given  beach.  Therefore, 
period  must  be  the  only  characteristic  wave  parameter  which  is  carried 
all  the  way  across  the  surf  zone.  The  wave  height  characteristics  are 
lost  in  the  process  of  forming  and  maintaining  the  equilibrium  nearshore 
profile. 

It  is  also  seen  that  the  sediment  fall  velocity  influences  the 
runup,  since  runup  increases  as  w  increases.  This  could  be  inter- 
preted in  the  following  manner.  As  the  grain-size  increases,  oj  in- 
creases and  the  beach  becomes  steeper,  as  was  seen  earlier  in  this 
chapter.  This  steeper  beach  profile  must  dissipate  slightly  less  wave 
energy  through  turbulence,  which  means  slightly  more  energy  must  be 
dissipated  in  the  final  wave  runup,  thus  giving  larger  values  for 

R  . 
v 

For  the  dune  erosion  experimental  series,  the  values  of  R  /uT 
have  been  calculated  and  are  shown  on  Figure  34.  A  certain  amount 
of  scatter  exists,  but  not  enough  to  distract  from  the  apparent 
trend  of  R  /luT  =  constant.  Three  factors  could  contribute  to  the 
scatter: 

1)  Measurement  error,  which  probably  does  not  amount  to  very  much 
since  the  profiling  instrument  calibration  is  quite  good. 
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However,  the  manual  placement  of  the  probe  could  create  errors 
on  the  order  of  ±0.1  inch.  This  error  translates  to  approxi- 
mately ±0.13  in  the  value  of  the  constant,  1.9. 

2)  It  is  very  likely  that  the  final  profile,  while  near  equilib- 
rium, had  not  reached  a  final  equilibrium  in  some  instances, 
thus  contributing  to  the  error. 

3)  While  the  profile  remained  quite  uniform  across  the  width  of 
the  wave  tank,  there  sometimes  appeared  to  be  minor  variations 
in  the  profile  as  a  result  of  the  side  wall  effects. 

Taking  these  effects  into  account,  it  would  appear  that  R  /uT 
does  remain  nearly  constant. 

The  average  for  all  data  points  shown  on  Figure  34,  excluding  the 
three  circled  points  which  fall  way  below  the  trend,  is  given  by 

^=1.9.  (63) 

The  calculated  RMS  error  was  0.267,  which  is  about  14  percent. 

An  examination  of  Ma's  (1979)  results  for  different  grain-sizes 
gave  values  for  R  /jT  in  the  neighborhood  of  1.9,  but  they  varied  con- 
siderably, due  mostly  to  the  difficulty  of  determining  R  from  the 
profiles. 

Due  to  the  above  physical  agreements  and  experimental  results,  it 
is  concluded  that  the  dimensionless  parameter  R  /coT  is  a  constant,  and 
represents  a  simple  means  for  the  determination  of  vertical  wave  runup 
for  either  prototype  beaches  or  movable-bed  models. 
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Off shore  Curve  Fit 

In  the  region  seaward  of  the  offshore  bar  the  bottom  equilibrium 
profile,  arising  from  storm  activity,  exhibits  features  similar  to  the 
nearshore  profile.  Comparisons  of  the  near-equilibrium  profiles  for 
the  experimental  runs  indicate  two  obvious  trends  which  can  aid  in 
the  analysis  of  this  portion  of  the  barred  profile,  the  profile  becomes 
progressively  steeper  as  the  depth  decreases,  and  the  overall  profile 
steepens  as  the  peak  surge  level  increases.  From  the  first  of  these 
observations  comes  the  fact  that  the  profile  can  best  be  described  as 
a  power  curve  rather  than  as  a  straight  line.  This  corresponds  with 
the  same  conclusion  found  for  the  nearshore  region.  The  second  obser- 
vation makes  it  necessary  to  incorporate  the  surge  level  into  the 
analysis. 

Making  use  of  some  known  physical  relationships  for  sediment  trans- 
port in  this  region,  combined  with  a  few  simplifying  assumptions,  it 
was  possible  to  derive  a  simple  power  curve  relationship  for  the  pro- 
file in  the  immediate  offshore  region,  given  by  the  expression 

h  .  A^(x  -  x.)2/^  ,  (64) 


where 


h  =  depth  relative  to  peak  surge  level  (feet), 

X  =  horizontal  range  (feet), 

x.  =  reference  range  where  the  curve  intersects  the  water 
level  (feet),  and 


A  =  shape  factor  (feet^^^). 


0 

The  complete  derivation  is  presented  in  Appendix  E.2  along  with  the 
assumptions  and  reasoning  implemented.  Equation  (54)  is  of  the  same 
form  as  equation  (56),  which  was  used  for  the  nearshore  profile.  This 
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implies  a  similarity  of  profiles  between  these  two  regions,  a  fact 
that  intuitively  seems  correct. 

Experimental  Results 

Using  the  method  of  least-squares  curve  fitting,  the  power  curve 
of  equation  (64)  was  applied  to  the  experiment  profiles  in  order  to 
find  the  values  for  A  .  As  can  be  seen  from  Figure  35a,  the  values  of 
A  appear  to  be  independent  of  wave  height  and  wave  period,  but  they 
are  definitely  a  function  of  surge  level.  While  not  proven  in  this 
research,  it  is  also  felt  that  A  will  have  some  dependence  on  grain- 
size,  as  was  the  case  for  the  nearshore  profile.  The  profile  should 
become  steeper  for  an  increased  value  of  sediment  fall  velocity. 

The  values  of  vertical  RMS  error  calculated  for  each  profile 
varied  between  0.136  inches  to  0.583  inches,  with  the  average  RMS  error 
for  all  profiles  being  0.304  inches.  This  average  corresponds  to  a 
value  of  7.6  inches  in  the  prototype. 

Figure  35b  illustrates  the  linear  trend  given  by  the  averages  of 
all  values  at  each  surge  level,  at  least  for  the  surge  range  tested. 
This  linear  trend  is  given  by  the  best-fit  of 

A„  =  0.167  +  0.058S  , 
0  m 

which  leads  to  the  offshore  profile  equation  for  the  model  given  as 


h^  =  (0.167  +  0.058S^)(x^  -  x.)2/3  ,  (55) 

where  S^  is  the  surge  level  in  inches,  and  h  ,  x  ,  and  x.  are  also  in 
m                           m   m      1 

inches.  The  equivalent  equation  for  the  prototype  is  scaled  up  as 
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hp  =  (0.132  +  0.022Sp)(Xp  -  x-)^/^  .  (66) 


Here 


h  =  depth  in  feet, 
x^  =  range  in  feet, 
X.  =  reference  range  in  feet,  and 
S  =  surge  level  in  feet. 
Unfortunately,  equations  (65)  and  (65)  have  some  rather  severe 
limitations  placed  on  their  usage.  These  are  discussed  below. 

Applications 

There  are  two  major  constraints  placed  upon  the  empirically  de- 
rived offshore  profile  equations  given  above. 

Because  the  experimental  test  series  did  not  attempt  to  investi- 
gate the  effects  of  different  prestorm  offshore  profiles,  the  results 
cannot  be  universally  applied.  During  the  storm  surge  rise,  sediment 
is  being  transported  offshore  and  finally  deposited  when  the  bottom 
velocities  are   sufficiently  decreased.  As  the  surge  rises,  some  sedi- 
ment deposited  at  lower  surge  levels  can  no  longer  be  moved  because  the 
limiting-depth  of  transport  has  been  exceeded.  So,  in  essence,  the 
process  of  profile  formation  is  actually  a  series  of  discrete  steps 
which  take  place  as  the  water  level  rises.  The  portion  which  developed 
at  lower  surge  levels  is  incremented  as  the  surge  level  increases,  thus 
making  the  profile  steeper.  This  explains  why  the  values  of  A  appear 
to  be  related  to  surge  level. 

In.  the  same  manner  it  can  be  seen  that  an  initially  steeper  pre- 
storm profile  will  give  a  higher  value  of  A  ,  indicating  a  steeper 
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offshore  storm  profile.  This  comes  about  because  the  eroded  sediment 
cannot  be  carried  as  far  offshore  before  it  is  dropped,  when  compared 
to  the  milder  prestorm  profile. 

The  second  constraint  is  the  prototype  grain-size.  Because  equa- 
tion (66)  was  scaled  up  using  the  modeling  relationships  derived 
earlier,  it  can  only  be  applied  to  profiles  having  a  similar  grain- 
size  of  around  0.26  mm.  This  ties  into  the  first  constraint  since  the 
initial  prestorm  profile  steepness  will  be  a  function  of  sediment  fall 
velocity.  In  addition,  it  should  be  expected  that  the  increased 
stability  of  larger  grain-sizes,  when  exposed  to  erosive  conditions, 
should  allow  a  steeper  equilibrium  offshore  storm  profile. 

On  the  brighter  side,  however,  the  experimental  test  series  were 
designed  to  approximate  the  average  Florida  beach  in  terms  of  grain- 
size  and  prestorm  profile;  thus,  the  results  from  the  offshore  curve 
fitting,  given  by  equation  (66),  could  be  used  with  a  higher  degree  of 
confidence  when  applied  to  the  Florida  coastline. 

To  reiterate,  the  empirically  derived  prototype  relationship  for 
the  offshore  portion  of  the  equilibrium  barred  storm  profile,  given  by 
equation  (66),  should  only  be  applied  to  beaches  exhibiting  the  same 
prestorm  underwater  configuration  as  that  given  by  the  expression 

h  =  0.15x^/^  , 

where  0.15  has  the  units  of  (feet)^  ,  which  is  typical  of  Florida 
beaches.  The  prototype  effective  grain-size  diameter  should  be  in  the 
neighborhood  of  0.26  mm,  and  it  is  doubtful  that  the  relationship  could 
be  extrapolated  much  past  a  peak  surge  value  of  16  feet  in  the  proto- 
type. 
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Summa ry 

This  chapter  has  been  concerned  with  the  parameterization  of  cer- 
tain characteristics  of  the  equilibrium  barred  storm  profile.  By 
determining  this  profile,  it  then  becomes  possible  to  determine  the 
maximum  amount  of  dune  erosion  and  recession  which  would  be  expected 
on  any  prestorm  beach-dune  configuration  for  a  given  set  of  storm 
parameters.  This  is  done  simply  by  shifting  the  calculated  storm  pro- 
file relative  to  the  prestorm  profile  until  a  balance  is  established 
between  eroded  and  deposited  sediment. 

It  was  demonstrated  that  the  distance  between  the  dune  toe  and  bar 
crest,  as  well  as  the  distance  between  dune  toe  and  bar  trough,  can  be 
found  as  a  function  of  breaking  wave  height,  wave  period,  and  sediment 
fall  velocity.  Physical  explanations  were  given  for  these  expressions. 

A  hypothesis  was  developed  for  plunging  distance,  but  only  to 
demonstrate  a  possible  direction  for  further  research. 

The  nearshore  profile  between  the  dune  toe  and  the  bar  trough  was 
found  to  fit  a  curve  given  by 

in  the  model,  where  h^  and  x  are  in  inches.  The  equivalent  prototype 
equation  for  the  nearshore  region  scales  up  to  be 


hp  =  0.39Xp2/3  ,  (58) 

where  h  and  x  are  given  in  units  of  feet.  Equation  (68),  however,  is 

restricted  to  use  on  beaches  having  an  effective  grain-size  diameter 
in  the  vicinity  of  0.26  mm. 
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An  investigation  Into  the  relationship  between  the  water  depths 
over  the  bar  and  crest  proved  inconclusive. 

A  simple  means  of  calculating  wave  runup  was  derived  and  empiri- 
cally calibrated  for  use  in  either  the  model  or  the  prototype.  It 
appears  that  the  vertical  runup  is  principally  a  function  only  of  wave 
period  and  sediment  size. 

Finally,  an  analysis  of  the  offshore  portion  of  the  barred  pro- 
file was  undertaken,  which  yielded  results  applicable  to  the  average 
Florida  beach  but  needs  to  be  further  analyzed  to  determine  the  effects 
of  grain  size  and  of  a  different  prestorm  profile. 

From  the  foregoing  analyses,  it  is  now  possible  to  construct  the 
equilibrium  barred  profile  from  the  incoming  wave  parameters. 


Sample  Prototype  Calculation 

Given  the  deep  water  wave  height,  H  ,  the  wave  period,  T,  the 
peak  surge  level,  S,  and  the  sediment  fall  velocity,  w,  the  determina- 
tion of  the  equilibrium  barred  storm  profile  proceeds  as  follows: 

1)  Calculate  H,  from  equation  (34). 

2)  Determine  L,  from  equation  (37). 

3)  Calculate  h,  from  equation  (32). 

4)  Calculate  the  vertical  runup,  R   from  equation  (63). 

5)  Using  equation  (68),  replace  h  with  the  value  for  R  and  solve 
for  X  . ,  This  is  the  distance  between  the  dune  toe  and  the 
intersection  of  the  still  surge  level  with  the  beach.  Using 
this  intersection  point  as  the  zero  reference,  plot  the  posi- 
tion of  the  dune  toe  as  (-x  •,R  ). 

^  pi '  v' 
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6)  Determine  the  distance  between  the  dune  toe  and  bar  trough 
using  equation  (44).  The  horizontal  range  of  the  bar  trough, 
relative  to  the  reference  point,  is  given  as  (x,  ,  -  x  •). 

7)  Using  equation  (58),  calculate  the  nearshore  profile  pairs 
(x  ,h  )  and  reference  them  to  the  chosen  coordinate  system 
through  the  transformations 


X  =  X 


„  -  x„.  and  y  =  R  -  h  . 
p    pi      ■^    V    p 


8)  The  shearface  of  the  eroded  dune  is  drawn  in  as  a  nearly 
vertical  line. 

9)  Calculate  the  bar  crest  distance  from  the  dune  toe  (x,  )  using 
equation  (39);  then  plot  the  bar  crest,  relative  to  the  co- 
ordinate system  as  (X|^^  -  x  .  ,-h|^).  Connect  the  bar  crest 
and  trough  with  a  straight  line. 

10)  Using  equation  (66)  for  the  offshore  portion,  replace  h  with 
the  value  for  h^^  and  solve  for  (x^^  -  x  •  -  x^-)  from  which  the 
value  of  X.  can  be  found.  Placing  this  value  for  x-  back 
into  equation  (66)  then  allows  the  calculation  of  the  offshore 
depths  relative  to  the  intersection  of  surge  level  and  pro- 
file. 
By  the  above  method,  the  equilibrium  barred  profile  can  be  cal- 
culated, and  ultimately  the  maximum  dune  recession  can  be  found  as 
discussed  earlier. 

This  method  was  used  to  calculate  the  profiles  for  runs  18  and  34 
of  the  experimental  series,  chosen  arbitrarily,  and  the  comparisons 
with  the  actual  measured  profiles  are  shown  in  Figure  36.  As  can  be 
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seen,  the  comparisons  are  very  favorable.  It  should  be  noted  that  the 
model-value  equations  (65)  and  (67)  were  used  in  this  case. 

The  next  step  is  to  determine  the  effect  of  peak  surge  duration, 
which  is  examined  in  the  following  chapter. 


CHAPTER  IX 
TRANSIENT  EROSION  EFFECTS 

The  quantitative  analyses  presented  in  Chapter  VIII  provide  a 
means  of  calculating  the  near-equilibrium  erosion  which  would  occur  for 
a  given  set  of  storm  conditions  acting  on  a  given  prestorm  profile. 
However,  the  resulting  erosion  values  represent  the  maximum  dune  re- 
cession which  would  occur  if  the  storm  conditions  persisted  indefi- 
nitely at  the  peak  surge  level.  While  a  severe  "northeaster"  lasting 
three  days  might  well  approximate  this  condition,  hurricanes  move  much 
faster,  with  the  surge  water  level  remaining  at  its  peak  value  often 
only  several  hours.  As  a  consequence,  the  erosion  experienced  by  the 
beach-dune  system  during  a  hurricane,  while  severe,  never  reaches  a 
full  equilibrium.  Hence,  there  remains  the  need  of  being  able  to 
estimate  the  effects  of  the  peak  surge  duration  and  then  to  adjust 
the  calculated  dune  recession  accordingly. 

Essentially,  what  is  required  in  determining  the  transient  erosion 
is  a  partitioning,  with  respect  to  time,  of  the  dune  erosion  which 
occurs  from  the  time  that  the  peak  surge  level  is  reached  until  the 
time  that  near-equilibrium  is  established.  By  doing  this,  it  then 
becomes  possible  to  determine  the  effect  of  peak  surge  duration,  since 
no  further  dune  recession  can  take  place  once  the  surge  level  begins 
to  drop. 
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A  further  advantage  of  developing  a  transient  dune  erosion  model 
is  that  future  refinements  are  possible  as  more  prototype  data  become 
available.  Verification  of  an  "equilibrium-type"  model,  such  as 
Edelman's,  requires  equilibrium  prototype  storm  erosion,  which  seldom 
occurs.  On  the  other  hand,  a  transient  model  can  be  adjusted  to  suit 
any  peak  storm  surge  duration  with  the  advantage  that,  every  time 
storm  data  are  collected,  further  verification  and  refinement  of  the 
model  will  be  possible. 

Theoretical  Development 

As  was  the  case  for  the  development  of  the  equilibrium  barred 
storm  profile,  the  actual  physics  involved  in  the  transient  dune  ero- 
sion phenomenon  are  beyond  our  grasp.  Therefore,  the  investigation 
must  become  a  parameterization  of  the  important  factors,  strengthened 
by  physical  reasoning  where  possible. 

Energy  Concepts 

During  normal  conditions,  a  beach  is  able  to  reach  a  state  of 
quasi-equilibrium  when  subjected  to  a  nearly  constant  incoming  wave 
field.  In  addition,  the  beach  is  able  to  rapidly  adjust  to  small  varia- 
tions in  the  incoming  waves  and  the  tidal  fluctuations  in  water  level. 
In  this  case  the  incoming  wave  energy  flux  is  totally  dissipated  over 
the  length  of  the  surf  zone  by  the  mechanisms  of  wave  breaking,  bottom 
friction,  wave  transformation,  turbulence,  and  to-and-fro  sediment 
transport,  with  the  final  small  amount  of  wave  energy  being  converted 
to  wave  uprush.  Essentially,  there  is  no  remaining  wave  energy  left 
over  at  the  uprush  limit. 
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However,  during  a  severe  storm,  the  water  level  and  the  incoming 
wave  energy  increase  dramatically  over  a  short  period  of  time,  so 
rapidly  in  fact  that  the  beach  is  unable  to  adjust  quickly  enough  to 
totally  dissipate  this  energy  increase.  The  result  is  an  "excess  of 
energy"  left  over  at  the  uprush  limit  which  can  only  be  dissipated  by 
further  wave  uprush  and  the  corresponding  erosion  of  the  dune  sand, 
now  exposed  to  the  wave  action  by  virtue  of  the  surge  level  rise.  In 
time,  the  beach  profile  will  be  eroded  to  the  point  that  equilibrium 
will  once  again  be  reached,  and  the  incoming  wave  energy  will  be 
totally  dissipated  over  the  length  of  the  new  surf  zone,  with  no  ex- 
cess energy  reaching  the  eroded  dune  toe.  In  actual  fact,  this  time 
never  arrives  because  the  storm  event  ends  before  equilibrium  can  be 
reached. 

Based  upon  the  above  description,  it  is  possible  to  write  a 
mathematical  statement,  given  as 


dE  (t)       dE  (t) 


where 


E|^  =  average  wave  energy  density  per  unit  width  at  the  breaking, 

E^  =  average  wave  energy  density  per  unit  width  arriving  at 
the  dune  toe,  and 

B  =  a  time  dependent  proportaional ity. 

Equation  (69)  states  that,  during  a  storm,  the  time-rate  of  wave 

energy  per  unit  width  reaching  the  dune  toe  is  proportional  to  the 

time-rate  of  wave  energy  per  unit  width  reaching  the  breaking  point. 

The  quantities  dE/dt  have  units  of  wave  power  per  unit  width,  while 

6(t)  is  dimensionless. 
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Since  the  "excess  energy"  is  such  a  difficult  parameter  to  work 
with,  it  is  necessary  to  recast  it  into  a  more  useful  form.  This  is 
done  by  reasoning  that  energy  is  a  form  of  work,  and  that  a  certain 
amount  of  work  will  be  required  in  order  to  erode  a  certain  amount  of 
sediment.  Under  the  assumption  that  the  time-rate  of  erosion  (per 
unit  width  of  beach)  occurring  at  the  dune  toe  is  directly  proportional 
to  the  excess  energy  available  at  the  dune  toe,  it  is  possible  to 
write 

dV     dE 

dT^^dT'  (70) 

where 

V^  =  eroded  volume  of  sediment  per  unit  width  of  beach,  and 
K  =  proportionality  constant  having  units  of  inverse  stress. 
It  is  reasonably  expected  that  the  rate  of  volumetric  erosion  should  be 
a  function  of  the  shear  stress  acting  on  the  sand  grains.  Thus,  the 
factor  K  should  be  a  function  of  the  sediment  characteristics,  such  as 
grain  size,  grain  shape,  and  specific  weight. 
Combining  equations  (69)  and  (70)  yields 

^^e       ^^b 

d/^K3(t)^.  (71) 

To  find  a  form  for  6(t),  consider  the  case  of  an  instantaneous  sea 
level  rise  with  a  constant  incoming  wave  energy.  For  this  case,  it  has 
been  well  observed  in  model  experiments  that  the  rate  of  erosion  is 
initially  large  and  then  decreases  in  time,  as  shown  in  Figure  37a. 
This  curve  exhibits  the  characteristics  of  an  exponential  decay. 
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Assuming  the  exponential  form  of 

Kgft)  =  A  +  Be""^ 


and  applying  the  boundary  condition  of  dV  /dt  approaching  zero  as  t 


approaches  infinity,  the  result  is 


6(t)=^.  (72) 


By  placing  equation  (72)  back  into  equation  (59),  and  by  then  assuming 
that,  at  time  t  =  0,  the  rate  of  energy  arriving  at  the  beach  is  a  frac- 
tion, n,  of  the  incoming  energy,  i.e., 

dE^(O)     dE^(O) 
~dt    '^       dt   ' 

it  is  then  seen  that  B  =  nK  and 

&(t)  =  ne-^t  .  (73) 

Finally,  equation  (73)  can  be  placed  into  equation  (71)  to  yield 

dV^(t)       .  dE, 

—dT-  -  ^^^    dT  '  (74) 

which  is  a  parametric  equation  for  the  rate  of  dune  erosion  at  a  fixed 
surge  level . 

Continuing  the  special  case  of  instantaneous  sea  level  rise  and 
constant  incoming  wave  energy,  equation  (74)  can  be  integrated,  if  K, 
n,  and  a  are  not  functions  of  time.  This  integration  gives 

V^(t)  =  ^  (1  -  e-^)  -^  ,  (75) 
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which  has  the  general  form  as  illustrated  in  Figure  37b.  As  shown, 
the  erosion  is  initially  very   rapid,  and  then  it  slows  down  as 
equilibrium  is  approached.  Finally,  the  cumulative  erosion  curve 
approaches  the  equilibrium  erosion  value  asymptotically,  as  would  be 
expected. 

While  the'form  of  equation  (75)  is  intuitively  correct  for  the 
case  of  an  instantaneous  surge  increase  with  constant  incoming  wave 
energy,  the  equation  cannot  be  used  to  describe  an  actual  event  where 
wave  energy  varies  with  time. 

Empirical  Coefficients 

The  factor,  K,  will  be  a  function  of  sediment  characteristics 
alone,  i.e.. 


K-U6^,y^,C,.. 


where 


d^  =  effective  grain  diameter, 
e  ^  ' 

Y  =  sediment  specific  weight,  and 

C  ,  =  grain  drag  coefficient. 
The  proportion  of  breaking  wave  energy  to  initially  reach  the  dune  toe, 
n,  has  a  value  between  zero  and  one,  and  it  is  probably  a  function  of 
surge  level,  wave  period,  wave  height,  and  beach-dune  configuration, 
i.e., 

n  =  n(S,T,H,, profile)  . 
Likewise,  it  should  perhaps  be  expected  that 
a  =  a(S,T,H,, profile)  . 
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Complete  Equation 

In  an  actual  storm,  the  surge  level,  wave  height,  and  wave 
period  are  all  functions  of  time,  which  make  n,  a,  and  dE,/dt  functions 
of  time.  Thus,  the  differential  equation  (74)  becomes 

dVJt)         .,,  ,  dE,  (t) 

-^  =  KMe-^(t)   t  _^  ^  (76) 

which  will  be  most  difficult  to  integrate  until  more  is  known  about 
the  storm  parameter  time-variations,  and  how  they  effect  n(t)  and  a(t). 
If  the  erosion  process  could  be  considered  as  passing  through  dis- 
crete steps  of  equilibrium  as  the  surge  level  rises,  it  might  be 
possible  to  find  functions  for  n  and  a  in  order  to  complete  a  numerical 
integration  of  equation  (76).  However,  comparisons  of  the  experimental 
test  results  strongly  indicate  that  this  is  not  the  case.  So,  for  the 
present,  it  will  not  be  possible  to  predict  the  transient  erosion 
taking  place  while  the  surge  level  rise  is  on  the  increase. 

Discussion 

The  development  just  completed  cannot  be  put  into  use  at  this 
time  due  to  the  difficulties  in  determining  the  empirical  coefficients. 
The  only  portion  of  the  present  experiments  where  portions  of  the 
analysis  might  hold  would  be  from  the  time  of  peak  surge  until  the 
time  of  near-equilibrium.  In  this  case,  the  arguments  used  in  arriving 
at  the  integrated  equation  (75)  apply  to  the  process.  Given  enough 
data,  it  would  be  possible  to  perform  a  curve-fit  in  order  to  determine 
T]   and  a.  However,  most  of  the  experiments  provided  only  3  data  points, 
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which  make  curve-fitting  somewhat  unreliable.  Attempting  to  record 
more  data  would  not  really  help  the  cause  since  the  erosion  occurs  in 
"chunks"  as  large  sections  of  the  dune  break  away  and  collapse  into 
the  water.  Thus,  the  actual  process  is  not  strictly  continuous. 
Another  factor  to  be  considered  is  the  continuous  starting  and  stopping 
of  the  wave  paddle  in  order  to  make  measurements.  This  must  certainly 
have  an  adverse  effect,  if  done  too  often. 

Physically,  the  preceding  development  seems  to  fit  in  with  the 
results  given  in  Chapter  VIII.  In  that  case,  it  was  seen  that  the  surf 
zone  length  was  determined  by  the  incoming  wave  power.  For  the  equilib- 
rium case,  this  power  was  dissipated  over  the  length  of  the  surf  zone. 
The  concept  of  "excess  energy"  implies  that  the  surf  zone  is  not  of 
sufficient  length  necessary  to  dissipate  all  of  the  incoming  wave 
power.  Thus,  during  the  transient  erosion,  it  would  be  expected  that 
the  distance  between  the  bar  crest  and  the  dune  toe  would  be  less  than 
the  final,  equilibrium  distance.  This  trend  can  be  observed  for  all 
the  time-dependent  experimental  profiles  shown  in  Appendix  C. 

Thus,  the  theoretical  development  presented  has  physical  roots, 
and  it  might  possibly  lend  itself  to  further  development  in  the  future. 
At  present,  though,  it  must  be  regarded  as  an  heuristic  model  at  best. 

Fortunately,  the  present  experimental  results  do  allow  a  simpli- 
fied development  which  can  provide  an  engineering  estimate  of  the 
effect  of  peak  surge  duration.  This  analysis  is  presented  in  the 
following  sections. 
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Experimental  Results 

As  previously  mentioned,  the  beach  profiles  were  recorded  for  each 
experiment  when  peak  surge  was  reached,  and  then  at  15  minute  inter- 
vals, maintaining  a  constant  peak  surge  level.  All  runs  were  continued 
at  peak  surge  until  a  near-equilibrium  appeared  to  be  established. 
This  typically  took  about  30  minutes,  which  represents  5.2  hours  in 
the  prototype.  Very  few  storms  will  have  a  peak  surge  duration  of 
this  length,  but  those  that  do  will  probably  reach  a  storm  equilibrium. 

The  profiles  obtained  for  each  experimental  run  are  given  in 
Appendix  C.  A  quick  inspection  of  these  profiles  gives  a  good  indi- 
cation of  the  large  amount  of  erosion  which  did  occur  by  the  time  the 
peak  surge  level  was  reached.  A  simple  computer  program  was  used  to 
calculate  the  volume  of  sediment  eroded  per  inch  of  beach  width. 
These  results  are  presented  in  Table  7. 

Preliminary  Data  Analysis 

Assuming  that  the  erosion  obtained  from  the  final  profile  for  each 
experiment  (involving  time-dependent  surges)  represents  the  equilibrium 
erosion,  the  erosion  values  given  at  peak  surge,  and  at  peak  surge  plus 
15  minutes,  can  be  expressed  as  a  percentage  of  the  equilibrium  erosion. 
For  the  33  runs  used  in  this  analysis,  it  was  found  that,  at  the  time 
of  peak  surge,  the  percentage  of  total  erosion  which  had  occurred 
ranged  between  60  percent  and  91  percent  with  the  overall  average 
being  79  percent.  Similarly,  the  percentages,  15  minutes  after  peak 
surge  had  been  reached,  ranged  between  80  percent  and  98  percent,  with 
an  overall  average  of  90  percent.  Table  8  presents  a  further  breakdown 
of  percentage  values  when  averaged  by  wave  period.  As  can  be  seen. 
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there  appears  to  be  no  distinguishable  trend,  and  the  values  vary  by 
only  7  percent  at  the  most. 

TABLE  8 
FRACTION  OF  EQUILIBRIUM  EROSION  FOR  WAVE  PERIODS 


1 

Period 
(sec) 

Average  at  Peak  Surge 

Average  at  Peak  Surge 
+  15  minutes 

0.970 
1.163 
1.357 

0.788 
0.774 
0.840 

0.904 
0.895 
0.930 

When  the  erosion  percentages  were  averaged  according  to  surge 
level,  however,  a  definite  trend  was  observed,  as  can  be  seen  in 
Table  9. 

TABLE  9 
FRACTION  OF  EQUILIBRIUM  EROSION  FOR  SURGE  LEVELS 


Surge  Level 
(inch) 

Average  at  Peak  Surge 

Average  at  Peak  Surge 
+  15  minutes 

3.84 
5.28 
6.72 

0.714 
0.798 
0.838 

0.863 
0.908 
0.930 

At  peak  surge,  there  was  a  12  percent  variation,  with  a  7  percent 
variation  given  for  the  values  obtained  15  minutes  later. 
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This  trend  indicates  that  a  higher  surge  level  promotes  a  larger 
percentage  of  erosion  during  the  surge  level  rise. 

Fortunately,  it  may  be  possible  to  neglect  this  dependency,  for 
the  present  analysis,  because  the  variations  given  in  Table  9  are 
within  7.5  percent  of  the  average  at  worst,  and  this  amount  of  error 
can  certainly  be  tolerated  in  a  time-dependent  erosion  solution. 

Engineering  Guidelines 

Using  the  overall  averages  calculated  for  the  percentage  of 
equilibrium  erosion  obtained  at  peak  surge  and  15  minutes  after  peak 
surge  had  been  reached,  it  was  possible  to  obtain  a  nearly  linear  rela- 
tionship as  shown  in  Figure  38. 

This  linear  approximation  is  given  as 


^i^  =   0.794  +  0.007  t  (77) 


for  the  model,  where  t  has  units  of  minutes.  The  corresponding  proto- 
type equation  becomes 

^^  =   0.794  +  0.041  t  ,  (78) 

e  ^ 

where  t  is  the  prototype  peak  surge  duration  expressed  in  hours.  In 

both  of  the  above  equations,  V(t)  is  the  transient  erosion  volume  per 

unit  width  of  beach,  while  V  is  the  equilibrium  erosion  volume  per 

unit  width  obtained  using  the  method  presented  in  Chapter  VIII. 
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10                               20 
MODEL  TIME  (MINUTES) 
J I \ \ I  I 


2  3  4 
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30 


PEAK     SURGE    DURATION 


FIGURE  38:  TIME-DEPENDENT  EROSION  FROM  PEAK  SURGE  TO  EQUILIBRIUM. 
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Discussion 


The  experimental  results  have  yielded  a  simple  means  of  adjusting 
the  dune  erosion  with  respect  to  peak  surge  duration.  While  there 
exists  some  variation  in  the  percentages  of  erosion  given  at  peak  surge, 
these  variations  dre   within  10  percent.  When  considering  the  range  of 
parameters  tested  in  the  model  series,  these  results  are  quite  satis- 
fying. The  important  point  to  notice  is  that  a  large  percentage  of 
erosion  had  occurred  by  the  time  the  storm  reached  peak  surge.  Thus, 
the  difference  in  erosion  between  a  fast-moving  and  a  slow-moving  storm 
is  only  about  15  percent  of  the  total  equilibrium  erosion. 

From  this  observation,  it  can  be  concluded  that  peak  surge  dura- 
tion plays  an  important  role  only  when  the  dune  heights  are  at  a  rela- 
tively low  elevation.  In  this  case,  a  greater  dune  recession  must  take 
place  in  order  to  erode  the  volume  of  sediment  necessary  to  achieve  an 
equilibrium.  This  is  readily  seen  in  the  Appendix  C  profiles  when  the 
surge  level  is  high  or  the  dune  elevation  is  low. 

This  whole  analysis  is  based  on  the  assumption  that  the  morpho- 
logical time  scaling  is  the  same  as  the  hydraulic  time  scaling.  Argu- 
ments supporting  this  assumption  appear  in  Chapter  IV,  and  the  verifi- 
cation results  of  Chapter  V  seem  to  support  this  claim.  This  aspect, 
however,  urgently  needs  to  be  further  verified,  and  as  such,  becomes  a 
top  priority  for  further  research. 

For  the  present,  it  is  felt  that  the  linear  partitioning  of  surge 
duration-dependent  erosion,  given  by  equation  (78),  represents  a  work- 
able method  of  estimating  dune  erosion  during  severe  storms. 
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Appli cation 

For  a  given  storm  for  which  the  peak  surge  duration  is  known  (along 
with  the  storm  parameters  needed  to  apply  the  results  of  Chapter  VIII), 
the  erosion  estimation  is  carried  out  in  the  following  manner: 

1)  Calculate  the  equilibrium  barred  storm  profile  using  the  pro- 
cedure outlined  in  Chapter  VIII. 

2)  Transpose  this  profile,  relative  to  peak  surge,  along  the 
prestorm  profile  until  a  balance  is  obtained  between  erosion 
and  deposition.  This  configuration  represents  the  maximum 
dune  erosion  expected  to  take  place. 

3)  Calculate  the  equilibrium  erosion  volume,  V  . 

'  e 

4)  For  the  given  peak  surge  duration  (less  than  5  hours),  calcu- 
late the  transient  erosion  volume,  V(t),  using  equation  (78). 

5)  Find  the  difference  between  V  and  V(t),  and  reduce  the  dune 
recession  distance  until  this  difference  in  erosion  volume  has 
been  "replaced."  This  will  depend  upon  the  prestorm  dune 
elevation. 

A  first  estimate  can  be  found  by  assuming  the  transient  beachface 
profile  is  the  same  as  the  near-equilibrium  profile.  While  not  correct, 
the  Appendix  C  profiles  indicate  that  this  assumption  is  fairly  reasonable. 

The  method  only  gives  an  estimation  of  the  transient  dune  re- 
cession. The  rest  of  the  calculated  profile  represents  equilibrium 
conditions. 

For  high  elevation  dunes,  the  method  should  provide  fairly  accu- 
rate results.  As  the  dune  elevation  decreases,  this  accuracy  will 
fall  off  somewhat. 


CHAPTER  X 
SUMMARY  AND  CONCLUSIONS 

The  ability, to  accurately  predict  the  amount  of  dune  recession 
during  extreme  storm  events,  such  as  hurricanes,  is  one  of  the  most 
important  problem  solutions  required  by  coastal  planners.  This  ability 
would  assist  coastal  engineers  in  providing  safe  recommendations  for 
the  increasing  coastal  construction  being  experienced  on  many  portions 
of  the  coastline.  This  research  has  been  directed  toward  the  establish- 
ment of  a  reliable  engineering  method  for  the  prediction  of  beach-dune 
erosion  under  any  specified  storm  condition. 

Model  Law 

Movable-bed  small-scale  modeling  of  coastal  erosion  has  been 
historically  a  very  inexact  science.  The  modeling  relationships  de- 
veloped in  Chapter  IV  were  derived  through  the  consideration  of  the 
inertial  forces,  represented  by  the  turbulent  shear  stress,  and  the 
gravity  force  in  the  nearly  horizontal  direction  of  the  principal  flow. 
This  resulted  in  a  dynamic  scaling  relationship  for  a  distorted  model. 
The  profile  similarity  criterion  was  based  on  the  preservation  of  the 
dimensionless  fall  velocity  parameter,  which  has  been  shown  to  be  an 
important  factor.  Combining  of  these  two  criteria  resulted  in  a  model 
law  which  is  similar  to  previous  results,  differing  primarily  in  the 
hydrological  time  scaling.  By  proper  selection  of  model  grain  size. 
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it  is  possible  to  have  an  undistorted  model,  but  this  will  rarely  be 
possible. 

It  is  concluded  that  the  modeling  relationships  presented  repre- 
sent a  plausible  means  of  scaling  beach-dune  erosion  in  the  small-scale 
physical  model.  The  results  of  the  model  verification  further  add  to 
the  strength  of  this  conclusion. 

Model  Verification 

The  verification  of  the  derived  modeling  relationships  was  per- 
formed in  the  wave  flume  by  the  close  approximation  of  an  actual  proto- 
type event.  By  careful  examination  of  the  limited  data  available,  it 
was  possible  to  reasonably  estimate  the  prototype  storm  parameters. 
Reproduction  of  these  parameters  in  the  wave  tank  gave  an  almost  exact 
reproduction  of  the  actual  erosion  experienced  in  the  prototype.  Further 
credibility  can  be  given  to  the  verification  by  the  fact  that  the  pro- 
cedure involved  a  time-dependent  storm  surge  hydrograph  and  an  increas- 
ing wave  height  as  the  storm  progressed.  To  the  author's  knowledge, 
this  was  the  first  time  that  these  features  had  been  incorporated  into 
an  erosion  model  verification. 

Attempts  to  verify  the  model  using  a  randomly  produced  irregular 
wave  train,  fitting  the  spectral  requirements,  proved  unsuccessful. 
Arguments  against  the  use  of  irregular  waves  have  been  given,  but  later 
developments  may  make  their  use  feasible. 

It  is  concluded  that  the  model  verification  has  been  successful 
in,  at  least,  one  instance  for  which  prototype  data  were  available. 
The  reproduction  of  the  storm  event,  to  a  degree  of  realism  never  be- 
fore attempted,  certainly  adds  weight  to  this  conclusion. 


■184- 


Experimental  Test  Series 

The  experimental  test  series  was  designed  around  a  typical  Florida 
beach  profile  with  sand  grain-size  in  the  vicinity  of  0.26  mm.  The 
series  looked  at  the  effects  in  varying  the  wave  period,  wave  height, 
peak  surge  elevation,  dune  height,  and  surge  level  rise  duration.  The 
parameters  were  chosen  to  represent  the  upper  and  lower  values  for  a 
severe  storm.  The  variable-surge  experiments  began  with  an  equilibrium 
profile  and  then  elevated  the  surge  level  over  a  typical  time  span  ex- 
pected to  occur.  Upon  reaching  the  peak  surge  level,  the  test  was 
continued  until  near-equilibrium  conditions  prevailed. 

These  time-dependent  surge  rise  experiments  represent  the  first  of 
this  type,  and  the  results  indicated  that  the  finite  surge  rise  is 
necessary  to  more  accurately  depict  the  dune  erosion  phenomenon. 

The  testing  procedure  was  carefully  monitored  and  controlled,  and 
the  conclusion  is  made  that  the  results  represent  the  most  valid  dune 
erosion  experiments  to  date. 

Qualitative  Profile  Analysis 

Qualitative  comparisons  between  the  resulting  near-equilibrium 
profiles  resulted  in  several  conclusions:  It  is  necessary  to  employ  a 
time-dependent  surge  level  increase,  there  appears  to  exist  an  equilib- 
rium barred  profile  which  is  characteristic  of  the  storm  parameters, 
erosion  can  be  found  by  shifting  the  barred  profile  to  obtain  a  sediment 
balance,  the  storm  surge  level  is  the  most  important  factor  in  dune 
erosion,  and  the  profile  between  the  dune  toe  and  the  offshore  bar  is 
independent  of  surge  level  and  prestorm  profile  configuration.  In 
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addition,  some  aspects  of  the  dune  erosion  process  were  discussed  in 
general  terms. 

Equilibrium  Barred  Storm  Profile 

In  Chapter  VIII  it  was  shown  that  several  features  of  the  equi- 
librium barred  storm  profile  could  be  expressed  in  terms  of  the  dimen- 
sionless  fall  velocity  parameter,  thus  adding  further  strength  to  the 
modeling  profile  similarity  criterion. 

Relationships  were  established  for  the  bar  crest  and  bar  trough 
distances  from  the  dune  toe  in  a  nondimensional  form,  which  is  pre- 
served by  the  model  law.  It  was  seen  that  these  relationships  could  be 
physically  explained  in  terms  of  incoming  wave  energy  flux,  per  unit 
width,  entering  the  surf  zone. 

The  nearshore  profile  was  found  to  fit  a  power  curve  representa- 
tion with  a  coefficient  based  solely  on  grain-size  fall  velocity.  It 
was  seen  earlier  that  this  form  also  was  characteristic  of  the  model 
1  aw. 

A  simple  nondimensional  parameter  was  developed  and  correlated 
for  the  vertical  runup  distance.  Having  nearly  the  same  form  as  the 
dimensionless  fall  velocity  parameter,  the  runup  parameter  is  also 
preserved  by  the  model  law. 

An  analysis  into  the  offshore  portion  of  the  barred  profile  yielded 
results  applicable  to  the  average  Florida  beach,  but  further  research 
needs  to  be  conducted  before  applying  the  method  elsewhere. 

Finally,  a  method  was  given  for  use  in  calculating  the  equilibrium 
barred  profile,  using  the  given  storm  parameters.  With  this  calculated 
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profile,  it  becomes  possible  to  determine  the  maximum  equilibrium 
erosion  expected  to  occur  for  any  specified  storm  on  any  given  beach 
profile. 

It  is  concluded  that  the  method  provides  results  superior  to  the 
previous,  simplified  methods.  This  was  the  first  analysis  to  include 
the  offshore  bar  feature. 

Transient  Dune  Erosion 

An  examination  of  the  experimental  profiles  revealed  that  between 
70  and  90  percent  of  the  final  equilibrium  erosion  had  taken  place  by 
the  time  the  peak  surge  level  had  been  reached, 

A  theoretical  development,  based  on  energy  concepts,  was  pre- 
sented, but  it  was  found  that  it  could  not  be  applied  in  this  instance. 
However,  a  simple  engineering  approximation  was  given  for  the  linear 
partitioning  of  the  erosion  which  occurs  during  peak  surge  duration. 
Using  this  method,  it  becomes  possible  to  predict  the  dune  recession 
as  a  function  of  peak  surge  duration,  making  it  possible  to  allow  for 
situations  where  equilibrium  could  not  be  reached  over  the  time  span 
of  the  storm  event. 

Application  and  Assumptions 

The  proposed  method  for  the  prediction  of  dune  erosion  is  limited 
in  its  application  by  the  following  assumptions: 

1)  Natural  beaches  with  sand  grain-size  distributions  having  a 
mean  diameter  of  around  0.25  mm.  While  the  relationships 
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contain  provision  for  different  grain-sizes,  they  have  not 
been  fully  proven,  as  yet. 

2)  Straight  beaches  with  fairly  uniform  offshore  depth  contours. 

3)  No  interactions  with  tidal  inlets  or  coastal  structures. 

4)  Onshore-offshore  sediment  transport  with  no  alongshore 
currents. 

5)  No  dune  overtopping. 

In  addition,  the  offshore  prestorm  profile  must  be  similar  to  that 
modeled  in  the  experimental  test  series.  While  small  variations  in 
this  profile  would  probably  have  little  effect,  major  differences 
would  certainly  effect  the  predictive  ability  of  the  method.  For- 
tunately, many  undeveloped  Florida  beaches  fall  into  the  above  stipu- 
lations; thus,  the  method  may  be  used  to  aid  in  the  analyses  of  the 
Coastal  Construction  Setback  Line. 

Summary 

It  is  felt  that  this  method  represents  an  advance  in  the  parameter- 
ization of  the  dune  erosion  process,  and  it  provides  an  improved  means 
of  determining  dune  erosion  resulting  from  severe  storms.  In  addition, 
several  surf  zone  characteristics  can  now  be  described  by  physically 
plausible  relationships. 

The  research  has  also  provided  a  new  small -scaling  movable-bed 
modeling  law  which  has  been  verified  under  realistic  conditions.  It 
is  hoped  that  further  research  will  add  greater  confidence  to  the 
modeling  relationships,  and  thus,  pave  the  way  for  routine  physical 
modeling  of  coastal  processes. 
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Recommendations  for  the  Future 


It  is  most  imperative  that  work  continue  into  the  verification  of 
the  modeling  relationships.  This  research  should  confirm  the  existence 
of  profile  similarity  between  models  of  different  scale  and  different 
grain-sizes.  It  would  also  be  beneficial  to  determine  how  small  the 
model  can  be  made. 

Further  verification  of  prototype  events  should  be  attempted  when- 
ever accurate  prototype  data  become  available.  This  would  further 
enhance  the  reliability  of  the  model  law  and  give  a  more  complete 
understanding  of  the  scaling  effects.  When  it  becomes  practical  to 
use  wave  spectra  in  a  movable-bed  model,  then  they  must  be  employed. 

More  testing  should  be  done  using  different  grain-sizes  and  dif- 
ferent offshore  profiles  to  confirm  the  empirically  derived  expressions 
used  in  the  calculation  of  the  equilibrium-barred  storm  profile,  and 
to  generalize  the  method  for  the  application  to  any  sandy  shoreline. 

After  the  above  research  has  been  conducted,  and  the  method  has 
been  improved  upon,  it  would  be  possible  to  examine  the  effects  of  waves 
overtopping  the  dune,  man-made  structures,  tidal  inlets,  and  three- 
dimensional  processes  involving  alongshore  currents. 


APPENDIX  A 
COMPUTER  PROGRAM  LISTINGS 


PROFILE  DATA  TAKING  PROGRAM 
C 

C  THIS  PROGRAM  WILL  RECORD  THE  DUNE  PROFILE  WHEN  THE  PROBE  IS 

C  MANUALLY  PLACED  AT  THE  RIGHT  ELEVATION  AND  THE  BUTTON  IS  PUSHED. 

C  THE  CART  MUST  THEN  BE  MOVED  TO  THE  NEXT  POINT,  AND  THE  PROCESS 

C  REPEATED. 
C 

C  IT  WILL  ALSO  RECORD  THE  UNDERWATER  PROFILE  AUTOMATICALLY. 
C 

C  THE  CART  SHOULD  ALWAYS  MOVE  FROM  THE  MOST  LANDWARD  PORTION  OF 

C  THE  DUNE  TOWARD  THE  WAVE  PADDLE  OFFSHORE. 

C   STEVE  HUGHES   3-JUNE-1980 
C 

LOGICAL*!  IVAL,Y 

DIMENSION  IPOINT(2048) 

INTEGER*2  CMONTH , CDATE , CYEAR, lEXP, IPR0F,HDRWD1 

BYTE  DIA(9) 

BYTE  FILENM(14) 

DATA  FILENM/ 'D','X','1',':','E','X','X','P','X','X', 

1' . ' , 'P' , ' R' , '0'/ 

DATA  Y/"131/ 

C   CHECK  FOR  DATE 
C 

CALL  I DATE (CMONTH, CDATE, CYEAR) 

IF ( CMONTH. EQ.O)  STOP  'ENTER  SYSTEM  DATE' 

C   IDENTIFY  PROGRAM  TO  USER 

TYPE  1 
1  FORMAT ( '      *******************************'   / 

1'      *  PROFILE  DATA  TAKING  PROGRAM  *',/, 
2'     *******************************'  j  is 


C   OUTPUT  CURRENT  SYSTEM  DATE 

CALL  DATE(DIA) 

TYPE  400, DIA 
400  FORMAT (3X, 9 Al) 

C 

C   GATHER  INFORMATION  FROM  OPERATOR 
C 

TYPE  2 

2  FORMAT ( 3 X, ' ENTER  CURRENT  EXPERIMENT  NUMBER:',?) 
ACCEPT  3 , lEXP 

3  FORMAT (12) 
TYPE  4 

4  FORMAT! 3X, ' ENTER  THE  CONSECUTIVE  PROFILE  NUMBER',/ 
13X,'F0R  THIS  EXPERIMENT:',?) 

ACCEPT  3, IPROF 
IEXP1=IEXP/10 
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IEXP2=IEXP- ( 10*IEXP1 ) 
IPROF1=IPROF/10 
IPR0F2=IPR0F- ( 10*IPROF1 ) 
FILENM( 6 )=IEXPl+4a 
F ILENM ( 7 ) = IEXP2+48 
FILENM(9)=IPROFl+48 
FILSlIM(iO)  =  IPROF2+48 
TYPE  5 

5  FORMAT ( 3 X, 'ENTER  REAL  VALUE  ( F6 . 2 )  OF  ELAPSED',/ 
13X,'TIME  IN  MINUTES  SINCE  START  OF  EXPERIMENT:',*) 
ACCEPT  6 , ELAPTM 

6  FORMAT (F6. 2) 

C   OUTPUT  FILE  NAME 

TYPE  7 , F ILENM 

7  FORMATOX,  'FILE  NAME  IS:  ',I4Ai) 
TYPE  7  00, lEXP, IPROF, ELAPTM 

700         FORMAT (3X, 'EXPERIMENT  #:', 12 , 2X, ' PROFILE  ?:', 12 , 2X, ' ELAPSED 

i,F6.2,iX, 'MINUTES' ,/) 
C 

C   OPEN  UNFORMATTED  OUTPUT  FILE 
C 

CALL  ASSIGN  ( 1 , FILENM, 14 , ' NEW' ) 

DEFINE  FILE  1  (1 , 2048 , U , INEXT ) 

HDRWD1=-1 

C   DEFINE  VERITICAL  CALIBRATION  CONSTANT 

CVERT=a.25 
C   DETERMINE  MEAN  SEA  LEVEL 

TYPE  10 
10  FORMAT (3X, 'PLACE  PROBE  AT  MEAN  SEA  LEVEL  POSITION') 

PAUSE  'TYPE  <CR>  WHEN  READY. ' 

JOE=0 

DO  45  1=1,5 
45  J0E=J0E+IADC(14) 

IFIRST=INT(JOE/5.0) 

C   DUNE  OR  OFFSHORE  PROFILE? 

L=l 
TYPE  90 

90  FORMAT ( 3X ,' ARE  YOU  TAKING  DUNE  PROFILE?  (Y  OR  N):',5) 
ACCEPT  91,IVAL 

91  FORMAT (Al) 
IFtlVAL.NE.Y)GO  TO  305 

C   DUNE  PROFILE  MEASUREMENTS 
C 

TYPE  9  2 
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92  FORMATOX, 'MOVE  CART  TO  INITIAL  DUNE  POSITION  AND  INPUT  INITIAL',/ 
1,3X, 'HORIZONTAL  CART  POSITION  (NEGATIVE  INCHES  F6 . 2 ) : ' , ? ) 
ACCEPT  6,  HORIZ 

TYPE  93 

93  format ( 3x, 'push  button  to  manually  record  dune  profile',/, 
13x, 'type  <cr>  on  completion') 

70  ;hi=o 

call  take ( ihi ) 

IF(IHI)72,73,74 

72  JOE=0 

DO  46  11  =  1.,  5 
46  JOE=JOE+IADC( 14) 

IPOINT(L)=INT( JOE/5.0) 
GO  TO  70 

73  L=L+1 
IPOINT(L)»0 
GO  TO  70 

74  CONTINUE 
C 

C   OFFSHORE  PROFILE? 

TYPE  31 
ai  FORMAT ( 3X ,' ARE  YOU  TAKING  OFFSHORE  PROFILE?  (Y  OR  N):',5) 
ACCEPT  9 1 , IVAL 
IF(IVAL.NE.Y)J=L 
IF(IVAL.NE.Y)GO  TO  40 
TYPE  32 

32  FORMAT t3X, 'MOVE  CART  TO  INITIAL  OFFSHORE  POSITION,  TYPE  <CR>  WHEN 
1  IN  POSITION' ) 

ICHAR=0 
36  IHI=0 

CALL  TAKE ( IHI ) 
IF(IHI)a3,34,aS 

33  GO  TO  36 

84  IPOINT(L)=0 

L=L+1 

GO  TO  36 
35  GO  TO  802 

C   BRANCH  FOR  NO  DUNE  PROFILE 

305         TYPE  301 

801  FORMAT ( 3X ,' MOVE  CART  TO  INITIAL  POSITION  AND  INPUT  HORIZONTAL' 

13X,  'POSITION  IN  INCHES  ( F6  .  2 )  :  '  , ? ) 

ACCEPT  6, HORIZ 

C   RECORD  UNDERWATER  PROFILE 
C 
302         J=L 

ICHAR=ITTINR( ) 
TYPE  12 
12  F0RMAT(3X, 'TURN  ON  PROFILER  AND  TYPE  <CR>  TO  BEGIN  PROFILE  DATA',/ 
13X, 'TAKING.   'rfHEN  COMPLETED,  STOP  CART  AND  TYPE  <CR> . ' ) 
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PAUSE  'TYPE  <CR>  TO  START' 

C 

IHI-0 
15  CAU,  TICK(IHI) 

IP0INT(J)=IADC(14) 
J=J+i 

IF(IHI.GT.0)GO  TO  40 
GO  TO  15 

C   STORE  HEADER  WORDS  AND  PROFILE  DATA  ON  DISK 

40  K=J-1 

M-2048-11-K 
C 
C   CALCULATE  HORIZONTAL  CALIBRATION  CONSTANT 

TYPE  560 
560         FORMAT ( 3X, ' INPUT  FINAL  HORIZONTAL  POSITION  { F6 . 2 )  IN  INCHES-'  S) 

ACCEPT  6, FINAL 

CH0RI2»(FINAL-H0RIZ)/K 

WRITE (I'D  HDRWDi , ELAPTM , CHORIZ , CVERT , HORIZ , IFIRST , HDRWDl , 

KIPOINTd)  ,I=1,K)  ,  (HDRWDl,  1  =  1,  M) 

CALL  CLOSE ( 1 ) 
550         TYPE  300 

300  FORMAT (3X, 'OPTIONAL  DATA  DISPLAY!',//, 
13X, "TYPE  [Ij  FOR  NUMERICAL  RESULTS',//, 

23X, 'TYPE  [2]  FOR  GRAPHICAL  DISPLAY',//, 
3  3X,'TYPE  <CR>  TO  EXIT  PROGRAM') 
ACCEPT  301, I OP 

301  FORMAT!  ID 
IFdOP.NE.DGO    TO    302 

CALL    DATA (K, CHORIZ, CVERT, HORIZ, IFIRST, IPOINT) 
30    TO    550 

302  IF( I0P.NE.2)G0    TO    350 

CALL  PLOT (K, CHORIZ, CVERT, HORIZ, IFIRST , IPOINT ) 
GO  TO  550 
3  50         STOP 
END 

SUBROUTINE  TICK(IHI) 
10  IHI=IADC(15) 

ICHAR=ITTINR( ) 

IFdCHAR.GT.OlaO    TO    30 

IF(IHI.LT.2460)GO    TO    10 
20    IHI=IADC(15) 

ICHAR=ITTINR( ) 

IF(ICHAR.GT.0)GO  TO  30 

IF(IHI.GT.2460)GO  TO  20 

IHI»0 

RETURN 
30  IHI=1 

RETURN 


■194- 


SUBROUTINE  TAKE(IHI) 

IYES=0 
10  IHI=IADC{15) 

IYES=IADC(13) 

ICHAR=ITTINR( ) 

IF( ICHAR.GT.0)GO  TO  30 

IF( IYES.aT.2458)GO  TO  40 

IF( IHI .LT.2460)GO  TO  10 
20  IHI=IADC(15) 

IYES=IADC(13) 

ICHAR=ITTINR( ) 

IF(ICHAR.GT.0)GO  TO  30 

IF(IYES.GT.2458)GO  TO  40 

IFdHI  .GT.2460)GO  TO  20 

IHi=0 

RETURN 
30  IHI=i 

RETURN 

40  IYES=IADC(13) 
IF(IYES.LT.2458)GO  TO  41 
GO  TO  40 

41  IHI=-1 
RETURN 
END 


SUBROUTINE  DATA( J , CHORIZ , OVERT , HORIZ , IFIRST, IPOINT) 
DIMENSION  IPOINT (2000) 


REQUEST  DO-LOOP  INCREMENT 


TYPE  5 
5   FORMAT ( 3X ,' ENTER  VALUE  OF  [N]  TO  PRINT  EVERY  N-TH  POINT  (12)') 

ACCEPT  6 , N 
5   FORMAT (12) 

TYPE  7 
7   F0RMAT(3X, ' DATA  POINT ', 3X ,' RANGE ( INCHES )', 3X ,' ELEV . WRT . MSL  (INCHES 

1)  '  ) 

DO  10  I=1,J,N 

RANGE=HORIZ+(CHORIZ*(I-l) ) 

IF(  IPOINTd  )  .EQ.0)GO  TO  10 

VERT= ( IFIRST-IPOINT ( I ) ) *CVERT/409 . 6 

TYPE  2  0, I, RANGE, VERT 
20  FORMATOX,  I3,9X,F7.2,  11X,F6.2) 

GO  TO  10 
10  CONTINUE 

RETURN 

END 
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SUBROUTINE  PLOT ( J , CHORIZ , CVERT , HORIZ , I FIRST, IPOINT) 
L0GICAL*1  STARS (65) 
DIMENSION  IPOINT (2000) 
DATA  STARS/65* ' *' / 

:       REQUEST    DO-LOOP     INCREMENT 

TYPE  1 

1  FORMAT ( 3X, ' ENTER  VALUE  OF  [N]  TO  PRINT  EVERY  N-TH  POINT  (12)') 
ACCEPT  2,N 

2  FORMAT (12) 

TYPE  OUT  LABELS 


TYPE  3 

3  FORMAT (2X, 'RANGE  ELEV .  MSL ' , 4X, ' ABS .  VALUE  OF  ELEV.(*=l/4  INCH)' 
1/,1X, ' (INCH)  (INCH)  .',6X,'(N)  DENOTES  NO  DATA  TAKEN  THIS  RANGE' 
TYPE  4 

4  F0RMAT(1X,7  5(1H-) ) 
C 

C   PLOT  DATA 
C 

DO  10  1=1, J, N 

RANGE=HORIZ+ (CHORIZ* (I-l) ) 

IFdPOINTd)  .EQ.O  )G0  TO  30 

VERT=( IFIRST- IPOINT ( I ) ) *CVERT/409 . 6 

IFT=INT( (ABS (VERT) *4.0) +0.5) 

TYPE  20, RANGE, VERT, ( STARS ( M ) ,M=1, IFT ) 
20  F0RMAT(IX,F6.2,1X,F6.2, lX,65Ai) 

GO  TO  10 
30  TYPE  3  5, RANGE 

3  5  FORMAT ( IX , F6 . 2 , 3X ,  '  —  '  , 3X ,  ' N '  ) 
10  CONTINUE 

RETURN 

END 
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C   THIS  PROGRAM  OPENS  DATA  FILES  COLLECTED  USING 

C   THE  PROFILE  PROGRAM  AND  ANALYSES  THEM. 

C 

C   STEVE  HUGHES    NOV.  24,  1980 

C 

INTEGER* 2  HDRWDl 

DIMENSION  IPOINT(2040) , RANGE (100) , VERT (100) ,E(100) 

DIMENSION  A(IOO) ,ATOT(100) 

TYPE  1 
1   FORMAT (3X, 'ENTER  DATA  FILE  SAME:',?) 

CALL  ASSIGN  ( 1 , ' XXXX . PRO ' , -8 , ' RDO ' , ' NC ' , 1 ) 

DEFINE  FILE  1  ( 1 , 2048 , U, INEXT ) 

INEXT=i 
C 

C   READ  IN  DATA  FILE 
C 

READ ( 1 ' INEXT ) HDRWDl , ELAPTM, CHORIZ , CVERT , HORIZ , IFIRST , HDRWDl 

KIPOINTd)  ,  1  =  1,2037) 
C 

C       CONVERT    TO    REAL    VALUES 
C 

Ria=409.6 

K=l 

DO    3    1=1,2037 

IFdPOINTd)  .EQ.-1)G0  TO  5 

IFdPOINTd)  .EQ.0)GO  TO  3 

RANGE(K)=HORIZ+( d-l)*CHORIZ) 

VERT(K)=(IFIRST-IPOINT(I) ) * ( CVERT/RIB ) 

K=K  +  1 
3   CONTINUE 

5  CONTINUE 
L=K-1 

C   DETERMINE  PRE-STORM  DUNE  HEIGHT 
C 

IF(VERTd)  .LT.8.0)GO  TO  28 

IF(VERT(1) .GT.20.0)GO  TO  28 

DUNHT=VERT(1) 

6  FORMAT ( F6 . 2 ) 
GO  TO  29 

23  DUNHT=VERT(2) 

VERT(1)=VERT{2) 
29  CONTINUE 

C   CALCULATE  PRE-STORM  PROFILE  ARRAY 

E3=0.0 

R3=0.0 

E2=3.36 

R2=-16.289 

E1=DUNHT 

R1=-(E1-E2)*0.969+R2 

DO  7  1=1, L 
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IF(RANGE(I)  .GT.RDGO    TO    10 

E(I)=E1 

GO    TO    7 

10  XF(RANGE(I ) .GT.R2)G0  TO  11 

E ( I)=- ( RANGE ( I ) -R2 ) /O . 969+E2 
GO  TO  7 

11  IF(RANGE(I) .GT.R3 )G0  TO  12 
E(I)=-{RANGE(I)-R3)/4.347+E3 
GO  TO  7 

12  E(I)=-O.19O3*(RANGE(I)**(2.0/3.O) ) 
7   CONTINUE 

C 

C   INPUT  VALUE  OF  HORIZONTAL  REFERENCE 

C 

TYPE  18 
18  FORMAT! 3X, 'ENTER  ELEVATION  FOR  NET  EROSION  (INCH . F6 . 2 ) : ' , 5 ) 

ACCEPT  6,ALEVEL 
C 

C   CALCULATE  EROSION  AREA 
C 

K=L-1 

QO  13  1=1, K 

C=RANGE ( I+l ) -RANGE ( I ) 

IF(RANGE(I) .LT.R1.AND.RANGE( I+l) .GT.RDGO    TO    14 

IF(RANGE(I) .LT.R2.AND.RANGE(I+1) .GT.R2)G0    TO    15 

IF(RANGE(I) .GT.R3)G0    TO    16 

ATOT(I)=(ECl)-VERT( I ) +E ( I+l ) -VERT ( l+l ) )*C/2.0 

IF(VERT(I+1)  .LT.ALEVEDGO    TO    17 

A{I)=ATOT(I) 

GO    TO    13 
17    A(I)=(E(I)+E(I+1)-2.0*ALEVEL)*C/2.0 

GO    TO    13 

14  CALL  DISCON(VERT, RANGE,E, Rl.El, ALEVEL,TOT,RTOT, I) 
ATOT { I ) =TOT 

A( I)=RTOT 
GO  TO  13 

15  CALL  DI SCON (VERT, RANGE, E, R2,E2,ALEVEL, TOT, RTOT, I ) 
ATOT ( I ) =TOT 

A(I)=RTOT 
GO  TO  13 

16  T1=(2.0*E(I)-VERT(I)-VERT(I+1) )*C/2.0 
T2=0.1142*(RANGE(I)**(5.0/3 . 0 ) -RANGE ( I+l ) ** ( 5 . 0/3 . 0 ) ) 
T3=E(I)*C 

AT0T(I)=T1+T2-T3 
A(I)=0.0 

13  CONTINUE 
ATOT(L)=0.0 
A(L)=0.0 
TOTPOS=0.0 
TOTNEG=0 . 0 
APOS=0.0 

DO  20  1=1, L 
IF(ATOT(I)  ) 2 1,2 1,2 2 
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2 1  TOTNEG=TOTNEG+ATOT ( I ) 
GO  TO  23 

22  TOTPOS=TOTPOS+ATOT( I) 

23  IF(A(I) .GT.0.0)APOS=APOS+A(I) 
20  CONTIlWE 

::   PRINT  RESULTS 

TYPE  30 

30  FORMAT ( 3X ,' TYPE  [1]  TO  SUPPRESS  INTERMEDIATE  VALUES:',?) 
ACCEPT  31, I WHAT 

31  FORMAT (II) 

IF( IWHAT.EQ. 1 )G0  TO  35 
TYPE  3  2,ALEVEL 

32  F0RMAT(14X, ' FINAL' ,4X, ' INITIAL   DISPLACED   EROSION  ABOVE ' /4X, 

1' RANGE' ,5X, 'ELEV. ' ,5X, ' ELEV . ' , 5X , 'AREA' ,5X,F5.2,1X, ' INCH.  ELEV . 

23X,'(INCH.)    (INCH.)    (INCH.)   ( SQ .  IN.)  (SQ.  IN.)'/) 

DO  8  1=1, L 

TYPE  9 , RANGE ( I  )  , VERT ( I )  , E ( I )  , ATOT ( I  )  ,  A  (  I  ) 
9   FORMAT (3X,F7. 2, 3X,F7 . 2 , 3X, F7 . 2/ 32X, F8 . 2 , 3X, F8 . 2 ) 
3   CONTINUE 

35  TYPE  36,TOTPOS,TOTNEG, ALEVEL,APOS 

36  FORMAT (3X,  'TOTAL  EROSION=  '  , F8  .  2 ,  '  SQ  .  IN.'/3X, 
1' TOTAL  DEPOSITION=  ' , F8 . 2 , '  SQ .  IN.'/3X, 

2 'TOTAL  EROSION  ABOVE  ELEV.  '  , F5 . 2 ,  '  =  '  , F8  .  2 ,  '  SQ .  IN.') 

CALL  CLOSE(l) 

STOP 

END 


SUBROUTINE  DISCON ( V, R, E, RR, EE , ALEVEL, TOT, RTOT, I) 

DIMENSION  V(IOO) ,R(100) ,E(100) 

V1=(V(I)-V(I+1) )*{R(I+1)-RR)/(R(I+1)-R(I) )+V(I+l) 

Tl= ( E ( I ) +EE-V ( I ) -VI ) * ( RR-R( I ) ) / 2 . 0 

T2=(EE+E(H-l)-Vl-V(I+i) )*(R(I+l)-RR)/2.0 

T0T=T1+T2 

IF(V(I  +  1)  .LT.ALEVEDGO    TO    2 

RTOT=TOT 

GO    TO    3 

Tl= (E(I)+EE-2.0* ALEVEL) *(RR-R( I ) )/2.0 

T2=(EE+E(I+l)-2.0*ALEVEL)*(R(I+l)-RR)/2.0 

RT0T=T1+T2 

RETURl^ 

END 


APPENDIX  B 
INSTRUMENT  CALIBRATIONS 


VERTICAL  ELEVATION 
(CENTIMETERS) 


FIGURE   39:      WAVE  GAUGE   CALIBRATION, 
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APPENDIX  E.l 

DEAN'S  DERIVATION  FOR  A  BEACH  PROFILE  DUE  TO  UNIFORM 
ENERGY  DISSIPATION  PER  UNIT  VOLUME 


The  energy  flux  within  the  surf  zone  can  be  expressed  as 

If  " -a^  -hD(d)  ,  (79) 


where 


E  =  wave  energy  per  unit  width 

=  wave  group  V( 
h  =  water  depth, 


C  =  wave  group  velocity, 


x'  -  coordinate  directed  shoreward  from  breaking  point,  and 

D(d)  =  rate  of  wave  energy  dissipation  per  unit  volume  as  a 
function  of  grain  size. 

Equation  (79)  states  that  all  of  the  incoming  wave  energy  per  unit 

width  is  completely  dissipated  uniformly,  per  unit  volume,  across  the 

surf  zone.  In  other  words,  it  is  assumed  that  a  sediment  particle  can 

tolerate  a  certain  level  of  wave  energy  dissipation  per  unit  volume 

and  still  remain  stable.  In  actual  fact,  the  sediment  is  nearly  always 

in  motion,  so  stability  in  this  case  refers  to  net  transport  of 

sediment. 

By  considering  steady  conditions  when  averaging  over  a  wave 
period,  i.e.. 


^=   0 
dt   ^ 
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and  replacing  dx'  with  -dx,  which  is  a  coordinate  system  directed  sea- 
ward from  the  beach,  equation  (79)  becomes 


d(EC 


dx 


31  = 


hD(d) 


(80) 


From  linear  wave  theory. 


E-^    and  Cg  =  (gh)l/2 


which  can  be  used  in  equation  (80).  By  further  assuming  that  the  bore 
is  decreasing  in  height  proportional  to  the  depth,  then 

H  =  Kh  . 

Making  these  substitutions  into  equation  (80)  yields 


^2^1/2  .^.5/2v 
8     dx     ^^^°'    ' 


which  can  be  integrated  to  give 


^f^  h3/2  =  D(d)x 


Rearranging  gives 


h  =  Ax 


2/3 


(8i: 


where 


„  _  r24D(d) 
"  L.2  1/2 


i2/3 


-SyK'g' 


(82: 
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Here  it  is  seen  that  "A"  is  a  function  only  of  sediment  size, 
and  it  is  reasonable  to  believe  that  the  grain  fall  velocity  plays 
the  predominant  part  in  dissipating  the  turbulent  energy. 

As  was  pointed  out  in  the  chapter  dealing  with  the  model  law,  the 
form  given  by  equation  (81)  is  preserved  by  the  model  law  when 

which  means  "A"  is  proportional  to  u  '    .     This  further  aids  in  explaining 
why  profile  similarity  exists  when  the  parameter  H/uT  is  scaled  in 
the  model . 


APPENDIX  E.2 
DERIVATION  OF  THE  OFFSHORE  EQUILIBRIUM  PROFILE  EQUATION 

Quite  a  few  papers  have  been  written  on  the  subject  of  sediment 
transport  rates  and  equilibrium  beach  conditions  in  the  nearshore 
vicinity.  While  all  papers  have  their  merit,  perhaps  the  most  widely 
quoted  on  the  subject  of  surf  zone  energetics  are  Bagnold  (1963),  Inman 
and  Bagnold  (1963),  and  more  recently,  Madsen  (1976).  All  three 
papers  deal  with  the  mechanics  of  nearshore  processes  in  a  thorough 
physical  manner. 

Inman  and  Bagnold  developed  a  model  for  the  equilibrium  beach 
slope  given  as 


tane  =  tand)  (y^)  ,  (83) 


where 


tanB  -   equilibrium  beach  slope, 

tantj)  =  coefficient  of  intergranular  friction,  and 

_  energy  dissipated  during  the  seaward  sediment  flow 
energy  dissipated  during  the  shoreward  sediment  flow  * 

Essentially,  equation  (83)  is  a  balance  between  gravity  which  tends  to 
move  sediment  offshore,  and  the  non-symmetrical  bottom  velocity  dis- 
tribution which  moves  sediment  onshore.  To  illustrate  this,  consider 
a  symmetrical  velocity  distribution  in  which  c  =  1.  Equation  (83)  then 
gives  tang  =0,  or  a  horizontal  bed.  Thus,  in  order  to  have  a  beach 
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slope,  it  is  necessary  to  have  an  oscillatory  bottom  velocity  distri- 
bution which  is  skewed  in  the  onshore  direction. 

In  the  nearshore  region,  seaward  of  the  breaking  point,  beach 
slopes  greater  than  0.1  (1  in  10)  would  be  highly  unusual.  For  small 
slopes  such  as  these, 

tang  Ti  B  in  radians  =  -^   . 

Making  this  substitution  into  equation  (83)  and  solving  for  "c"  yields 

^   L  dh 

tani})  dx 

As  an  approximation  for  the  ratio  of  energy  dissipation,  c,  Inman 
and  Frautschy  (1965)  proposed  that 


U    3 
Voff 


u 

m-on 


(85) 


where 


U^_Qff  =  maximum  offshore  component  of  the  orbital  velocity  and 

U„  ^^  -  maximum  onshore  component  of  the  orbital  velocity, 
m-on  •' 

More  recently,  Bailard  and  Inman  (1981)  demonstrated  that  equation  (83) 
could  be  more  accurately  expressed  as 

tane  =  (tan^)(S)  ,  (86) 

where 

S  =  skewness  of  the  oscillatory  velocity  distribution, 
defined  as 
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S  = 


Y    /  U-^(t)dt 

t 


t+T 

/ 
t 


1  ^"^1      3 
J    /  |U(t)|  dt 


Since  they  state  that  equation  (83)  can  be  shown  to  be  approxi- 
mately equal  to  equation  (86),  when  "c"  is  given  by  equation  (85),  it 
is  perhaps  wise  to  use  equation  (83),  in  view  of  the  almost  total  lack 
of  knowledge  regarding  the  time  variation  of  the  bottom  oscillatory 
velocity  field. 

Equating  equations  (84)  and  (85)  yields 


m-off 


U 
m-on 


1   dh 


tand)  dx 


1  +-L_dh 
tan*  dx 


1/3 


(87) 


Since  the  beach  slope  must  always  be  less  than  the  angle  of  repose, 


n    1   dh   ^ 
^  ^  tand)  dx  ^  ^  ' 


then  the  right-hand  side  of  equation  (87)  can  be  expanded  in  terms  of 
two  binomial  series,  as  shown  below. 


m-off 


U 


m-on 


1       dh  _       1         ,  dh 

Stand)  dx  "  „+     2,    Mx 
^  9tan  d) 


)  + 


i+..^4^+    1 


_  ^ ,dh , 

3tan(i)  dx  "^  „,  2,  ^dx' 
9tan  i> 


Using  just  the  first  two  terms  of  the  expansions,  it  is  seen  that 


U     ,   1  ._i_dh 
m-off  ^    3tanij)  dx 

!Von  I  "  1+^^^ 
3tand)  dx 
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To  test  this  approximation,  the  probable  steepest  beach  slope  of 
dh/dx  =  0.1  was  chosen  along  with  an  angle  of  repose  of  29°,  taken  from 
page  463  of  Henderson  (1966),  which  shows  <j)  :^  29°  over  the  beach  sedi- 
ment range  of  .25  mm  to  2.0  mm.  Using  these  values  in  equation  (88) 
produces 


m-off 


m-on 


0.8866 


while  the  exact  equation  (87)  gives 

U„ 


m-off 


U 
m-on 


=  0.8855  . 


Thus  it  can  be  concluded  that  equation  (88)  is  an  accurate  approxima- 
tion of  equation  (87) . 

Experimental  measurements  of  the  maximum  bottom  velocity  asymmetry 
carried  out  by  Hamada  (1951)  indicated  that  U  r.f^  varied  between 

0.78  and  0.86,  which  is  in  good  agreement  with  the  above  values. 

It  is  now  desirable  to  find  an  expression  for  the  near-bottom 
maximum  velocity  asymmetry. 

Using  the  horizontal  particle  velocity  under  the  wave  crests  as 
the  single  most  important  feature,  LeMehaute  et  al .  (1958)  compared 
12  shallow  water  wave  theories.  Among  their  conclusion  was  the  sur- 
prising fact  that  simple  1st  order  linear  (Airy)  theory  provides  a 
reasonable  approximation  of  the  bottom  velocity  which  the  higher  order 
Stokes  theories  never  approach.  Wave  theories  providing  better 
approximations  than  linear  theory  include  Solitary  and  Cnoidal  wave 
theory.  Since  it  is  desirable  to  express  the  bottom  maximum  velocities 
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as  a  function  of  the  depth,  h,  an  analytical  form,  such  as  linear 
theory,  is  preferred  over  the  numerical  Cnoidal  theory  or  Solitary 
theory,  which  give  only  onshore  horizontal  bottom  velocity. 

Goda  (1964)  observed  that  maximum  horizontal  velocities  under  the 
wave  crest  for  nearshore,  unbroken  waves  were  slightly  larger  than 
predicted  by  linear  theory. 

By  superimposing  a  small  onshore-directed  current,  u",  on  the 
oscillatory  wave  bottom-velocity  from  linear  theory,  given  as 


U(t)  =  U^  cos  tut  , 


where 

ttH 


m   T  sinh(kh)  ' 

the  onshore  maximum  velocity  magnitude  can  be  expressed  as  U  +  "u ; 

m 

the  maximum  offshore  velocity  magnitude  becomes  U  -  u",  which  gives 

m  ^ 


m-off 

U 
m-on  ' 


U„  -  u 
m 

1  -  ^ 
m 

U     +  u 
m 

''I 

m 

(89) 


where  u  <<  U  . 
m 


Equating  equation  (88)  and  equation  (89)  gives 


1   dh   1  -  ^ 

Stand)  dx       m 


1+   1   dh 


Stancj)  dx   1+7]- 
m 


from  which  it  is  immediately  seen  that 

u_  ^   1   dh 

U^   3tan(})  dx 
m      ^ 


:9o) 
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Assuming  that  u/U^  can  be  written  as  a  function  of  depth,  h,  only, 
equation  (90)  becomes 


U 


dh  =  Stand)  dx  . 


(91) 


With  the  assumption  that  only  small  amounts  of  total  wave  energy  are 
lost  by  movement  of  sediment  to  and  fro,  the  local  wave  height  may  be 
expressed,  approximately,  by  the  linear  shoaling  relationship 


H  -  H. 


0 


2cosh  kh 


il/2 


sinh  2kh  +  2kh 


(92) 


for  waves  propagating  normal  to  the  beach. 

As  pointed  out  by  Eagleson  and  Dean  (1966),  experiments  indicate 
that  equation  (92)  predicts  smaller  wave  heights  than  measured  on 
sloping  bottoms  in  wave  tank  experiments.  However,  Guza  and  Thornton 
(1980)  measured  wave  spectra  on  a  line  extending  from  a  10  meter  depth 
to  inside  the  surf  zone  at  Torry  Pines  Beach,  California,  and  compared 
shallow  water  spectra  with  those  calculated  using  linear  shoaling 
theory  applied  to  the  10  meter  depth  spectrum.  Comparisons  were  very 
favorable,  particularly  in  the  dominant  frequency  band,  which  indi- 
cates that  equation  (92)  is  perhaps  a  reasonable  assumption  to  make. 

Substitution  of  equation  (92)  into  equation  (91),  making  use  of 
the  relationships 


yields 


ttH  j     1        qT 

'm  ^  Tsinh   kh     ^"^     L  ^  4r 


gHj 


2L       sinh  2kh  +  2kh 


tanh  kh   , 


1/2 


(93) 
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For  shallow  water  waves,  sinh  2kh  ^^   2kh.  Noting  that  k  =  2irl,   equation 
(93)  becomes 

Eagleson  and  Dean  (1966)  also  point  out  that  wave  celerity  is  ade- 
quately predisted  by  linear  theory  for  shoaling  waves  up  to  breaking. 
Thus  for  shallow  water, 

C  =  (gh)l/2  =^    or     I  -   (gh)^/^^  ^ 

and  equation  (94)  then  becomes 


3/4^  yl/2 


"™  -  ^77?^  •  ''=' 


Placing  equation  (95)  into  equation  (91)  yields 


3/4   1/2 
9   "o^    dh 


— 1  /p  _ T7T  ^  3tan(|)  dx  .  (96) 

477^/"^  u  h^^^ 


It  is  now  necessary  to  determine  a  form  for  u"  as  a  function  of  h, 
H^,  and  T.  For  the  wave  velocity  asymmetry  to  increase  as  depth  de- 
creases, it  is  necessary  for  the  ratio  U/U  to  increase.  Since 
^1^  -^  1/h  '  ,  it  can  be  stated  in  general  that  u" -v  1/h"  where  n  >  3/4. 
Thus, 

u    h3/4     1 


Um    h"    h"-3/4 
By  this  logic,  it  can  be  said  that 
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.(h^,t; 


where  a  is  an  arbitrary  function  and  n  >  3/4.  Equation  (96)  can  now 
be  written  and  integrated  as 

3/4„  t1/2  u 

wp  f  h""^^^dh  =  Stan*  f  dx  , 

where  x^.  is  a  reference  location,  yielding  a  power  equation  for  the 
equilibrium  beach  profile  seaward  of  the  breaking  point  given  as 


12n  +  3)TT-^/^tan 


g3/4,  ,1/2 


H„j; 


4n+l 


(x-x.)4"^l.        (97; 


To  find  an  upper  limit  for  the  exponent,  n,  assume  u"  is  equal  to  the 
mass  transport  velocity  obtained  from  2nd  order  Stokes  theory,  given 
by  Eagleson  and  Dean  (1966)  at  the  bottom  as 


u 


'lM^2    C 

2sinh^  kh 


(98) 


Applying  the  linear  wave  shoaling  equation  (92)  and  the  shallow  water 
wave  length  expression  given  before,  this  expression  becomes 


_   gTH/   a(T,H^) 


32TTh^     h^ 


(99) 


Unfortunately,  2nd  order  Stokes  theory  gives  very   poor  bottom  velocity 
results  in  shallow  water,  as  stated  before,  with  maximum  velocities 
being  substantially  larger  than  actual  measurements  indicate.  For  this 
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reason,  as  depth  decreases,  the  theoretical  mass  transport  velocity  at 
the  bottom  increases  too  rapidly;  as  a  consequence,  the  exponent  n  =  2 
in  equation  (99)  represents  an  upper  limit. 

Writing  equation  (97)  in  shorter  notation  as 

h  -  A(x')'",  where  m  =  ^^\  ^  and  3/4  <  n  <  2  , 

it  can  be  seen  that  m  varies  between  4/9  and  1. 

It  has  been  decided  to  choose  m  =  2/3,  which  lies  almost  exactly 
between  4/9  and  1  and  corresponds  to  n  =  5/4  for  the  following  reasons: 

1)  This  gives  the  uniform  onshore  drift  current  as 


,5/4  ■ 

which  seems  very  reasonable  in  view  of  the  analysis  done  using 
2nd  order  Stokes  theory. 

2)  When  Hughes  (1978)  applied  the  curve  h  =  Ax"^  to  464  actual 
beach  profiles  from  the  east  coast  of  Florida  extending  to 
depths  of  about  20  feet,  he  found  that  m  =  2/3  provided  the 
"best  fit"  in  a  majority  of  the  cases.  Greater  error  was 
found  when  the  curve  was  fitted  to  depths  of  about  32  feet. 

3)  Dean  (1977)  was  able  to  derive  an  equation  of  the  form 

2/3 
h  =  Ax    for  the  beach  profile  inside  the  surf  zone  based  on 

the  concept  of  uniform  dissipation,  per  unit  water  volume,  of 

incoming  wave  energy.  He  also  employed  linear  wave  theory  in 

his  derivation.  It  is  reasonable  to  believe  that  there  exists 

a  type  of  "profile  similarity"  between  the  immediate  nearshore 

region  and  the  surf  zone,  since  many  of  the  surf  zone 
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characteristics  can  be  expressed  in  terms  of  Battjes'  (1974) 
surf  similarity  parameter. 

4)  As  pointed  out  in  Chapter  IV,  the  equilibrium  profile  given 

2/3 
by  h  =  Ax  ^   is  preserved  in  the  model  law  derived  for  these 

experiments. 

5)  Bruun  (1954)  also  proposed  an  equation  of  the  form  h  =  Ax^'^"^ 
for  the  region  just  seaward  of  the  breaker  line  and  found  good 
correspondence  when  applied  to  actual  beach  profiles. 

So  finally,  the  proposed  equation  for  the  bottom  profile  just 
seaward  of  the  breaker  line  is  given  as 


1/2 

L,  _  I  Utt   tan*    /n  X' 

'   -  ',3/4,  ^1/2'  •  ^'(^o'T. 
0 


i2/3 


.2/3 


(100) 


At  this  point,  little  can  be  done  about  determining  the  function 
for  a(H^,T).  Dimensionally,  a  has  the  units  of 


[=] 


9/4 


but  dimensional  analysis  cannot  be  employed  because  a  may  also  contain 
the  constants  of  gravity,  fluid  density,  etc. 
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